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ABSTRACT 

An  investigation  of  the  interaction  effects  of  a  lower  heated  tube  on  pool  boiling  of 
pure  R-124  from  an  upper  horizontal  tube  was  conducted  at  a  saturation  temperature  of 
2.2  °C.  The  test  tubes  used  were:  (1)  smooth  tubes  and  (2)  deformed  surface  (TURBO-B) 
enhanced  tubes.  The  effects  of  tube  spacing/configuration  and  lower  tube  heat  flux  on  the 
heat  transfer  performance  of  the  upper  tube  were  investigated. 

For  both  tube  arrays,  the  enhancing  effect  of  bubbles  from  a  lower  tube  was 
dramatic.  This  enhancement  increased  as  lower  tube  heat  flux  increased.  However,  when 
upper  tube  heat  fluxes  were  greater  than  20  kW/m2,  all  enhancement  disappeared. 

For  a  smooth  tube  array  in  natural  convection,  the  effect  of  a  lower  heated  tube  on 
the  heat  transfer  from  an  upper  tube  was  small.  In  nucleate  boiling,  a  P/D  of  1.8  gave  the 
best  upper  tube  heat  transfer  performance  and  a  vigorously  nucleating  lower  tube 
eliminated  upper  tube  hysteresis.  With  the  lower  tube  unheated  and  an  upper  tube  heat 
flux  of  greater  than  3  kW/m2,  the  performance  using  R-124  was  generally  better  than  for 
R-114.  With  a  nucleating  lower  tube  (at  10  kW/m2),  again  the  performance  of  R-124  was 
better,  but  only  for  upper  tube  heat  fluxes  of  greater  than  40  kW/m2. 

For  a  TURBO-B  tube  array,  a  30  degree  offset  of  the  upper  tube  reduced  the  upper 
tube  heat  transfer  performance  (compared  to  the  in-line  configurations).  This  may 
indicate  bubbles  depart  TURBO-B  tubes  differently  than  smooth  tubes. 
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NOMENCLATURE 

This  nomenclature  reflects  that  which  is  found  within 
programs  SETUPJY  (Appendix  A)  and  DRPJY  (Appendix  B) .  Note 
that  the  upper  and  lower  test  tube  dimensions  are  the  same 
since  like  test  tubes  are  used. 

SYMBOL   UNITS        NAME/DESCRIPTION 

Area 

Tube  active  boiling  surface  area 

Tube  cross-sectional  area 

Specific  heat  of  upper  tube 

Specific  heat  of  lower  tube 

Diameter 

Tube  outside  diameter 

Tube  inside  diameter 
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Upper  tube  heat  flux 
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°c 
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°c 
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°c 
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AV 
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8b 

K 
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K 
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Upper  tube  average  wall  temperature  at  a 
thermocouple  location 

Lower  tube  average  wall  temperature  at  a 
thermocouple  location 

Critical  temperature 

Upper  tube  film  temperature 
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Upper  tube  outer  wall  temperature 
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Au      m2/s        Kinematic  viscosity  associated  with  the 
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p       kg/m3       Density  associated  with  the  upper  tube 

Ap      kg/m3       Density  associated  with  the  lower  tube 
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I .   INTRODUCTION 

The  United  States  Navy  intends  to  eliminate  its  dependence 
on  chlorof luorocarbons  (CFC's)  as  soon  as  feasible. 
Currently,  the  Navy  is  using  CFC-114  (Refrigerant  or  R-114)  in 
most  of  it's  chilled  water  air  conditioning  (AC)  plants. 
These  shipboard  AC  plants  provide  the  chilled  cooling  water  to 
cool  the  vital  weapons  systems  and  electronic  components 
besides  providing  habitability  comfort  for  the  crew. 

The  Navy  has  identified  a  hydrochlorof luorocarbon,  HCFC- 
124  or  R-124,  as  a  possible  mid-term  replacement  for  R-114  in 
naval  shipboard  AC  systems.  With  an  ozone  depletion  potential 
(ODP)  of  only  2%  compared  to  70%  for  R-114,  R-124  is  a  much 
more  tolerable  fluid  to  use  until  an  acceptable  alternative 
refrigerant  can  be  found  that  has  zero  ODP. 

Although  both  refrigerants  operate  near  the  same 
evaporating  temperature  (2.2  °C)  ,  the  characteristics  of  R-124 
under  similar  conditions  as  R-114  are  quite  different.  For 
example,  the  refrigeration  components  (i.e.,  heat  exchangers, 
piping,  etc.)  must  be  able  to  withstand:  (1)  the  higher 
operating  pressures  of  R-124  (11  psig  at  2.2  °C  compared  to  0 
psig  for  R-114)  and  (2)  the  higher  vapor  pressure  after  the 
system  is  secured  (in  a  main  engineering  space  at  100  °F,  R- 
124  would  reach  66  psig  compared  to  31  psig  for  R-114)  .  Given 
the  system  components  can  operate  safely  at  these  higher 


pressures  of  R-124,  then  the  hardware  modifications  to  the 
system  should  be  limited  to: 


1.  Modifying  the  compressor  to  maintain  rated  capacity  and 
use  a  miscible  alkylbenzene  oil  for  lubrication. 

2.  Reset/replace  the  system  safety  devices  to  account  for 
the  higher  pressures  of  R-124. 

3.  Replace  the  system  seals  with  neoprene  seals,  which 
resist  deterioration  when  in  contact  with  R-124,  to 
maintain  leak  tight  integrity. 


Besides  hardware  modifications,  another  main  concern  is 
how  does  the  heat  transfer  characteristics  of  R-124 
compare  to  R-114.   Table  1,  from  Bertsch  [Ref.  1],  provides 

TABLE  1.   R-114/R-124  PHYSICAL  PROPERTIES  AT  2 . 2  °C 


Property 

R-114 

R-124 

Psat  (kPa) 

96 

177 

p(liq)  (kg/m2) 

1526 

1429 

p (vapor)  (kg/m2) 

7.5 

11.3 

Cp(liq)  (J/kg-K) 

932 

1060 

hfq  (kJ/kg) 

134.7 

158.9 

k(liq)  (W/m-K) 

0.0701 

0.0744 

u(liq)  (g/m-s) 

0.4449 

0.3375 

o  (N/m) 

0.0136 

0.0129 

the  physical  properties  of  R-114  and  R-124  at  a  saturation 
temperature  of  2.2  °C.  As  seen  from  Table  1,  R-124  has 
higher   latent  heat  of  vaporization,   vapor  density  and 


thermal  conductivity  than  R-114  at  2.2  °C  and  a  considerably 
lower  liquid  viscosity.  Therefore,  it  is  reasonable  for  one 
to  assume  that  R-124  would  have  better  heat  transfer 
performance  at  2.2  °C  than  R-114.  Bertsch  [Ref.  1]  evaluated 
and  compared  the  heat  transfer  performance  of  R-114  and  R-124 
for  a  single  horizontal  smooth  tube  and  various  enhanced  tubes 
using  pure  refrigerant  and  refrigerant/oil  mixtures.  One  of 
his  conclusions  was: 


For  smooth  and  finned  tubes,  significant  improvements  in 
heat  transfer  (>50%)  were  obtained  in  pure  R-124  compared 
with  pure  R-114  at  the  same  saturation  temperature.  This 
was  attributed  to  the  higher  saturation  pressure  of  R-124 
activating  more  nucleation  sites  and  the  higher  latent 
heat  of  vaporization,  vapor  density  and  thermal 
conductivity  of  R-124  over  R-114. 


Perry  [Ref.  2]  investigated  the  effect  of  oil,  pressurization 
and  subcooling  on  the  onset  of  nucleate  boiling  of  R-124  from 
a  single  horizontal  tube.  When  Perry  [Ref.  2]  compared  his 
results  to  Bertsch  [Ref.  1]  for  pure  R-124,  he  found  a  larger 
incipient  wall  superheat  was  being  obtained,  but  almost 
identical  data  trends  were  observed  at  high  heat  fluxes  and 
decreasing  heat  fluxes.  He  attributed  the  larger  incipience 
to  his  degassing  procedure  of  the  refrigerent  system  after 
reassembly. 

Now,  what  if  multiple  tubes  are  arranged  in  a  bundle 
configuration?  How  do  the  lower  tubes  influence  the  heat 
transfer  performance  of  the  upper  tubes?  Lake  [Ref.  3]  in 
search  of  an  answer  to  these  questions,  used  a  configuration 


of  two  horizontal  tubes  (smooth  and  enhanced  tubes)  in 
vertical  alignment,  submerged  in  a  pool  of  R-114,  separated  by 
a  pitch-to-diameter  (P/D)  ratio.  After  using  various 
combinations  of  P/D  ratios  and  lower  tube  heat  fluxes,  Lake 
[Ref.  3]  found  that  a  nucleating  lower  tube  significantly 
enhanced  the  heat  transfer  performance  of  an  upper  tube, 
particularly  at  low  heat  fluxes.  Also,  a  P/D  ratio  of  1.8 
gave  the  best  heat  transfer  performance  for  smooth  tubes  at 
high  heat  fluxes  (no  P/D  effect  was  found  using  the  enhanced 
tubes) . 

To  fully  understand  the  nucleate  pool  boiling 
characteristics  of  R-124  and  to  directly  compare  the  heat 
transfer  performance  of  R-124  to  R-114,  multiple  tube  pool 
boiling  data  for  R-124  is  needed.  Therefore,  the  following 
thesis  objectives  were  established: 


1.  Collect  natural  convection  and  pool  boiling  heat  transfer 
data  of  pure  R-124  using  a  single  smooth  and  enhanced 
tube  and  check  the  repeatability  with  the  existing  R-124 
data  of  Bertsch  [Ref.  1]  and  Perry  [Ref.  2] . 

2.  Collect  natural  convection  and  pool  boiling  heat  transfer 
data  of  R-124  over  a  range  of  increasing  and  decreasing 
heat  fluxes  using  two  tubes  while  varying  the  spacing 
between  the  tubes  and  the  heat  flux  on  the  lower  tube. 

3.  Compare  the  R-124  data  to  the  R-114  data  obtained  by  Lake 
[Ref.  3]  . 


II.   LITERATURE  SURVEY 

A.   POOL  BOILING 

Pool  boiling  occurs  when  a  heated  surface  is  submerged 
below  the  free  surface  of  a  liquid  pool,  at  a  temperature 
sufficiently  above  the  local  saturation  temperature  to  cause 
bubbles  to  form  from  the  surface.  Figure  2.1,  from  Bar-Cohen 
[Ref.  4],  is  a  pool  boiling  curve  for  a  highly  wetting  liquid 
with  several  regimes  displayed.  Only  the  natural  convection 
(or  single  phase  convection) ,  boiling  incipience  and  nucleate 
boiling  regimes  are  discussed. 

1 .  Natural  Convection 

The  regime  (a-b)  in  Figure  2.1,  is  where  heat  transfer 
occurs  because  of  the  liquid  flow  caused  by  buoyancy  forces 
due  to  density  differences  in  the  liquid  pool.  In  this 
regime,  the  heat  flux  from  the  heated  surface  is  proportional 
to  some  power  of  wall  superheat  (i.e.,  q"  «  ATW1,25  or  ATW1,33). 
This  proportionality  remains  valid  until  the  wall  superheat 
reaches  a  value  where  vapor  bubbles  first  form  on  the  heated 
surface. 

2 .  Boiling  Incipience 

The  onset  of  nucleate  boiling  (ONB)  or  incipient 
boiling  is  a  random  phenomenon  (You  et  al.  [Ref.  5])  which 
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Figure  2.1   Pool  Boiling  Curve 


is  influenced  by  factors  such  as  surface  condition  of  the 

heated  surface,  liquid  properties  and  dissolved  gases  in  the 

liquid  pool.   Bar-Cohen  and  Simon  [Ref.  6]  and  Tong  et  al. 

[Ref.  7]  investigated  boiling  incipience  of  highly-wetting 

liquids  using  refrigerants  and  dielectrics.   They  found  that 

these  liquids  effectively  flood  the  micro-cavities  on  a  tube 

surface,   causing   fewer   active   embryos   (micro-cavities 

containing  trapped  vapor) .   Ultimately,  these  highly-wetting 

liquids  require  larger  wall  superheats  to  initialize  bubble 

growth  in  the  flooded  micro-cavities.   As  a  result  of  the 

larger  wall  superheats  (up  to  100  °C,  Bar-Cohen  and  Simon 

[Ref.  6]),  nucleation  at  (b)  ,  in  Figure  2.1,  occurs  in  an 

explosive  manner.   This  explosive  boiling  from  the  active 

micro-cavities   causes   the   larger   inactive   cavities   to 

activate,  rapidly  decreasing  the  required  wall  superheat  to 

sustain  boiling,  noted  as  ATex  in  Figure  2.1.   This  drift  from 

(b)  to  (c)  only  occurs  during  increasing  heat  flux  and  is 

referred  to  as  thermal  overshoot  or  superheat  excursion. 

3 .   Nucleate  Boiling 

The  regime  (c-d)  in  Figure  2.1  is  the  fully  developed 
nucleate  boiling  region.  At  (c)  ,  small  isolated  bubbles  are 
leaving  the  heated  surface  at  low  heat  fluxes.  As  the  heat 
flux  is  increased,  these  small  bubbles  become  more  numerous  as 
more  active  sites  are  created.  At  high  heat  fluxes,  these 
bubbles  begin  to  coalesce,  forming  slugs  of  vapor,  which  grow 


and  depart  from  the  surface.  According  to  Thome  [Ref.  8],  in 
this  regime,  heat  removal  is  accomplished  by  sensible  heat  in 
the  form  of  superheated  liquid  or  as  latent  heat  in  vapor 
bubbles.  Also,  the  principal  mechanisms  of  heat  transfer  are: 
(1)  liquid-phase  convection,  enhanced  by  the  stripping  of  the 
thermal  boundary  layer  by  the  departing  vapor  bubbles,  and  (2) 
evaporation  of  the  thin  liquid  micro-layer  trapped  between  a 
growing  vapor  bubble  and  the  heated  surface. 

B.   SINGLE  TUBE  HEAT  TRANSFER  EXPERIMENTS  WITH  R-124  AT  NAVAL 
POSTGRADUATE  SCHOOL 

Bertsch  [Ref.  1]  and  Perry  [Ref.  2]  independently  studied 

the  heat  transfer  performance  of  single  smooth  and  enhanced 

tubes  submerged  in  a  liquid  pool  of  R-124  at  2.2  °C.   Perry 

[Ref.   2]   determined   that   in   natural   convection   the 

experimentally  based  correlation,  equation  (1),  developed  by 

Churchill  and  Chu  [Ref.  9]  fairly  well  predicted  the  heat 

transfer  performance  of  a  smooth  tube. 


Nu  = 


,    12 

0.6+0.387/ 


(1+  (0.559/Pr)9/16)16/9 


(1) 


where   Nu  =  Nusselt  number  =  hDo/k 

Ra  =  Rayleigh  number  =  (g-p-{Twall  -  Tsat }  -Do3)  / (u -a, 


During  nucleate  boiling,  most  correlations  available  to 
predict  the  heat  transfer  performance  of  a  smooth  tube  are 
expressed  as  heat  flux  as  a  function  of  tube  wall  superheat. 
Bertsch  [Ref.  1]  used  a  best  fit  correlation  through  his 
experimental  nucleate  boiling  data  to  express  the  heat  flux 
from  a  smooth  single  tube  in  pure  R-124  as: 

q"  =  (17.0)ATW3   (W/m2)        (2) 

In  the  above  equation,  the  wall  superheat  ATW  is  expressed  in 
degrees  Celsius. 

For  a  single  enhanced  tube,  Bertsch  [Ref.  1]  concluded 
that  nucleate  pool  boiling  enhancement  (i.e.,  improvements  in 
heat  transfer  compared  to  a  single  smooth  tube)  were  as  much 
as  6.2  times  for  re-entrant  cavity  surfaces  (TURBO-B  and  HIGH 
FLUX).  Thome  [Ref.  10]  attributes  these  enhancements  to  the 
numerous  re-entrant  cavities  that  trap  vapor,  which  create  a 
large  boiling  site  density,  reducing  the  wall  superheat 
required  to  sustain  nucleation. 

C.   SMALL  TUBE  ARRAY  HEAT  TRANSFER  PERFORMANCE 

The  effect  of  a  lower  tube  on  the  heat  transfer 
performance  of  an  upper  tube  is  very  important  when  designing 
flooded  evaporators.  It  is  essential  to  understand  the 
mechanisms  that  affect  the  heat  transfer  performance  in  small 


tube  arrays  in  natural  convection  and  nucleate  boiling  for 
both  smooth  and  enhanced  surfaces. 
1 .   Natural  Convection 

In  natural  convection,  the  interaction  effects  of  a 
lower  heated  tube  on  the  heat  transfer  performance  of  an  upper 
tube  are  the  result  of  two  offsetting  mechanisms.  With  a 
lower  tube  heated,  a  rising  heated  buoyant  plume  impinges  on 
and  flows  around  an  upper  tube,  which  subjects  the  upper  tube 
to  local  forced  convection  effects.  This  forced  convection 
increases  the  local  fluid  velocity  around  the  upper  tube  and 
raises  the  upper  tube  heat  transfer  coefficient. 
Additionally,  the  heated  plume  from  the  lower  tube  increases 
the  local  ambient  temperature  of  the  fluid  impinging  on  the 
upper  tube.  This  reduces  the  natural  convection  buoyancy 
driving  force  for  this  tube,  decreasing  the  upper  tube  heat 
transfer  coefficient.  Therefore,  the  presence  of  a  lower 
heated  tube  has  opposing  influences  on  the  heat  transfer 
performance  of  an  upper  tube.  Which  of  these  influences  is 
dominant  depends  on  the  tube  spacing,  local  heat  flux  and  the 
thermophysical  properties  of  the  fluid. 

Sparrow  and  Niethammer  [Ref .  11]  studied  the  effect  of 
vertical  separation  between  two  horizontal  cylinders  on  the 
heat  transfer  performance  of  an  upper  cylinder  in  natural 
convection  using  air.   It  was  discovered  that  heat  transfer 
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enhancement1  increased  as  the  distance  between  the  tubes 
increased  with  a  maximum  enhancement  of  approximately  1.3  for 
7  <  P/D  <  9.  Also,  for  a  P/D  <  3,  the  enhancement  decreases 
below  1.0.  Al-Alusi  and  Bushnell  [Ref.  12]  and  Tokura  et  al. 
[Ref .  13]  also  working  with  a  small  array  of  horizontal  smooth 
cylinders  in  air,  found  generally  the  same  results.  As  the 
distance  between  the  cylinders  increased,  the  heat  transfer 
enhancement  increased,  reaching  a  maximum,  then  returning  to 
a  lower  value  of  enhancement. 

Sparrow  and  Boessneck  [Ref.  14]  investigated  the 
effect  of  transverse  misalignment  of  a  lower  cylinder  on  the 
heat  transfer  performance  of  an  upper  cylinder,  for  2  £  P/D  £ 
9,  using  air.  For  the  perfectly  aligned  upper  and  lower 
cylinders,  the  heat  transfer  enhancement  of  the  upper  cylinder 
increased  as  the  P/D  ratio  increased.  Again,  the  largest 
enhancement  was  found  with  a  P/D  of  9  and  the  P/D  of  2 
produced  the  lowest  enhancment.  As  misalignment  was 
increased,  the  enhancement  began  to  merge  to  a  common  value  of 
about  1.03  at  sufficiently  large  offsets  (about  2.5  cylinder 
diameters) .  Sparrow  and  Boessneck  concluded  that  this 
behavior  is  caused  by  the  rising  plume  from  the  lower  cylinder 
presenting  a  pre-heated  quasi-forced-convection  flow  to  the 


Enhancement  is  defined  for  a  small  tube  array  as  the 
upper  cylinder  heat  transfer  coefficient  when  both  cylinders 
are  heated  divided  by  the  upper  cylinder  heat  transfer 
coefficient  when  the  lower  tube  is  unheated. 
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upper  cylinder.  At  smaller  P/D  ratios,  the  increased  local 
ambient  temperatures  of  the  fluid  impinging  on  the  upper 
cylinder  dominates  and  decreases  the  upper  cylinder  heat 
transfer  performance.  At  larger  P/D  ratios,  the  forced 
convection  effects  are  more  dominant  and  increase  the  heat 
transfer  from  the  upper  cylinder. 

Inagaki  and  Komori  [Ref.  15]  studied  natural 
convection  heat  transfer  from  two  vertically  aligned 
horizontal  cylinders  in  water.  As  discovered  with  air,  a 
decrease  in  enhancement  which  falls  below  1.0,  is  found  for 
P/D  <  1.0.  A  maximum  enhancement  of  1.25  occurs  between  2.5 
£  P/D  £  3.0  and  at  larger  P/D  ratios,  the  enhancement 
decreases  to  1.0.  Inagaki  and  Komori  attribute  this  behavior 
to  the  same  two  opposing  influences  discussed  previously. 
2 .   Nucleate  Boiling 

For  small  tube  arrays  in  nucleate  boiling,  the  most 
significant  mechanism  of  enhancement  within  the  array,  is  the 
well  documented  bundle  effect2.  Cornwell  and  associates 
[Refs.  16-19]  studied  fluid  flow  in  tube  bundles  and  found 
significiant  enhancements  in  heat  transfer  in  the  upper 
portions  of  the  bundles,  concluding  that  forced  convection  and 
bubbly  flow  from  the  lower  tubes,  creating  sliding  bubbles 


2Bundle  effect  is  the  enhancement  of  an  upper  tube  due  to 
the  activation  (i.e.,  the  nucleation  of  bubbles)  of  lower 
tubes  within  the  bundle. 
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around  the  upper  tube,  cause  significant  enhancement  at  low 
heat  flux  conditions. 

Fujita  et  al.  [Ref.  20]  conducted  tests  with  an  array 
of  three  horizontal  smooth  tubes  in  R-113  at  0.1  MPa.  For 
heat  fluxes  <  20  kW/m2,  the  bundle  effect  enhanced  the  upper 
tube  heat  transfer  performance.  Fujita  et  al .  [Ref.  20]  also 
used  an  array  of  two  smooth  tubes,  where  the  lower  tube  heat 
flux  was  held  constant  and  the  upper  tube  heat  flux  was 
decreased  uniformly.  The  upper  tube  enhancement  increased  as 
the  lower  tube  heat  flux  was  raised.  Lake  [Ref.  3] 
investigated  the  effect  of  P/D  ratio  and  lower  tube  heat  flux 
on  the  heat  transfer  performance  of  an  upper  tube,  using  an 
array  of  two  smooth  tubes  in  R-114  at  1  atmosphere.  For  these 
tests,  the  bottom  tube  heat  flux  was  held  constant  and  the 
upper  tube  heat  flux  was  raised  and  lowered  steadily.  Lake 
found  that  for  an  upper  tube  heat  flux  of  <  20  kW/m2,  the 
enhancement  of  the  upper  tube  was  due  to  a  bundle  effect  and 
the  amount  of  enhancement  increased  with  increased  lower  tube 
heat  flux,  which  is  consistent  with  Fujita  et  al .  [Ref.  20]. 
Also,  Lake  found  that  for  upper  tube  heat  fluxes  of  >  20 
kW/m2,  all  bundle  effects  faded  and  that  the  P/D  ratio  of  1.8 
provided  the  best  upper  tube  heat  transfer  performance. 

For  a  small  array  of  enhanced  tubes,  several  authors 
report  findings  similar  to  those  found  for  the  smooth  tube 
arrays.  Stephan  and  Mitrovic  [Ref.  21],  using  a  small  array 
of  three  horizontal  T-shaped  finned  tubes  in  R-114  at  0.2  MPa 
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and  R-12  at  0.15  MPa,  found  bundle  effects  present  at  heat 
fluxes  <  20  kW/m2,  with  these  effects  disappearing  as  the  heat 
flux  was  raised.  Hahne  and  associates  [Refs.  22-24],  using 
two  different  fin  arrays  in  R-ll  at  0.1  MPa,  also  found 
significant  bundle  effects  for  heat  fluxes  of  <  20  kW/m2, 
which  subside  at  higher  heat  fluxes.  Hahne  et  al.  [Ref.  24] 
also  investigated  the  influence  of  P/D  ratio  on  upper  tube 
heat  transfer.  For  heat  fluxes  <  20  kW/m2,  the  maximum  P/D 
ratio  tested  (P/D  of  3.0)  produced  the  best  heat  transfer 
performance  for  both  finned  arrays.  Above  20  kW/m2,  neither 
P/D  ratio  nor  tube  fin  pitch  affected  the  upper  tube  heat 
transfer.  Lake  [?ef.  3],  in  conjunction  with  his  testing  of 
smooth  tube  performance  above,  found  no  effect  of  P/D  ratio  on 
HIGH  FLUX  tube  (porous  coated  tube)  performance  at  any  heat 
flux  and  that  there  was  no  effect  from  a  lower  heated  tube  on 
the  heat  transfer  performance  of  a  nucleating  upper  HIGH  FLUX 
tube.  Fujita  et  al.  [Ref.  20],  like  Lake,  also  tested  porous 
coated  tubes  using  the  same  testing  conditions  as  for  the 
smooth  tubes  and  he  concluded  that  the  lower  tube  heat  flux 
setting  had  no  influence  on  the  heat  transfer  coefficient  of 
the  upper  tube,  which  agrees  with  the  conclusions  of  Lake 
[Ref.  3] . 
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III.    EXPERIMENTAL  APPARATUS 

A.   SYSTEM  DESCRIPTION 

A  schematic  diagram  of  the  experimental  apparatus  is  shown 
in  Figure  3.1.  Details  of  the  apparatus  design,  construction 
and  calibration  are  discussed  by  Karasabun  [Ref.  25],  Reilley 
[Ref.  26]  and  Sugiyama  [Ref.  27].  The  apparatus  is  very 
similar  to  the  one  used  by  Lake  [Ref.  3]  to  conduct  pool 
boiling  of  pure  R-114  with  a  two  tube  configuration.  The  only 
major  change  being  the  replacement  of  the  Pyrex-glass  Tees 
used  as  the  condenser  and  evaporator.  Bertsch  [Ref.  1] 
replaced  the  Pyrex-glass  Tees  with  fiberglass  reinforced 
Pyrex-glass  Tees  enabling  the  system  to  withstand  the 
increased  operating  pressure  of  R-124. 

As  shown  in  Figure  3.1,  the  present  system  contains  the 
following  components:  (1)  a  reinforced  Pyrex-glass  Tee,  used 
as  the  evaporator  to  boil  the  refrigerant;  (2)  a  reinforced 
Pyrex-glass  Tee,  used  as  the  condenser  to  condense  the 
refrigerant  vapor;  (3)  a  liquid  refrigerant  storage 
reservoir;  (4)  an  oil  storage  and  mixing  subsystem  (not 
used);  (5)  a  vacuum  pump;  (6)  an  ethylene-glycol/water 
mixture  cooling  subsystem;  (7)  a  condenser  pump; 
(8)  aluminum  framework  with  plexiglass  siding  which  houses 
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Figure  3 . 1   Schematic  of  Experimental  Apparatus 
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components  (1)  through  (4);   and  (9)  a  data  acquisition  and 
instrumentation  system  (not  shown  in  Figure  3.1). 

The  system  operates  in  a  closed  loop  with  vapor  being 
generated  by  the  heated  tubes  submerged  in  the  R-124  liquid 
pool  in  the  evaporator.  The  vapor  travels  to  the  condenser, 
via  a  two  inch  diameter  aluminum  L-shaped  pipe,  where  it's 
condensed  on  a  copper  coil  cooled  by  the  chilled  ethylene- 
glycol/water  mixture.  The  condensed  liquid  refrigerant  is 
gravity  fed  back  to  the  evaporator  via  the  return  line. 

B.   BOILING  TEST  SECTION 
1 .   Evaporator 

The  evaporator  is  a  Corning  Pyrex-glass  Tee  section 
(three  inch  nominal  inner  diameter)  coated  with  a  continuous 
winding  of  fiberglass  impregnated  filaments  bonded  with  a 
modified  epoxy  resin.  It  was  installed  horizontally  with  the 
side-arm  of  the  Tee  vertical  and  hydrostatically  safety  tested 
to  75  psig  by  Bertsch  [Ref.  1],  Aluminum  endplates,  210  mm 
diameter  and  12.7  mm  thick,  house  the  thermocouple  wells,  an 
oil  entry  fitting  (not  used)  and  the  Teflon  inserts,  which 
hold  the  test  tubes  in  position.  The  endplates  are  mated  to 
the  evaporator  using  cast-iron  flanges.  At  the  very  bottom  of 
the  evaporator  sits  the  refrigerant  liquid  return  line  from 
the  condenser.  A  copper  deflector  plate  was  installed  just 
above  the  return  line  to  prevent  any  vapor  bubbles  that  might 
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rise  off  the  return  line  from  impinging  on  the  test  tubes.  The 
liquid  refrigerant  level  in  the  evaporator  is  approximately  10 
mm  below  the  top  of  the  evaporator  and  approximately  10  mm 
above  the  top  of  the  upper  tube.  Figure  3.2  is  a  sketch  of 
the  evaporator  and  a  cast-iron  flange. 
2 .   Evaporator  Tubes 

Figure  3.3  is  a  schematic  of  a  boiling  tube.  Two 
different  tube  types  were  used  in  this  study:  a  smooth  tube 
and  a  TURBO-B  enhanced  surface  tube.  Table  2  provides  the 
dimensions  of  each  tube  type  while  Figure  3.4  shows  an 

TABLE  2.   DIMENSIONS  OF  TEST  TUBES 


Tube 

Dl 
mm 

D2 
mm 

Di 
mm 

Do 
mm 

L 
mm 

Lu 
mm 

kc 
W/m-K 

smooth 

12.4 

15.9 

12.7 

15.9 

203.2 

76.2 

344 

TURBO-B 

11.6 

13.8 

11.8 

15.9 

203.2 

76.2 

398 

expanded  view  of  the  enhanced  surface  of  the  TURBO-B  tube. 
Each  tube  is  heated  by  a  1  kW,  240  volt  stainless-steel 
cartridge  heater.  The  heater  measures  6.35  mm  outside 
diameter,  203.2  mm  overall  length  with  a  heated  length  of  190 
mm  and  is  centered  inside  the  test  tube.  The  tubes  are 
instrumented  with  eight  thermocouples  at  various  locations  to 
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Corning  Pyrex  Glass  Evaporator  CD  x  Q.  =  "+02x51  mm, 
L  =  17  3  mm,  I4  =  127  mm) 


Cast  Iron  Flange  (d,  =  190  mm,  d0  =  210  mm) 


Figure  3.2   Sketch  of  Pyrex  Glass  Evaporator 
and  Cast-iron  Flange 
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Figure  3 . 4   Surface  Profile  of  Enhanced  (TURBO-B)  Tube 
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monitor  the  tube  wall  temperatures.  Bertsch  [Ref.  1]  provides 
a  detailed  assembly  procedure  for  the  manufacturing  and 
instrumentation  of  a  test  tube. 

The  test  tubes  are  secured  in  the  evaporator  by  Teflon 
inserts  and  sealed  with  o-rings  (Figure  3.5).  Strong-backs, 
shown  in  Figure  3.6,  are  used  to  anchor  the  Teflon  inserts  in 
the  endplates. 

3 .   Tube  Configurations 

Figure  3.7  shows  the  six  tube  configurations  that  were 
used.  Configuration  CO  is  the  single  tube  configuration  used 
by  Bertsch  [Ref.  1]  and  Perry  [Ref.  2].  In  configurations  CI 
through  C4,  the  top  tube  is  fixed  relative  to  the  evaporator 
while  the  lower  tube  position  is  varied.  Configuration  C5 
consists  of  rotating  the  Teflon  insert  used  for  C4  through  an 
angle  of  30°.  Note  that  when  rotation  is  completed,  the  upper 
tube  position  has  moved  off  the  evaporator  centerline  and  the 
lower  tube  position  remains  unchanged. 

C.   CONDENSER  SECTION 

The  condenser  is  also  a  f iberglass-epoxy  coated  Corning 
Pyrex-glass  Tee  section  identical  in  size,  shape  and  substance 
to  the  evaporator  shown  in  Figure  3.2.  It  is  mounted 
vertically  (see  Figure  3.1)  with  aluminum  endplates  bolted  to 
cast-iron  flanges,  similar  to  the  evaporator  arrangement.  The 
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Figure  3.5   Sketch  of  an  Endplate  and  Teflon  Insert 
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Figure  3.6  Sketch  of  a  Strongback 


24 


u 


CO 
H 

Q 

a* 


a 

o 

m 

a 
m 

t 

in  § 


o 


O 


Q 

Q 

W 

04 

U 

<«_<■ 

O 

3 

tn 

J 

c 

CU 

0 
•H 
■M 

w 

CO 

Z3 

6* 

03 
M 

3 

^, 

o 

<4-l 

c 

CM 

2 

o 

• 

CJ 

~  «-* 

1 

O 

~»*. 

3 

a* 

E-« 

>—' 

CQ 

r- 

3 

• 

&* 

en 

1 

u 

^ 

3 

•H 
Cm 

o 
u 


U 


us 

Q 

a* 


25 


top  endplate  is  fitted  with  connections  for  refrigerant 
charging/vacuum  operations  and  the  ethylene-glycol/water 
coolant  inlet  tube.  The  bottom  endplate  is  fitted  with 
connections  for  the  gravity  feed  return  line  to  the  evaporator 
and  the  ethylene-glycol/water  coolant  outlet  tube.  The 
refrigerant  vapor  rises  from  the  evaporator,  passing  through 
the  L-shaped  aluminum  pipe  and  enters  the  condenser  where  it 
is  condensed  on  a  3/8  inch  outside  diameter  copper  tube  which 
is  helically  coiled  inside  the  condenser. 

D.   REFRIGERANT  STORAGE  RESERVOIR 

The  refrigerant  storage  reservoir  is  an  aluminum  cylinder, 
230  mm  in  diameter  and  254  mm  in  height.  It  is  located  within 
the  apparatus  enclosure  above  the  evaporator  and  below  the 
condenser  (see  Figure  3.1).  The  reservoir  enables  the 
operator  to  store  the  liquid  R-124  instead  of  releasing  it 
into  the  atmosphere,  prior  to  changing  out  the  test  tubes. 
The  installed  sight  glass  allows  monitoring  of  the  liquid 
level  during  replenishing  from  an  outside  source  (i.e.,  145 
lb.  storage  cylinder)  or  refilling  the  evaporator  after  a  test 
tube  changeout .  The  reservoir  is  isolated  from  or  aligned  to 
the  system  by  means  of  valves  V6,  V7,  V12  and  V13. 
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E.  OIL  SUBSECTION 

This  study  was  devoted  to  pure  R-124.  In  order  to  ensure 
that  no  oil  contaminated  the  system,  the  oil  line  from  the  oil 
cylinder  to  the  evaporator  was  removed  (see  Figure  3.1)  and 
the  access  on  the  evaporator  endplate  was  capped.  Perry 
[Ref.  2]  describes  the  oil  subsection  components  and  the  oil 
addition  procedure  for  a  refrigerant/oil  mixture  study. 

F.  COOLING  SUBSYSTEM 

1 .  Ethylene-Glycol/Water  Coolant  Sump 

At  the  base  of  the  apparatus  sits  a  0.15  cubic  meter 
plexiglass  tank  containing  approximately  30  gallons  of 
ethylene-glycol/water  mixture  producing  a  freezing  point  of 
-25 °C.  The  sump  rests  on  a  wood  platform  and  is  wrapped  with 
22  mm  thick  foam  rubber  sheet  for  insulation.  Supply  and 
discharge  connections  for  the  condenser  and  sump  pump  as  well 
as  a  thermocouple  to  monitor  the  sump  temperature  are 
installed. 

2 .  Refrigeration  Plants 

A  1/2  ton  R-502  and  1/4  ton  R-12  refrigeration  plants 
are  used  simultaneously  to  maintain  the  ethylene-glycol/water 
coolant  temperature  at  about  -15°  C.  Both  plants  are  located 
close  to  the  sump  and  contain  the  following  components:  (1) 
hermetically  sealed  compressor  assembly;  (2)  air-cooled 
condenser;    (3)  receiver;  (4)  filter-dryer;   (5)  pressure 
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regulator;     (6)   temperature   control   switch;     and   (7) 
thermostatic  expansion  valve.   The  R-502  plant  (Figure  3.8) 
uses  a  counter  flow  heat  exchanger  as  the  system  evaporator 
while  the  R-12  plant  uses  a  9 . 5  mm  diameter  coiled  copper 
tubing  submerged  in  the  sump  tank  as  the  system  evaporator 
(Figure  3.1). 

3 .   Pumps  and  Control  Valve  (VC) 

A  positive  displacement  sump  pump,  rated  at  8  gpm,  is 
connected  to  the  discharge  connection  of  the  coolant  sump  with 
a  one  inch  diameter  PVC  pipe.  It  circulates  coolant  from  the 
sump  through  the  R-502  evaporator  and  back  to  the  sump.  A 
similar  pump,  called  the  condenser  pump,  circulates  coolant 
from  the  sump  to  the  condenser,  via  the  control  valve  VC,and 
back  to  the  sump  (see  Figure  3.1) .  Valve  VC  is  used  by  the 
operator  to  control  the  refrigerant  pool  saturation 
temperature.  A  manually  operated  bypass  valve  (V9,  Figure 
3.1)  is  installed  upstream  of  valve  VC  to  prevent  overloading 
the  condenser  pump  when  valve  VC  is  closed. 

G.   APPARATUS  ENCLOSURE 

A  welded  aluminum  frame  measuring  214  cm  high,  61  cm  wide 
and  51  cm  deep  is  used  to  contain  most  of  the  system 
components.  The  upper  half  of  the  frame  contains  the 
evaporator,  condenser,  oil  subsystem,  and  refrigerant  storage 
reservoir  (see  Figure  3.1).   The  lower  half  of  the  frame 
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contains  the  coolant  sump  tank.  For  operator  safety  and  to 
reduce  heat  leakage  into  the  evaporator  and  condenser,  all 
sides  of  the  upper  half  of  the  frame  are  covered  with  13  mm 
thick  plexiglass,  hinged  on  both  sides  to  allow  access  to  the 
system  components.  Holes  in  the  front  plexiglass  cover  allow 
operation  of  valves  VI  through  V8  while  maintaining  the 
integrity  of  the  enclosed  upper  frame. 

H .   INSTRUMENTATION 

1.  Power  Measurement  System 

Figure  3.9  is  a  schematic  of  the  power  measurement 
system.  Each  boiling  tube  heater  is  powered  by  a  240  volt  AC 
source  controlled  by  a  variac.  The  variac  output  ranges  from 
0-220  volt  and  0-5  amp,  adjustable  by  the  operator  to  obtain 
the  desired  heat  flux  at  the  tube  surface.  The  current  and 
voltage  sensors  convert  the  AC  input  signal  into  DC  voltages 
which  are  used  as  input  by  the  Data  Acquisition  Unit. 

2 .  Temperature  Measurement 

Temperatures  are  measured  using  30  gauge  (0.25  mm 
diameter)  copper-constantan  (type  T)  thermocouples.  For  each 
test  tube,  eight  thermocouples  are  spaced  over  the  heater 
sleeve,  as  shown  in  Figure  3.10,  to  measure  the  wall 
temperatures.  Three  thermocouples  are  inserted  into  housings 
that  penetrate  the  evaporator  endplate  and  are  submerged  in 
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the  refrigerant  liquid  pool  at  different  locations  (see  Figure 
3.1).  Each  housing  is  constructed  of  a  low  conductivity 
stainless  steel  body  and  a  high  conductivity  copper  tip 
(Figure  3.11).  This  minimizes  the  errors  in  pool  temperature 
measurements  caused  by  heat  conduction  through  the  housing. 
A  single  thermocouple  is  located  in  the  coolant  sump  tank  to 
measure  the  sump  temperature. 

Thermocouples  were  read  by  the  a  Hewlett-Packard  (HP) 
3497A  Data  Acquisition  Unit  which  is  controlled  by  an  HP  326 
Computer  System. 

I.   DATA  ACQUISITION  AND  REDUCTION  PROGRAM 

Data  acquisition  and  reduction  program  DRP72,  used  by  Lake 
[Ref.  3],  was  modified  to  incorporate  the  thermophysical 
properties  of  saturated  R-124.  Bertsch  [Ref.  1]  discusses  the 
generation  of  the  equations  used  in  the  program  calculations. 
The  modified  DRP72  program  was  renamed  DRPJY.  Appendix  B 
contains  the  program  DRPJY.  The  program  DRPJY  directs  the  HP 
362  Computer  System  and  the  HP  3497A  Data  Acquisition  Unit  to 
read  and  process  all  thermocouple,  current  and  R.M.S.  voltage 
readings.  Using  this  information,  DRPJY  performs  the 
following  calculations  (from  Bertsch  [Ref.  1]): 

1.  Converts  voltage  readings  into  tube  temperatures  and 
power. 

2.  Computes  heat  transfer  rate  from  cartridge  heaters. 
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3.  Computes  average  test  tube  wall  temperatures  and  wall 
superheat  (Twall  -  Tsat) . 

4.  Computes  physical  properties  of  R-124  using  property 
correlations  from  RefProp  [Ref.  28]  at  the  film 
temperature  (Twall  +  Tsat)/2. 

5.  Computes  heat  transfer  coefficient  of  R-124  from  the 
unheated  ends  and  heated  length  of  the  test  tubes. 

6.  Computes  heat  losses  from  the  unheated  ends  of  the  test 
tubes. 

7.  Computes  heat  flux  from  heated  length  of  test  tubes. 


The  processed  output  data  was  printed  by  an  HP  Inkjet 
printer  and  stored  on  the  computer  hard  drive.  As  a 
precaution,  each  output  data  file  was  copied  to  a  3.5  inch 
floppy  disk  in  the  event  the  computer  hard  drive  becomes 
unreadable.  Table  3  lists  the  channel  assignments  that  are 
sensor  inputs  to  the  Data  Acquisition  Unit. 
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TABLE  3.   HP-3497A  CHANNEL  ASSIGNMENTS 


Channel 

Channel  Assignments 

00-07 

Upper  test  tube  wall  temperature 

(T1-T8) 

08-15 

Lower  test  tube  wall  temperature 

(T9-T16) 

16 

Liquid  pool  thermocouple  (Tldl) 

17 

Liquid  pool  thermocouple  (Tld2) 

18 

Liquid  pool  thermocouple  (Tv) 

19 

Sump  coolant  temperature 

20 

R.M.S.  voltage  of  tube  heaters 

21 

Upper  tube  current  sensor 

22 

Lower  tube  current  sensor 
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IV.    EXPERIMENTAL  PROCEDURE 

A.   EVAPORATOR  DISASSEMBLY  PROCEDURE 

The  following  procedure  is  used  to  transfer  the  liquid 
refrigerant  from  the  evaporator  to  the  reservoir  to  facilitate 
test  tube  changeouts: 


1.  Operate  both  R-502  and  R-12  refrigeration  systems  to 
cool  the  sump  tank  temperature  to  approximately  -15  °C . 

2.  Start  the  condenser  pump  and  using  control  valve  VC, 
adjust  the  coolant  supply  pressure  to  1  psig. 

3.  Close  valve  V5  to  isolate  the  evaporator  from  the 
condensing  refrigerant  vapor  (see  Figure  3.1). 

4.  Align  the  refrigerant  storage  reservoir  to  the  system  by 
opening  valves  V6,  V7,  V12  and  V13. 

5.  Position  and  energize  the  high  volume  air  dryer  (set  for 
high  heat)  so  the  heated  air  stream  is  blowing  on  the 
center  of  the  evaporator.  Adjust  valve  VC  to  ensure  the 
refrigerant  vapor  is  condensing. 

6.  After  1-2  hours,  the  evaporator  should  be  dry  with  the 
liquid  refrigerant  in  the  reservoir  and  a  small  amount 
in  the  return  line. 

7.  Close  valves  V6,  V7,  V12  and  V13,  isolating  the 
reservoir  from  the  system. 

8.  Secure  and  remove  the  air  dryer. 

9.  Stop  the  condenser  pump. 

10.  Open  valves  V8  and  V10  to  reduce  to  system  pressure  to 
atmospheric.   Close  valve  V8  and  V10  when  the  system 
pressure  reads  0  psig  on  the  system  pressure  gage. 

11.  Disconnect  the  test  tube  heater  and  thermocouple  leads. 
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12.  Remove  strongbacks. 

13.  Remove  Teflon  inserts  and  test  tubes. 

B.   EVAPORATOR  ASSEMBLY  PROCEDURE 

1 .  Evaporator  Assembly 

Before  installing  smooth  test  tubes,  rinse  with 
acetone  then  wipe  clean  with  lint-free  wipes  or  if  installing 
TURBO-B  test  tubes,  brush  the  enhanced  surfaces  with  a  soft 
bristle  brush  (i.e.,  toothbrush)  using  ethyl  alcohol  then 
repeat  the  brushing  using  acetone  and  let  the  tubes  air  dry. 
Check  the  Teflon  insert  o-rings  for  damage  and  replace  as 
necessary.  Before  installing  the  test  tubes,  wipe  the 
evaporator  interior  clean  with  lint-free  wipes.  Install  the 
test  tubes  in  the  evaporator  using  the  Teflon  inserts  to 
secure  the  tubes  in  place.  Snug  the  strongbacks  up  to  the 
Teflon  inserts  and  connect  the  thermocouple  leads.  Prior  to 
connecting  the  heater  leads,  check  the  resistance  between  the 
heater  leads  and  the  outer  tube  wall  using  a  multimeter.  One 
megohm  resistance  or  greater  is  satisfactory. 

2.  Leak  Test 

After  assembling  the  evaporator  and  before  charging 
the  system  with  refrigerant,  perform  a  vacuum  test  to  verify 
the  integrity  of  the  system.  Evacuate  the  system  by  starting 
the  vacuum  pump  and  opening  valves  VI 1  and  V8  (see  Figure 
3.1).   After  evacuating  the   system  to   approximately  29 
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in.  Hg,  close  valves  Vll  and  V8,  then  stop  the  vacuum  pump. 
If  no  drop  in  vacuum  is  detected  (i.e.,  no  leak)  within  30 
minutes,  charge  the  system  with  R-124. 

If  a  leak  is  detected,  open  valves  V8  and  Vll  and 
pressurize  the  system  with  air  to  20  psig.  Apply  a  soap/water 
solution  to  the  system  sealing  surfaces  to  identify  the  leak 
location  (bubbles  will  be  seen  from  the  leak  source) .  After 
identifying  and  repairing  the  leak,  retest  the  system  for 
leaks . 

3.   Charging  System  with  Refrigerant 

Once  the  system  integrity  is  assured,  fill  the 
evaporator  with  R-124,  using  one  of  the  following  procedures, 
to  a  level  just  Jbelov  the  top  of  the  evaporator  (NOTE:  some 
refrigerant  will  be  lost  during  the  degas  procedure, 
therefore,  it  is  important  to  fill  the  evaporator  above  the 
level  needed  to  conduct  data  runs) : 

a.  From  Storage  Reservoir 

Open  valves  V12  and  V7  to  equalize  the  pressure 
between  the  storage  reservoir  and  the  system.  Open  valves 
V13,  V6  and  V5  to  gravity  fill  the  evaporator  with  liquid 
refrigerant.  When  the  evaporator  is  filled  to  the  proper 
level,  close  valves  V6,  V13,  V7  and  V12  to  isolate  the  storage 
reservoir  from  the  system. 
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Jb.  From  Refrigerant  Cylinder 

Connect  the  cylinder  to  valve  V10  using  the 
charging  hose  (Figure  3.1).  Open  valves  V5,  V8  and  V10  and 
then  open  the  cylinder  liquid  valve  to  release  the  refrigerant 
from  the  cylinder  into  the  system.  Close  valves  V8,  V10  and 
the  cylinder  liquid  valve  when  the  evaporator  is  filled  to  the 
proper  level.  Disconnect  the  charging  hose  from  valve  V10  and 
stow  the  cylinder. 

4 .   Degassing  Procedure 

This  degas  procedure  is  similar  to  Perry's  [Ref.  3] 
degassing  procedure  B.  Degas  the  system  by  boiling  the 
refrigerant,  using  a  upper  test  tube  heat  flux  of  85  KW/m2  and 
a  lower  test  tube  heat  flux  of  20  KW/m2,  for  30  minutes. 
Ensure  the  condenser  pump  is  running  and  control  valve  VC  is 
open  to  prevent  overpressurizing  the  system.  During  the 
boiling  process,  air  and  non-condensible  gases  rise  in  the 
system  and  collect  in  the  condenser.  After  the  30  minute 
boiling  time  has  elapsed,  secure  the  heat  flux  to  the  test 
tubes  and  stop  the  condenser  pump.  Let  the  system  set  for  30 
minutes  allowing  the  condensed  refrigerant  to  drain  to  the 
evaporator.  After  the  30  minute  set  time  has  elapsed,  start 
the  vacuum  pump  and  open  valve  VI 1,  then  open  valve  V8  for  20 
seconds  to  remove  the  accumulated  air  and  non-condensible 
gases.  After  closing  valve  V8,  stop  the  vacuum  pump  and  close 
valve  Vll. 
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5.   Data  Acquisition  Channel  Check 

After  all  evaporator  thermocouples  (test  tubes  and 
pool)  are  connected,  load  and  run  program  SETUPJY  to  verify 
the  output  of  each  thermocouple  is  correct.  Any  output  that 
is  inconsistent  or  suspect  is  either  fixed  or  noted  and  later 
excluded  from  the  calculations  in  program  DRPJY. 

C.   STANDARD  OPERATING  PROCEDURE 

The  following  procedure  is  used  to  conduct  each  data 
collection  run: 


1.  Operate  the  R-502  and  R-12  refrigeration  systems  to  cool 
the  sump  to  approximately  -15 °C. 

2.  Energize  the  HP-362  Computer  System,  Data  Acquisition 
Unit  and  the  variac  control  panel. 

3.  Load  and  run  program  SETUPJY.  Verify  the  thermocouple 
outputs  are  consistent. 

4.  Verify  control  valve  VC  is  closed,  then  start  the 
condenser  pump. 

5.  Slowly  open  valve  VC  until  the  coolant  pressure  gage 
reads  1/8  psig.  This  pressure  setting  will  ensure  that 
the  refrigerant  pool  will  be  slowly  cooled. 

6.  Monitor  the  thermocouple  outputs  with  SETUPJY  until 
confident  the  thermocouples  are  operating  properly. 
Note  which  thermocouples  are  inconsistent  and  must  be 
excluded  from  program  DRPJY. 

7.  When  the  refrigerant  pool  temperature  is  approaching  5 
°C,  load  and  run  program  DRPJY.   Select  'TAKE  DATA' 
option. 

8.  At  the  prompt,  enter  the  filename  that  the  data  is  to  be 
stored  under  (page  38  explains  the  file  naming  scheme 
used) . 
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9.  At  the  prompt,  enter  the  number  of  defective  upper  tube 
thermocouples  that  were  identified  by  SETUP JY. 

10.  At  the  prompt,  enter  the  defective  upper  tube 
thermocouple  locations  that  were  identified  by  SETUPJY. 

11.  At  the  prompt,  enter  the  number  of  defective  lower  tube 
thermocouples  that  were  identified  by  SETUPJY. 

12.  At  the  prompt,  enter  the  defective  lower  tube 
thermocouple  locations  that  were  identified  by  SETUPJY. 

13.  Select  the  type  of  test  tubes  installed. 

14.  From  the  main  menu,  select  'SET  Tsat'  .  Monitor  the  pool 
thermocouples  and,  if  required,  adjusted  valve  VC  to 
maintain  the  desired  pool  saturation  temperature  of  2.22 
±  0.1  °C. 

15.  From  the  main  menu,  select  'SET  UPPER  TUBE  HEAT  FLUX', 
then  adjust  the  main  variac  to  the  desired  heat  flux. 
Initial  heat  flux  setting  is  500-600  KW/m2 . 

16.  Press  the  keyboard  button  fl  to  return  to  the  main  menu 
and  select  'SET  LOWER  TUBE  HEAT  FLUX'.  Adjust  the  aux 
variac  to  the  desired  heat  flux. 

17.  From  the  main  menu,  select  'SET  Tsat'.  Adjust  valve  VC 
to  maintain  the  desired  pool  saturation  temperature 
while  monitoring  the  pool  thermocouples. 

18.  Once  the  pool  saturation  temperature  is  stable,  return 
to  the  main  menu  and  select  'TAKE  DATA'. 

19.  For  each  heat  flux  setting  steps  15-18  are  repeated. 


D.   DATA  FILES 

Data  were  filed  using  the  following  convention: 

Example  (D0904ASC1)  : 

1.  First  letter  D  indicates  the  file  is  a  data  file. 

2.  The  following  four  characters  indicate  the  date;  in  this 
example  September  04,  1993. 
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3.  The  next  letter  indicates  the  order  in  which  the  data  run 
was  completed  on  that  day/  no  letter  -  1st  run,  A  -  2nd 
run,  B  -  3rd  run,  etc. 

4.  The  next  one  or  two  letters  indicate  the  tubes  used; 
S  -  smooth,  TB  -  TURBO-B. 

5.  The  last  letter  and  number  combination  represent  the  tube 
configuration  used. 


Table  4  is  a  summarized  list  of  the  data  runs  completed 
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TABLE  4 


LIST  OF  DATA  RUNS 


Data  File 

Tube  Type 

Purpose 

D0728ASC0 

Smooth 

Repeatability;  single  tube 

D0804ASC0 

Smooth 

Repeatability;  single  tube 

D0816SC0 

Smooth 

Repeatability;  single  tube 

D0520SC1 

Smooth 

Repeatability;  P/D  =  2.0;  no 
power  on  lower  tube 

D0714SC1 

Smooth 

Repeatability;  P/D  =  2.0;  no 
power  on  lower  tube 

D0831SC1 

Smooth 

Data;  repeatability;  P/D  =  2.0; 
no  power  on  lower  tube 

D0901SC1 

Smooth 

Data;  P/D  =2.0;  1  kW/m2  on 
lower  tube 

D0902SC1 

Smooth 

Data;  P/D  =  2.0;  5  kW/m2  on 
lower  tube 

D0902BSC1 

Smooth 

Data;  P/D  =  2.0;  20  kW/m2  on 
lower  tube 

D0903SC1 

Smooth 

Data;  P/D  =  2.0;  q"  on  both 
tubes  inc/dec  together 

D0904SC1 

Smooth 

Data;  P/D  =  2.0;  5  kW/m2  on 
upper  tube 

D0904ASC1 

Smooth 

Data;  P/D  =  2.0;  10  kW/m2  on 
upper  tube 

D0905SC1 

Smooth 

Data;  P/D  =  2.0;  20  kW/m2  on 
upper  tube 

D0906SC2 

Smooth 

Data;  P/D  =  1.8;  no  power  on 
lower  tube 

D0906ASC2 

Smooth 

Data;  P/D  =  1.8;  1  kW/m2  on 
lower  tube 

D0907SC2 

Smooth 

Data;  P/D  =  1.8;  5  kW/m2  on 
lower  tube 

D0907ASC2 

Smooth 

Data;  P/D  =  1.8;  10  kW/m2  on 
lower  tube 
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TABLE  4.   LIST  OF  DATA  RUNS  (continued) 


Data  File 

Tube  Type 

Purpose 

D0908SC2 

Smooth 

Data; 
lower 

P/D  = 
tube 

1.8; 

20  kW/m2  on 

D0908ASC2 

Smooth 

Data; 
tubes 

P/D  =  1.8; 
inc/dec  to 

q"  on  both 

gether 

D0908BSC2 

Smooth 

Data; 
upper 

P/D  = 
tube 

1.8; 

5  kW/m2  on 

D0909SC2 

Smooth 

Data; 
upper 

P/D  = 
tube 

1.8; 

10  kW/m2  on 

D0909ASC2 

Smooth 

Data; 

upper 

P/D  = 
tube 

1.8; 

2  0  kW/m2  on 

D0910SC3 

Smooth 

Data; 
lower 

P/D  = 
tube 

1.5; 

no  power  on 

D0910ASC3 

Smooth 

Data; 
lower 

P/D  = 
tube 

1.5; 

1  kW/m2  on 

D0917SC3 

Smooth 

Data; 
lower 

P/D  = 
tube 

1.5; 

5  kW/m2  on 

D0919SC3 

Smooth 

Data; 
lower 

P/D  = 
tube 

1.5; 

10  kW/m2  on 

D0919ASC3 

Smooth 

Data; 
lower 

P/D  = 
tube 

1.8; 

20  kW/m2  on 

D0920SC3 

Smooth 

Data; 

tubes 

P/D  =  1.8; 
inc/dec  to 

q"  on  both 
gether 

D0920ASC3 

Smooth 

Data; 
upper 

P/D  = 

tube 

1.5; 

5  kW/m2  on 

D0920BSC3 

Smooth 

Data; 

upper 

P/D  = 
tube 

1.5; 

10  kW/m2  on 

D0920CSC3 

Smooth 

Data; 
upper 

P/D  = 

tube 

1.5; 

20  kW/m2  on 

D0921SC4 

Smooth 

Data; 
lower 

P/D  = 

tube 

1.2; 

no  power  on 

DQ921ASC4 

Smooth 

Data; 

lower 

P/D  = 
tube 

1.2; 

1  kW/m2  on 

D0921BSC4 

Smooth 

Data; 
lower 

P/D  = 
tube 

1.2; 

5  kW/m2  on 
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TABLE  4.   LIST  OF  DATA  RUNS  (continued) 


Data  File 

Tube  Type 

Purpose 

D0921CSC4 

Smooth 

Data;  P/D  =  1.2;  10  kW/m2  on 
lower  tube 

D0922SC4 

Smooth 

Data;  P/D  =  1.2;  20  kW/m2  on 
lower  tube 

D0922ASC4 

Smooth 

Data;  P/D  =  1.2;  q"  on  both 
tubes  inc/dec  together 

D0922BSC4 

Smooth 

Data;  P/D  =  1.2;  5  kW/m2  on 
upper  tube 

D0922CSC4 

Smooth 

Data;  P/D  =  1.2;  10  kW/m2  on 
upper  tube 

D0922DSC4 

Smooth 

Data;  P/D  =  1.2;  20  kW/m2  on 
upper  tube 

D0930SC5 

Smooth 

Data;  P/D  =1.2  rotated  30°;  no 
power  on  lower  tube 

D0930ASC5 

Smooth 

Data;  P/D  =1.2  rotated  30°;  1 
kW/m2  on  lower  tube 

D0930BSC5 

Smooth 

Data;  P/D  =1.2  rotated  30°;  5 
kW/m2  on  lower  tube 

D0930CSC5 

Smooth 

Data;  P/D  =1.2  rotated  30°;  10 
kW/m2  on  lower  tube 

D0929CSC5 

Smooth 

Data;  P/D  =1.2  rotated  30°;  20 
kW/m2  on  lower  tube 

D0929DSC5 

Smooth 

Data;  P/D  =1.2  rotated  30°;  q" 
on  both  tubes  inc/dec  together 

D0929BSC5 

Smooth 

Data;  P/D  =1.2  rotated  30°;  5 
kW/m2  on  upper  tube 

D0929ASC5 

Smooth 

Data;  P/D  =1.2  rotated  30°;  10 
kW/m2  on  upper  tube 

D0929SC5 

Smooth 

Data;  P/D  =1.2  rotated  30°;  20 
kW/m2  on  upper  tube 

D1012TBC0 

TURBO-B 

Repeatability;  single  tube 

D1012ATBC0 

TURBO-B 

Repeatability;  single  tube 

D1012BTBC0 

TURBO-B 

Repeatability;  single  tube 
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TABLE  4.   LIST  OF  DATA  RUNS  (continued) 


Data  File 

Tube  Type 

Purpose 

D1018TBC1 

TURBO-B 

Repeatability;  P/D  =  2.0;  no 
power  on  lower  tube 

D1018ATBC1 

TURBO-B 

Repeatability;  P/D  =  2.0;  no 
power  on  lower  tube 

D1018BTBC1 

TURBO-B 

Data;  repeal 
no  power  on 

lability;  P/D  =  2.0; 
lower  tube 

D1018CTBC1 

TURBO-B 

Data; 
lower 

P/D  = 
tube 

2.0;  1  kW/m2  on 

D1018DTBC1 

TURBO-B 

Data; 

lower 

P/D  = 
tube 

2.0;  5  kW/m2  on 

D1019TBC1 

TURBO-B 

Data; 
lower 

P/D  = 
tube 

2.0;  10  kW/m2  on 

D1019ATBC1 

TURBO-B 

Data; 
lower 

P/D  = 
tube 

2.0;  20  kW/m2  on 

D1019BTBC1 

TURBO-B 

Data; 
tubes 

P/D  =  2.0;  q"  on  both 
decreased  together 

D1019CTBC1 

TURBO-B 

Data; 
upper 

P/D  = 
tube 

2.0;  5  kW/m2  on 

D1019DTBC1 

TURBO-B 

Data; 
upper 

P/D  = 
tube 

2.0;  10  kW/m2  on 

D1019ETBC1 

TURBO-B 

Data; 
upper 

P/D  = 
tube 

2.0;  20  kW/m2  on 

D1020TBC2 

TURBO-B 

Data; 
lower 

P/D  = 
tube 

1.8;  no  power  on 

D1020ATBC2 

TURBO-B 

Data; 
lower 

P/D  = 

tube 

1.8;  1  kW/m2  on 

D1020BTBC2 

TURBO-B 

Data; 
lower 

P/D  = 
tube 

1.8;  5  kW/m2  on 

D1020CTBC2 

TURBO-B 

Data; 
lower 

P/D  = 
tube 

1.8;  10  kW/m2  on 

D1020DTBC2 

TURBO-B 

Data; 
lower 

P/D  = 
tube 

1.8;  20  kW/m2  on 

47 


TABLE  4.   LIST  OF  DATA  RUNS  (continued) 


Data  File 

Tube  Type 

Purpose 

D1020ETBC2 

TURBO-B 

Data; 

tubes 

P/D  =  1.8; 
decreased 

q"  on  both 
together 

D1021TBC2 

TURBO-B 

Data; 
upper 

P/D  = 
tube 

1.8; 

5  kW/m2  on 

D1021ATBC2 

TURBO-B 

Data; 
upper 

P/D  = 
tube 

1.8; 

10  kW/m2  on 

D1021BTBC2 

TURBO-B 

Data; 
upper 

P/D  = 
tube 

1.8; 

20  kW/m2  on 

D1021CTBC3 

TURBO-B 

Data; 
lower 

P/D  = 

tube 

1.5; 

no  power  on 

D1021DTBC3 

TURBO-B 

Data; 

lower 

P/D  = 
tube 

1.5; 

1  kW/m2  on 

D1021ETBC3 

TURBO-B 

Data; 
lower 

P/D  = 
tube 

1.5; 

5  kW/m2  on 

D1021FTBC3 

TURBO-B 

Data; 
lower 

P/D  = 
tube 

1.5; 

10  kW/m2  on 

D1022TBC3 

TURBO-B 

Data; 
lower 

P/D  = 
tube 

1.5; 

2  0  kW/m2  on 

D1022ATBC3 

TURBO-B 

Data; 
tubes 

P/D  =  1.5; 
decreased 

q"  on  both 
together 

D1022BTBC3 

TURBO-B 

Data; 
upper 

P/D  = 
tube 

1.5; 

5  kW/m2  on 

D1022CTBC3 

TURBO-B 

Data; 
upper 

P/D  = 

tube 

1.5; 

10  kW/m2  on 

D1022DTBC3 

TURBO-B 

Data; 
upper 

P/D  = 
tube 

1.5; 

20  kW/m2  on 

D1022ETBC4 

TURBO-B 

Data; 
lower 

P/D  = 
tube 

1.2; 

no  power  on 

D1022FTBC4 

TURBO-B 

Data; 

lower 

P/D  = 

tube 

1.2; 

1  kW/m2  on 

D1025TBC4 

TURBO-B 

Data; 
lower 

P/D  = 
tube 

1.2; 

5  kW/m2  on 

D1025ATBC4 

TURBO-B 

Data; 
lower 

P/D  = 
tube 

1.2; 

10  kW/m2  on 
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TABLE  4.   LIST  OF  DATA  RUNS  (continued) 


Data  File 

Tube  Type 

Purpose 

D1025BTBC4 

TURBO-B 

Data; 
lower 

P/D  =  1.2;  20  kW/m2  on 

tube 

D1025CTBC4 

TURBO-B 

Data; 
tubes 

P/D  =  1.2;  q"  on  both 
decreased  together 

D1025DTBC4 

TURBO-B 

Data; 
upper 

P/D  =  1.2;  5  kW/m2 
tube 

on 

D1025ETBC4 

TURBO-B 

Data; 
upper 

P/D  =  1.2;  10  kW/m2  on 
tube 

D1025FTBC4 

TURBO-B 

Data; 
upper 

P/D  =  1.2;  20  kW/m2  on 
tube 

D1026TBC5 

TURBO-B 

Data; 
power 
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V.    DISCUSSION  OF  RESULTS 

A.   REPRODUCIBILITY/REPEATABILITY 

To  verify  that  the  equipment  and  operating  procedure  being 
used  were  consistent  with  those  used  by  Bertsch  [Ref.  1]  and 
Perry  [Ref.  2],  six  reproducibility/repeatability  runs  were 
performed.  Program  DRPGP  (from  Perry's  single  tube  work 
[Ref.  2])  was  used  for  the  runs  conducted  with  a  single  test 
tube  and  DRPJY  was  used  for  all  runs  conducted  with  two  test 
tubes.  Figures  5.1  and  5.2  are  plots  of  increasing  and 
decreasing  heat  flux  using  a  smooth  tube  in  configuration  CO 
(which  is  identical  to  the  configuration  used  by  Perry 
[Ref.  2]).  In  Figure  5.1,  the  three  runs  completed  are  in 
very  good  agreement  with  each  other.  At  lower  heat  fluxes 
(  <  4  kW/m2)  some  scatter  is  present  which  disappears  at  mid- 
range  heat  fluxes  (4-8  kW/m2)  .  Perry  [Ref.  2]  found  the  same 
phenomenon  and  attributed  the  scatter  to  convection  regime 
thermal  gradients  where  mixing  is  inefficient  within  the 
liquid  pool,  after  discovering  large  local  variations  in  the 
tube  wall  and  pool  thermocouples  of  up  to  5%.  In  this  study, 
local  variations  in  the  tube  wall  and  pool  thermocouples  were 
found  to  be  up  to  7%  at  lower  heat  fluxes  and  decreasing  as 
heat  flux  increased,  which  supports  Perry's  theory.   In 
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Figure  5.1   Increasing  Heat  Flux  Repeatability  Runs  using 
Smooth  Tube  Configuration  CO 
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natural  convection,  the  three  data  runs  started  to  the  left  of 
the  Churchill  and  Chu  correlation  [Ref.  9]  (shown  as  a  solid 
line)  at  low  heat  fluxes  and  converged  toward  the  Churchill- 
Chu  line  as  the  heat  flux  was  increased.  This  slight  shift  to 
the  left,  although  consistent  during  the  three  runs,  was 
within  the  uncertainty  in  the  region.  Good  repeatability  was 
also  found  at  the  onset  of  nucleate  boiling  (indicated  by 
arrows)  and  in  the  nucleate  boiling  region.  Comparing  the 
repeatability  runs  to  the  representative  data  set  of  Perry 
[Ref.  2],  at  low  heat  fluxes  the  data  vary  slightly  but  are 
within  the  uncertainty  band.  The  data  merge  up  to  the  point 
of  incipience,  but  why  Perry's  run  nucleated  earlier  than  the 
present  data  is  unclear  since  the  runs  were,  in  essence, 
identical.  In  the  boiling  region,  Perry's  data  is  slightly  to 
the  left  of  the  present  data  but  again  lie  within  the 
uncertainty  band  for  that  region. 

In  Figure  5.2,  good  repeatability  is  again  shown  for  all 
heat  fluxes.  Comparing  the  three  present  repeatability  runs 
to  the  representative  data  set  of  Perry  [Ref.  2],  again  shows 
good  reproducibility  throughout.  For  heat  fluxes  <  1  kW/m2, 
the  slopes  of  the  data  curves  change,  which  indicates  the 
transition  from  nucleate  boiling  back  to  natural  convection. 
Based  on  Figures  5.1  and  5.2,  data  reproducibility  is 
considered  satisfactory  with  the  single  smooth  tube  work 
carried  out  by  Perry  [Ref.  2]  on  the  same  apparatus. 
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Figure  5.3  is  a  comparison  of  decreasing  heat  flux 
repeatability  runs  using  an  enhanced  tube  (TURBO-B)  in 
configuration  CO  with  the  single  tube  data  of  Bertsch 
[Ref.  1].  For  the  present  data,  good  repeatability  is  again 
obtained  throughout.  In  comparing  the  current  data  with  the 
data  collected  by  Bertsch  [Ref.  1],  different  boiling  curves 
are  easily  seen.  This  was  a  result  of  two  different  TURBO-B 
test  tubes  being  used.  After  examining  both  test  tubes  under 
a  microscope  (using  50x  and  lOOx  magnification) ,  the  TURBO-B 
test  tube  used  in  this  study  had  wider  notches  between  the 
'mushroom-like  pedestals'  (see  Figure  3.4),  indicating  that  a 
different  fabrication  process  was  used  to  make  the  tube. 
Above  20  kW/m2  on  the  plot,  the  current  data  curves  are  to  the 
left  of  Bertsch's  data  (about  2  °C  less  wall  superheat), 
indicating  better  heat  transfer.  It  is  believed  this  is  due 
to  the  wider  notches  in  the  tube  surface  allowing  more  bubbles 
to  depart  the  tube  surface  at  higher  heat  fluxes,  removing 
more  heat  than  the  tube  used  by  Bertsch.  At  about  10  kW/m2, 
the  trend  is  reversed,  the  current  data  crosses  and  continues 
to  fall  below  Bertsch's  data  (about  0.3  °C  more  wall 
superheat) .  At  low  heat  fluxes,  the  wider  notches  allow  the 
refrigerant  to  flood  back  into  the  boiling  channels  (see 
Figure  3.4)  and  reduce  the  number  of  active  sites,  reducing 
the  heat  transfer  from  the  tube. 

Figure  5.4  is  a  comparison  of  increasing  heat  flux 
repeatability  runs  using  smooth  tube  configuration  CI  with  the 
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TURBO-B  Tube  Configuration  CO 
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Smooth  Tube  Configuration  CI  with  the  Lower  Tube  Unheated 
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lower  tube  unheated  and  a  single  smooth  tube  data  set 
(D0728ASC0)  from  Figure  5.1.  The  plot  shows  some  scatter  of 
data  at  low  heat  fluxes  that  is  attributed  to  the  liquid  pool 
temperature  gradients,  as  discussed  previously.  Again,  good 
repeatability  is  found  in  the  mid-range  heat  flux  and  boiling 
regions.  In  comparison  to  the  single  tube  data  (configuration 
CO)  and  the  Churchill  and  Chu  correlation  [Ref.  9],  the  upper 
tube  heat  transfer  performance  of  the  CI  configuration  is 
clearly  enhanced  in  natural  convection.  This  enhancement  is 
probably  due  to  the  liquid  pool  surface  being  much  closer  to 
the  upper  test  tube  (approximately  10  mm)  when  it  is  position 
CI,  increasing  the  circulation  (i.e.,  local  fluid  velocity) 
around  the  upper  tube  and  ultimately  increasing  heat  transfer 
performance  of  the  tube.  Lake  [Ref.  3],  in  his  comparison  of 
single  tube  R-114  data  to  his  CI  configuration  data,  found  the 
same  enhancement.  He  concluded  that  the  lower  tube  must  be 
modifying  the  flow  over  the  upper  tube.  In  actuality,  it  is 
probably  a  combination  of  both  the  decreased  distance  to  the 
pool  surface  and  the  lower  tube  modifying  the  flow  over  the 
upper  tube  that  causes  the  enhancement.  From  this  plot,  it 
was  concluded  that  the  program  DRPJY  was  performing  properly. 
Figure  5.5  is  the  corresponding  decreasing  heat  flux  data 
runs  for  Figure  5.4.  Good  repeatability  is  again  found  for 
the  CI  configuration  runs.  Comparing  the  single  tube  data 
(D0728ASC0  from  Figure  5.2)  to  the  Cl  configuration  data, 
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indicates  that  in  nucleate  boiling  with  decreasing  heat 
fluxes,  both  tube  configurations  perform  equally  well. 

Figure  5.6  is  the  decreasing  heat  flux  repeatability  runs 
using  TURBO-B  tubes  in  configuration  CI  with  the  lower  tube 
unheated.  Good  repeatability  is  evident.  Comparing  the 
single  TURBO-B  data  set  (D1012TBC0)  from  Figure  5 . 3  to  the  CI 
data,  at  high  heat  fluxes  both  tube  configurations  perform 
equally  well,  as  was  for  the  smooth  tube  case  (Figure  5.5). 
At  low  heat  fluxes  (<  8  kW/m2)  ,  the  CI  configuration  data  is 
to  the  left  of  the  CO  configuration  data,  indicating  an 
enhancement  due  to  the  addition  of  the  lower  tube  and  moving 
the  test  tube  closer  to  the  surface  of  the  liquid  pool;  this 
was  not  seen  for  the  smooth  tube  case  (Figure  5.5) .  As  will 
be  discussed  later,  the  addition  of  a  lower  tube  heat  flux, 
shifts  the  CI  data  curve  even  more  towards  the  left. 

Based  on  Figures  5.1  through  5.6,  the  repeatability  of  the 
current  data  and  reproducibility  of  previous  single  tube  work 
performed  by  Perry  [Ref.  2],  is  considered  satisfactory. 

B.   INFLUENCE  OF  THE  LOWER  TUBE 

To  investigate  the  effects  of  a  lower  heated  tube  on  an 
upper  tube,  five  different  tube  configurations  and  fixed  lower 
tube  heat  flux  settings  were  used.  The  five  tube 
configurations  selected  were  P/D  ratios  of  2.0,  1.8,  1.5,  1.2 
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Figure  5 . 6   Decreasing  Heat  Flux  Repeatability  Runs  using 
TURBO-B  Tube  Configuration  CI  with  the  Lower  Tube  Unheated 
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and  1.2  (rotated  30°).  P/D  of  2.0  was  selected  because  it 
matches  not  only  the  vertical  pitch  used  in  naval  shipboard 
flooded  evaporators,  but  also  the  pitch  used  in  the  flooded 
evaporator  bundle  apparatus  at  the  Naval  Postgraduate  School 
Heat  Transfer  Laboratory.  The  fixed  lower  tube  heat  flux 
settings  (LTHFS)  used  were  0,  1,  5,  10  and  20  kW/m2.  In 
addition  to  these  runs,  data  runs  for  each  tube  configuration 
were  also  done  using  a  lower  tube  heat  flux  setting  that  was 
equal  to  the  upper  tube  heat  flux  setting  (UTHFS) .  Finally, 
at  a  fixed  UTHFS  of  5,  10  and  20  kW/m2,  the  LTHFS  was  varied 
to  investigate  the  influence  of  the  strength  of  the  convection 
effects  from  below.  For  each  data  run,  the  same  tube  surface 
was  used  for  the  upper  and  lower  tube;  mixing  of  the  tube 
types  was  not  allowed.  For  the  smooth  tube  runs,  the  upper 
tube  was  started  in  the  natural  convection  region  at  low  heat 
fluxes,  was  run  up  and  down  through  various  heat  flux 
settings,  and  was  returned  to  the  natural  convection  region. 
For  the  TURBO-B  tubes,  only  decreasing  upper  tube  heat  flux 
runs  were  recorded  to  prevent  overlapping/repeating  data  (the 
TURBO-B  tubes  did  not  enter  the  natural  convection  region  at 
low  heat  fluxes) . 

To  evaluate  the  effect  of  a  lower  heated  tube  on  an  upper 
tube,  the  data  is  presented  by  the  effect  of  changing  the  P/D 
ratio  for  a  common  LTHFS,  the  effect  of  varying  the  LTHFS  for 
a  common  P/D  ratio  and  the  effect  of  varying  the  lower  tube 
heat  flux  for  a  fixed  upper  tube  heat  flux  setting. 
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1 .   The  Effect  of  Tube  Spacing 

a.  Smooth   Tubes   \rith   Lover   Tube   Unheated 

Figure  5.7  shows  the  increasing  heat  flux  data 
taken  for  the  five  tube  configurations  (C1-C5)  with  the  lower 
tube  unheated  and  the  corresponding  Churchill  and  Chu 
correlation  [Ref .  9] .  As  seen  in  the  repeatability  runs,  at 
low  heat  fluxes  the  data  is  shifted  to  the  left  of  the 
Churchill-Chu  line  and  as  the  heat  flux  is  increased,  the  data 
tends  to  approach  the  Churchill-Chu  line  up  to  the  incipient 
boiling  point  (indicated  by  arrows) .  This  implies  that  the 
circulation  occurring  within  the  evaporator  enhances  the 
natural  convection  that  occurs  from  a  single  tube  immersed  in 
an  infinite  expanse  of  fluid.  Following  nucleation  of  the 
upper  tube,  the  boiling  curve  for  P/D  of  1.8  is  to  the  left  of 
the  others  and  appears  to  be  yielding  the  best  heat  transfer 
performance.  It  is  unclear  why  this  is  so,  because  it  is  so 
repeatable  and  the  same  smooth  test  tubes  were  used.  Lake 
[Ref.  3]  also  found  that  P/D  of  1.8  gave  the  best  heat 
transfer  performance. 

Figure  5.8  shows  the  corresponding  decreasing  heat 
flux  data  for  Figure  5.7.  With  decreasing  heat  flux,  a  P/D  of 
1.8  again  provides  the  best  heat  transfer  performance  until 
the  lowest  heat  flux  region,  as  did  Lake's  data  [Ref.  3]. 
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Figure  5 . 7  Comparison  of  Smooth  Tube  P/D  Ratios  for 
Increasing  Heat  Flux  with  the  Lower  Tube  Unheated 
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b.    Smooth    Tubes   -with   LTHFS   of   1   kW/m2 

Figure  5.9  shows  the  increasing  heat  flux  data 
taken  for  the  five  tube  configurations  (C1-C5)  with  a  heat 
flux  of  1  kW/m2  on  the  lower  tube.  Both  the  upper  and  lower 
tubes  for  all  configurations  started  in  the  natural  convection 
region.  At  low  heat  fluxes,  some  scatter  is  present  which 
becomes  smaller  as  the  data  approaches  the  incipient  boiling 
points  and  tends  towards  the  Churchill  and  Chu  correlation 
[Ref.  9],  as  seen  previously.  Note  that  the  data  for  a  P/D  of 
1.2  crosses  the  Churchill-Chu  line  and  appears  to  nucleate 
earlier  than  the  others.  For  this  small  spacing,  the  heated 
plume  from  the  lower  tube  can  actually  reduce  natural 
convection  from  the  top  tube  since  the  local  ambient  pool 
temperature  is  increased.  This  creates  a  larger  wall 
temperature  and  a  larger  wall  superheat  (Twall  -  Tsat)  .  After 
nucleation,  a  P/D  of  1.8  again  provides  better  heat  transfer 
in  the  boiling  region  than  the  other  P/D  ratios.  This  trend 
suggests  that  in  the  boiling  region  the  enhancement  provided 
by  a  P/D  of  1.8  is  due  to  the  tube  spacing  itself  and  not  from 
heating  the  lower  tube. 

Figure  5.10  is  the  decreasing  heat  flux  data  for  to 
Figure  5.9.  A  P/D  of  1.8  provides  the  best  heat  transfer 
performance  until  approximately  4  kW/m2,  which  is  where  the 
boiling  curves  seem  to  merge.  With  a  LTHFS  of  1  kW/m2,  it 
appears  to  force  the  decreasing  boiling  curves  to  merge  at  a 


65 


10' 


CM 


E10' 


S.  10 

a. 
13 


10' 


+ 

P/D 

2.0 

X 

P/D 

1.8 

0 

P/D 

1.5 

a 

P/D 

1.2 

A 

P/D 

1.2  (rotated 

30 

deg) 

-Churchill-Chu 

<—  Incipient  Boiling 


xfi- 

x£}- 


10' 


10w  10' 

Wall  Superheat/(K) 


10' 


Figure  5 . 9  Comparison  of  Smooth  Tube  P/D  Ratios  for 
Increasing  Heat  Flux  with  1  JcW/m*  on  the  Lower  Tube 


66 


10' 


CM 


CO 
CD 


CO 

-2 


2.103 

CL 


10' 


E  10*     - 


-    P/D   2.0 

*9- 

*    P/D    1.8 

o    P/D    1.5 

*9- 

a    P/D    1.2 

a    P/D    1.2 

(rotated 

30   deg) 

§ 


A 

a 


? 


X 


<k 


Mill 


i        I 


10" 


10*  10' 

Wall  Superheat/(K) 


10' 


Figure  5 . 10   Comparison  of  Smooth  Tube  P/D  Ratios  for 
Decreasing  Heat  Flux  with  1  kW/m2"  on  the  Lower  Tube 


67 


higher  heat  flux  (4  kW/m2)  than  with  the  lower  tube  unheated 
(<  1  kW/m2)  . 

c.    Smooth    Tubes   with   LTHFS   of  5  kW/m2 

Figure  5.11  shows  the  increasing  heat  flux  data 
taken  for  the  five  tube  configurations  (C1-C5)  with  5  kW/m2  on 
the  lower  tube.  Both  upper  and  lower  tubes  for  all  the 
configurations  started  in  the  natural  convection  region.  At 
low  heat  fluxes,  the  data  is  to  the  left  of  the  Churchill  and 
Chu  correlation  [Ref.  9],  as  seen  in  Figures  5.7  and  5.9.  As 
the  heat  flux  is  increased,  some  interesting  things  happen. 
The  data  for  a  P/D  of  1.2  (rotated  30°)  quickly  crosses  the 
Churchill-Chu  line  and  begins  to  parallel  it,  then  nucleates 
at  a  low  heat  flux  of  about  3  kW/m2  .  This  is  believed  to  be 
an  anomaly  and  not  an  effect.  Of  the  remaining  P/D  ratios, 
1.8  is  the  only  one  to  tend  toward  the  Churchill-Chu  line  with 
increasing  heat  flux.  The  P/D  ratios  of  2.0,  1.5  and  1.2  seem 
to  parallel  the  Churchill-Chu  line,  probably  due  to  the  LTHFS 
intensifying  the  forced  convection  around  the  upper  tube.  In 
all  cases,  the  lower  tube  remained  in  natural  convection  until 
the  upper  tube  nucleated.  Once  the  upper  tube  nucleated,  the 
lower  tube  also  nucleated.  In  the  boiling  region,  a  P/D  of 
1.8  again  gave  the  best  enhancement. 

Figure  5.12  is  the  decreasing  heat  flux  data  for 
Figure  5.11.  Above  10  kW/m2,  P/D  of  1.8  is  clearly  providing 
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Figure  5 . 11  Comparison  of  Smooth  Tube  P/D  Ratios  for 
Increasing  Heat  Flux  with  5  kW/m2  on  the  Lower  Tube 
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Figure  5.12   Comparison  of  Smooth  Tube  P/D  Ratios  for 
Decreasing  Heat  Flux  with  5  kW/m2  on  the  Lower  Tube 
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the  best  heat  transfer.  Below  10  kW/m2,  the  upper  tube  is 
returning  to  natural  convection,  but  the  lower  tube  is  still 
nucleating.  The  bubbles  departing  the  lower  tube  impinge  upon 
the  upper  tube  causing  a  significant  increase  in  heat  transfer 
performance,  from  bubble  sweeping,  compared  to  the  data  for  a 
LTHFS  of  0  and  1  kW/m2 .  Also,  a  LTHFS  of  5  kW/m2  eliminates 
any  upper  tube  hysteresis  effect. 

d.  Smooth  Tubes  with   LTHFS  of  10  kW/m2 

Figure  5.13  shows  the  increasing  heat  flux  data 
taken  for  the  five  tube  configurations  (C1-C5)  with  10  kW/m2 
on  the  lower  tube.  For  this  LTHFS,  the  lower  tube  nucleated 
at  the  start  of  the  run  for  a  P/D  of  1.5,  but  the  upper  tube 
remained  in  natural  convection.  The  effect  of  the  bubbles 
from  the  lower  tube  impacting  and  sweeping  over  the  upper  tube 
is  seen  to  significantly  enhance  the  heat  transfer 
performance  of  the  upper  tube  in  the  low  heat  flux  regions. 
For  the  other  P/D  ratios,  the  incipient  boiling  points  were 
very  random;  at  the  largest  spacing,  nucleation  occurred 
quickly,  whereas  at  the  smallest  spacing,  nucleation  was 
delayed.  Possibly,  the  largest  spacing  (P/D  of  2.0)  allowed 
the  heated  plume  from  the  lower  tube  to  thin  the  boundary 
layer  around  the  upper  tube,  but  the  closest  spacing  (P/D  of 
1.2)  enabled  the  heated  plume  to  'skirt'  around  the  upper  tube 
and  not  disturb  the  upper  tube  boundary  layer  as  much, 
delaying  the  point  of  incipience.  Again,  once  the  upper  tube 
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nucleated,  the  lower  tube  also  nucleated  shortly  thereafter. 
In  the  fully  developed  nucleate  boiling  region,  a  P/D  of  1.8 
still  gives  the  best  performance. 

Figure  5.14  is  the  decreasing  heat  flux  data  for 
Figure  5.13.  A  P/D  of  1.8  provided  the  best  heat  transfer 
until  the  heat  flux  was  lowered  to  about  15  kW/m2,  where  no 
gain  over  the  other  P/D  ratios  is  noticed.  This  concurs  with 
the  findings  of  Lake  [Ref.  3].  Below  10  kW/m2,  the  nucleating 
lower  tube  is  definitely  enhancing  the  upper  tube  performance 
for  all  P/D  ratios,  although,  the  P/D  of  1.2  (rotated  30°)  is 
slightly  to  the  right  of  the  other  curves.  This  is  probably 
caused  by  the  angular  offset  of  the  upper  tube,  not  allowing 
the  bubbles  from  the  lower  tube  to  directly  impinge  upon  the 
upper  tube.  Therefore,  a  30°  rotation  slightly  decreases  the 
performance  of  the  upper  tube  compared  to  an  in-line  tube 
configuration. 

e.  Smooth  Tubes  vith   LTHFS  of  20  kW/w2 

Figure  5.15  shows  the  increasing  heat  flux  data 
taken  for  the  five  tube  configurations  (C1-C5)  with  20  kW/m2 
on  the  lower  tube.  In  each  run,  the  lower  tube  was  nucleating 
at  the  beginning  of  the  run  and  remained  nucleating 
throughout.  The  only  effect  from  P/D  spacing  is  noticed  above 
20  kW/m2  where  a  P/D  of  1.8  again  provides  the  best 
performance . 
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Figure  5.14   Comparison  of  Smooth  Tube  P/D  Ratios  for 
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Figure  5 . 15   Comparison  of  Smooth  Tube  P/D  Ratios  for 
Increasing  Heat  Flux  with  20  kW/m2-  on  the  Lower  Tube 
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Figure  5.16  is  the  decreasing  heat  flux  data  for 
Figure  5.15.  Again,  below  20  kW/m2,  only  a  P/D  of  1.2 
(rotated  30°)  drops  slightly  below  the  other  curves  at  mid- 
range  heat  fluxes.  It  appears  that  this  LTHFS  enables  more 
bubbles  from  the  lower  tube  to  impinge  upon  the  upper  tube, 
compared  to  a  LTHFS  of  10  kW/m2 . 

f.    Smooth   Tubes   with   UTHFS  equal   LTHFS 

Figure  5.17  shows  the  increasing  heat  flux  data 
taken  for  the  five  tube  configurations  (C1-C5)  with  the  UTHFS 
equal  to  the  LTHFS  throughout  the  entire  length  of  the  run. 
At  the  beginning  of  all  the  runs,  both  tubes  were  in  the 
natural  convection  region.  In  the  natural  convection  region, 
a  P/D  of  1.2  (rotated  30°)  is  the  only  configuration  that  is 
very  close  to  the  Churchill  and  Chu  correlation  [Ref .  9]  . 
This  was  not  expected  nor  can  be  explained.  Again,  a  random 
incipient  point  was  obtained  and  in  the  boiling  region,  a  P/D 
of  1.8  was  the  best  performer. 

Figure  5.18  is  the  decreasing  heat  flux  data  for 
Figure  5.17.  Above  20  kW/m2,  a  P/D  of  1.8  still  provides  the 
best  performance  and  below  20  kW/m2  it  appears  that  as  the  P/D 
decreases,  so  does  the  enhancement  of  the  upper  tube.  Between 
9  and  5  kW/m2,  it's  believed  that  the  tubes  shift  from  the 
boiling  region  back  to  natural  convection.  This  may  cause  the 
decrease  of  enhancement  because,  as  speculated  earlier,  as  the 
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Figure  5.16  Comparison  of  Smooth  Tube  P/D  Ratios  for 
Decreasing  Heat  Flux  with  20  JcW/m*  on  the  Lower  Tube 
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P/D  is  decreased  so  too  are  the  convection  effects  from  the 
lower  tube. 

g.    TUBBO-B   Tubes   with   Lower   Tube   Unheated 

Figure  5.19  shows  the  decreasing  heat  flux  data 
taken  for  the  five  tube  configurations  (C1-C5)  with  lower  tube 
unheated.  At  the  start  of  each  run,  an  undetermined  amount  of 
nucleation  was  occurring  from  the  lower  tube  due  to  the 
vigorous  boiling  action  off  the  upper  tube  at  the  run  start 
heat  flux  (100  kW/m2)  .  A  very  consistent  boiling  curve  above 
10  kW/m2  indicates  that  there  is  no  enhancement  from  the  lower 
tube  at  any  P/D  ratio.  Below  5  kW/m2,  the  rotated 
configuration  data  agrees  well  with  the  single  tube  data  in 
Figure  5.3.  This  indicates  that  the  lower  tube  is  not 
modifying  the  flow  around  the  upper  tube  and  the  upper  tube  is 
performing  as  a  single  tube  which  provides  less  heat  transfer 
than  the  in-line  configurations. 

h.  TUBBO-B    Tubes   with   LTHFS  of  1   kW/m2 

Figure  5.20  shows  the  decreasing  heat  flux  data 
taken  for  the  five  tube  configurations  (C1-C5)  with  1  kW/m2  on 
the  lower  tube.  At  the  start  of  each  run,  both  the  upper  and 
lower  tubes  were  nucleating.  Very  little  change  is  noticed 
between  Figures  5.19  and  5.20.  Above  20  kW/m2,  the  plots  are 
identical.  Below  10  kW/m2,  the  LTHFS  of  1  kW/m2  moved  the  data 
curves  slightly  to  the  left.  This  indicates  that  the  greater 
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Figure  5.19   Comparison  of  TURBO- B  Tube  P/D  Ratios  for 
Decreasing  Heat  Flux  with  the  Lower  Tube  Unheated 
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the  LTHFS,  the  more  enhanced  the  mid-range  and  lower  heat  flux 
region  becomes.  As  seen  in  the  previous  figure,  a  P/D  of  1.2 
(rotated  30°)  is  again  less  enhanced  than  the  other  P/D 
ratios . 

i.  TUKBO-B    Tubes    with   LTHFS   of  5   kW/m2 

Figure  5.21  shows  the  decreasing  heat  flux  data 
taken  for  the  five  tube  configurations  (C1-C5)  with  5  kW/m2  on 
the  lower  tube.  At  the  start  of  each  run,  both  tubes  were 
vigorously  boiling.  Above  20  kW/m2,  the  data  is  consistent 
with  the  two  previous  plots,  but  below  10  kW/m2,  a  definite 
separation  between  a  P/D  of  1.2  (rotated  30°)  and  the  other 
P/D  ratios  is  seen.  Also  for  heat  fluxes  <  10  kW/m2,  the  in- 
line configurations  provide  better  heat  transfer  than  the  same 
configurations  seen  in  Figure  5.20.  This  indicates  that  as 
LTHFS  increases,  more  bubbles  depart  the  lower  tube  and 
increase  the  upper  tube  heat  transfer  performance  for  the  in- 
line configurations.  For  the  rotated  configuration,  the 
bubbles  must  be  departing  the  lower  TURBO-B  tube  differently 
(compare  to  a  lower  smooth  tube)  and  are  not  impinging  upon 
the  upper  TURBO-B  tube.  This  explains  the  approximately  40  - 
45%  decrease  in  heat  transfer  performance  for  the  rotated 
configuration  compared  to  the  in-line  configurations. 

j.    TURBO-B   Tubes   with   LTHFS  of  10  kW/m2 

Figure  5.22  shows  the  decreasing  heat  flux  data 
taken  for  the  five  tube  configurations  (C1-C5)  with  10  kW/m2 
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Figure  5.21   Comparison  of  TURBO-B  Tube  P/D  Ratios  for 
Decreasing  Heat  Flux  with  5  JcW/ra2  on  the  Lower  Tube 
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Figure  5.22  Comparison  of  TURBO-B  Tube  P/D  Ratios  for 
Decreasing  Heat  Flux  with  10  JcW/m2"  on  the  Lower  Tube 
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on  the  lower  tube.  Again,  both  tubes  were  nucleating  at  the 
beginning  of  each  run.  There  are  no  obvious  changes  above  20 
kW/m2  from  the  previous  plots.  Below  10  kW/m2,  all  the  curves 
seem  to  have  shifted  to  the  left  compared  to  the  previous 
case . 

k.    TUKBO-B   Tubes  with   LTHFS   of  20   kW/m2 

Figure  5.23  shows  the  decreasing  heat  flux  data 
taken  for  the  five  tube  configurations  (C1-C5)  with  20  kW/m2 
on  the  lower  tube.  Figures  5.22  and  5.23  show  similar 
results,  therefore,  increasing  the  LTHFS  from  10  to  20  kW/m2 
provides  a  small  change  in  upper  tube  heat  transfer 
performance . 

I.  TUKBO-B   Tubes   with   UTHFS  equal   LTHFS 

Figure  5.24  shows  the  decreasing  heat  flux  data 
taken  for  the  five  tube  configurations  (C1-C5)  with  the  UTHFS 
equal  to  the  LTHFS.   There  is  not  much  change  in  this  plot 
relative  to  the  trend  set  by  the  previous  plots. 
2 .   The  Effect  of  a  Fixed  Lower  Tube  Heat  Flux 

a.  Smooth  Tubes  with  P/D  of  2.0 

Figure  5.25  shows  the  increasing  heat  flux  data 
taken  for  a  P/D  of  2.0  at  five  different  heat  flux  settings  on 
the  lower  tube.  The  Churchill  and  Chu  correlation  [Ref .  9]  is 
also  shown  for  comparison.  For  a  LTHFS  of  20  kW/m2,  the  lower 
tube  had  nucleated  immediately  at  the  beginning  of  the  run  and 
remained  nucleating  for  the  length  of  the  run.   All  other 
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Figure  5.23   Comparison  of  TURBO-B  Tube  P/D  Ratios  for 
Decreasing  Heat  Flux  with  20  JcW/m1  on  the  Lower  Tube 
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Figure  5.24   Comparison  of  TURBO-B  Tube  P/D  Ratios  for 
Decreasing  Heat  Flux  with  the  UTHFS  equal  to  the  LTHFS 


88 


10s 


eg 

E 


cd 


104     -J 


- 


2-10°   _| 


10' 


q" 
q" 

q" 

q" 

q" 


lower  tube  -  0  kW/m*2 

lower  tube  -   1  kW/m*2 

lower  tube  -  5   kW/m*2 

lower  tube  -   10  kW/m*2 

lower  tube  -  20  kW/m*2 


Churchill-Chu 


< —  Incipient   Boiling 


x  0+ 


I     I    I  I  l 


Mill 


10' 


10w  10 

Wall  Superheat/(K) 


10' 


Figure  5.25   Comparison  of  Lower  Tube  Heat  Flux  Settings  on 
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LTHFS  data  runs,  started  in  natural  convection.  In  the 
natural  convection  region,  the  data  are  shifted  to  the  left  of 
the  Churchill  and  Chu  correlation  [Ref.  9],  but  approach  it  as 
heat  flux  is  increased.  This  is  probably  caused  by  strong 
circulation  patterns  within  the  small  evaporator  enhancing  the 
heat  transfer  mentioned  previously.  In  natural  convection, 
the  effect  of  a  LTHFS  had  little  influence  on  the  upper  tube 
heat  transfer  behavior;  in  fact,  this  is  repeated  for  all  the 
P/D  ratios.  The  point  of  incipience  was  very  repeatable,  with 
the  exception  of  a  LTHFS  of  10  kW/m2 .  A  heat  flux  of  10  kW/m2 
is  approximately  the  heat  flux  that  'trips'  the  nucleation 
sites  causing  boiling  to  occur  from  a  smooth  tube,  so  shortly 
after  the  run  started,  the  lower  tube  nucleated  and  caused  the 
upper  tube  to  nucleate  quickly.  At  high  heat  fluxes,  all  the 
data  falls  on  one  nucleate  boiling  curve  indicating  that  the 
LTHFS  has  little  effect  on  the  heat  transfer  from  the  upper 
tube  in  this  region. 

Figure  5.26  shows  the  decreasing  heat  flux  data  for 
the  data  run  for  Figure  5.25.  Important  to  note  that  the 
upper  tube  starts  and  finishes  each  run  in  natural  convection, 
but  the  lower  tube  might  not.  For  a  LTHFS  of  0  and  1  kW/m2, 
the  lower  tube  starts  and  finishes  in  natural  convection;  for 
a  LTHFS  of  20  kW/m2,  it  starts  and  finishes  in  nucleate 
boiling;  and  for  a  LTHFS  of  5  and  10  kW/m2,  it  starts  in 
natural  convection  and  finishes  in  nucleate  boiling.  At  high 
heat  fluxes,  all  the  data  lay  on  one  nucleate  boiling  curve  as 
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Figure  5.26  Comparison  of  Lower  Tube  Heat  Flux  Settings  on 
Decreasing  Heat  Flux  of  Smooth  Tubes  with  P/D  of  2.0 
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expected.  Below  20  kW/m2,  five  separate  curves  are  seen.  The 
curves  with  a  LTHFS  of  0  and  1  kW/m2  return  to  the  natural 
convection  region,  indicating  a  deactivation  of  the  nucleation 
sites  as  the  upper  tube  heat  flux  is  decreased,  showing  the 
hysteresis  effect  of  the  upper  tube.  At  a  heat  flux  below  5 
kW/m2,  the  curves  with  a  LTHFS  of  5,  10  and  20  kW/m2  begin  to 
parallel  the  natural  convection  curve  found  on  the  increasing 
heat  flux  plot  and  have  significant  enhancement  over  the  0  and 
1  kW/m2  curves  at  low  heat  fluxes.  Also,  the  hysteresis 
effect  is  eliminated  by  the  vigorously  nucleating  lower  tube. 
This  same  phenomena  was  found  by  Lake  [Ref .  3]  .  He  attributed 
this  enhancement  to  strong  convection  effects  on  the  upper 
tube  from  the  bubbles  nucleating  off  the  lower  tube. 
Therefore,  significant  enhancement  of  heat  transfer  of  an 
upper  tube  is  possible  when  the  lower  tube  is  nucleating  and 
very  little  enhancement  is  noticed  when  both  tubes  are  in 
natural  convection. 

Jb.  Smooth   Tubes   with  P/D  of  1.8 

Figure  5.27  shows  the  increasing  heat  flux  data 
taken  for  a  P/D  of  1.8  at  five  lower  tube  heat  flux  settings. 
For  a  LTHFS  of  20  kW/m2,  the  lower  tube  had  nucleated 
immediately  at  the  beginning  of  the  run  and  remained 
nucleating  for  the  entire  length  of  the  run.  All  other  LTHFS 
started  in  natural  convection.  In  the  natural  convection 
region,  the  data  are  shifted  to  the  left  of  the  Churchill  and 
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Figure  5 . 27   Comparison  of  Lower  Tube  Heat  Flux  Settings  on 
Increasing  Heat  Flux  of  Smooth  Tubes  with  P/D  of  1.8 
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Chu  correlation  [Ref.  9],  but  approach  it  as  heat  flux  is 
increased,  as  was  the  case  with  a  P/D  of  2.0.  Again,  the 
effect  of  the  LTHFS  had  little  influence  of  the  upper  tube 
heat  transfer  behavior.  The  point  of  incipience  was  again, 
very  repeatable.  This  time,  unlike  for  a  P/D  of  2.0,  the 
lower  tube  remained  in  natural  convection  until  the  upper  tube 
nucleated.  Again,  at  high  heat  fluxes,  all  the  data  falls  on 
one  nucleate  boiling  curve  indicating  that  the  LTHFS  has 
little  effect  upon  the  heat  transfer  in  this  region. 

Figure  5.28  shows  the  decreasing  heat  flux  data  for 
the  data  runs  for  Figure  5.27.  The  data  trend  is  similar  to 
that  found  in  Figure  5.26.  Therefore,  significant  enhancement 
of  heat  transfer  from  an  upper  tube  is  again  seen  when  the 
lower  tube  is  nucleating  with  larger  than  5  kW/m2  on  the  lower 
tube. 

c.  Smooth  Tubes  vlth  P/D  of  1.5 

Figure  5.29  shows  the  increasing  heat  flux  data 
taken  for  a  P/D  of  1.5  at  five  lower  tube  heat  flux  settings. 
For  a  LTHFS  of  10  and  20  kW/m2,  the  lower  tube  had  nucleated 
immediately  at  the  beginning  of  the  run  and  remained 
nucleating  for  the  length  of  the  run.  All  other  LTHFS  started 
in  natural  convection.  This  plot  shows  similar  results  seen 
on  Figure  5.27  and  most  of  the  comments  still  apply. 
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Figure  5 . 28   Comparison  of  Lower  Tube  Heat  Flux  Settings  on 
Decreasing  Heat  Flux  of  Smooth  Tubes  with  P/D  of  1 . 8 
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Figure  5.29   Comparison  of  Lower  Tube  Heat  Flux  Settings  on 
Increasing  Heat  Flux  of  Smooth  Tubes  with  P/D  of  1.5 


96 


Figure  5.30  shows  the  decreasing  heat  flux  data  for 
the  data  runs  shown  in  Figure  5.29.  Again,  similar  results 
are  obtained  as  before.  Significant  enhancement  in  the  heat 
transfer  from  an  upper  tube  is  obtained  when  the  lower  tube  is 
nucleating  with  a  heat  flux  larger  than  5  kW/m2 . 

d.  Smooth   Tubes   vith  P/D  of  1.2 

Figure  5.31  shows  the  increasing  heat  flux  data 
taken  for  a  P/D  of  1.2  at  five  lower  tube  heat  flux  settings. 
For  a  LTHFS  of  20  kW/m2,  the  lower  tube  had  nucleated 
immediately  at  the  beginning  of  the  run  and  remained 
nucleating  for  the  length  of  the  run.  All  other  LTHFS  started 
in  natural  convection.  This  plot  shows  similar  results  seen 
on  Figure  5.29  and  most  of  the  comments  still  apply. 

Figure  5.32  shows  the  decreasing  heat  flux  data  for 
the  data  runs  shown  in  Figure  5.31.  This  plot  shows  similar 
results  seen  on  Figure  5.30  and  again,  most  of  the  comments 
still  apply. 

e.  Smooth   Tubes   vith  P/D  of  1.2    (Rotated  30°) 
Figure  5.33  shows  the  increasing  heat  flux  data 

taken  for  a  P/D  of  1.2  (rotated  30°),  at  five  lower  tube  heat 
flux  settings.  For  a  LTHFS  of  20  kW/m2,  the  lower  tube  had 
nucleated  immediately  at  the  beginning  of  the  run  and  remained 
nucleating  for  the  length  of  the  run.  All  other  LTHFS  started 
in  natural  convection.  The  increasing  data  for  this 
configuration  is  more  concentrated  than  the  others  seen 
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Figure  5.30   Comparison  of  Lower  Tube  Heat  Flux  Settings  on 
Decreasing  Heat  Flux  of  Smooth  Tubes  with  P/D  of  1.5 
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Figure  5.31  Comparison  of  Lower  Tube  Heat  Flux  Settings  on 
Increasing  Heat  Flux  of  Smooth  Tubes  with  P/D  of  1.2 


99 


10' 


CM 


*E     10"     - 


X 

3 


at 
© 


© 


© 


10°     - 


10' 


-t- 

q" 

lower 

tube  -  0   kW/m"2 

X 

q" 

lower 

tube  -  1  kW/m"2 

0 

q" 

lower 

tube  -  5   kW/m"2 

a 

q" 

lower 

tube  -  10  kW/rrf  2 

A 

q" 

lower 

tube  -  20  kW/m*2 

& 


AQ  + 


A  □ 


a 


aQ 


0 

X 


A 

a 


10' 


10w  10 

Wall  Superheat/(K) 


10' 


Figure  5.32  Comparison  of  Lower  Tube  Heat  Flux  Settings  on 
Decreasing  Heat  Flux  of  Smooth  Tubes  with  P/D  of  1.2 
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Increasing  Heat  Flux  of  Smooth  Tubes  with  P/D  of  1.2 
(Rotated  30°) 
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previously.   This  might  be  due  to  the  offset  orientation  of 
the  upper  tube  causing  the  data  to  be  less  scattered. 

Figure  5.34  shows  the  decreasing  heat  flux  data  for 
the  data  runs  shown  in  Figure  5.33.  At  high  heat  fluxes,  all 
the  data  lay  on  one  nucleate  boiling  curve.  Below  10  kW/m2, 
separate  curves  are  formed  and  as  seen  before.  Significant 
enhancement  of  heat  transfer  from  an  upper  tube  is  obtained 
when  the  lower  tube  is  nucleating  with  a  heat  flux  larger  than 
5  kW/m2. 

f.  TURBO-B   Tubes   with  P/D  of  2.0 

Figure  5.35  shows  the  decreasing  heat  flux  data  for 
the  TURBO-B  tubes  taken  at  five  lower  tube  heat  flux  settings. 
At  the  start  of  each  run,  a  very  consistent  boiling  curve  is 
seen  above  20  kW/m2,  indicating  that  above  a  UTHFS  of  20  kW/m2, 
there  is  no  enhancement  from  increasing  the  LTHFS.  Below  10 
kW/m2,  raising  the  LTHFS  does  increase  the  heat  transfer 
performance  of  the  upper  tube.  At  an  upper  tube  heat  flux  of 
2  kW/m2,  an  increase  in  performance  of  2.4  is  found  going  from 
an  unheated  lower  tube  to  a  LTHFS  of  20  kW/m2,  compared  to  an 
increase  of  3.6  for  the  smooth  tube  array  (Figure  5.26). 

g.  TURBO-B   Tubes   with  P/D  of  1.8 

Figure  5.36  shows  the  decreasing  heat  flux  data  for 
the  TURBO-B  tubes  taken  at  five  lower  tube  heat  flux  settings. 
Similar  results  are  found  compared  to  the  P/D  of  2.0  case.  At 
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Figure  5.34   Comparison  of  Lower  Tube  Heat  Flux  Settings  on 
Decreasing  Heat  Flux  of  Smooth  Tubes  with  P/D  of  1 . 2 
(Rotated  30°) 
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Figure  5.35   Comparison  of  Lower  Tube  Heat  Flux  Settings  on 
Decreasing  Heat  Flux  of  TURBO-B  Tubes  with  P/D  of  2 . 0 
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Figure  5.36  Comparison  of  Lower  Tube  Heat  Flux  Settings  on 
Decreasing  Heat  Flux  of  TURBO-B  Tubes  with  P/D  of  1 . 8 
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an  upper  tube  heat  flux  of  2  kW/m2,  an  increase  in  performance 
of  1.6  is  found  going  from  a  unheated  lower  tube  to  a  LTHFS  of 
20  kW/m2,  compared  to  an  increase  of  3.7  for  the  smooth  tube 
array  (Figure  5.28)  .  Unlike  the  smooth  tube  case  where  a  P/D 
of  1.8  produced  the  best  heat  transfer  performance  in  nucleate 
boiling,  no  optimal  P/D  ratio  was  found  for  the  TURBO-B  tube 
array. 

h.    TURBO-B   Tubes   with  P/D  of  1.5 

Figure  5.37  shows  the  decreasing  heat  flux  data  for 
the  TURBO-B  tubes  taken  at  five  lower  tube  heat  flux  settings. 
Similar  results  are  found  compared  to  the  P/D  of  1.8  case.  At 
an  upper  tube  heat  flux  of  2  kW/m2,  an  increase  in  performance 
of  2.2  is  found  going  from  a  unheated  lower  tube  to  a  LTHFS  of 
20  kW/m2,  compared  to  an  increase  of  3.5  for  the  smooth  tube 
array  (Figure  5.30). 

i.  TUKBO-B   Tubes   vith  P/D  of  1.2 

Figure  5.38  shows  the  decreasing  heat  flux  data  for 
the  TURBO-B  tubes  taken  at  five  lower  tube  heat  flux  settings. 
Similar  results  are  found  compared  to  the  P/D  of  1.5  case.  At 
an  upper  tube  heat  flux  of  2  kW/m2,  an  increase  in  performance 
of  2.2  is  found  going  from  a  unheated  lower  tube  to  a  LTHFS  of 
20  kW/m2,  compared  to  an  increase  of  5.4  for  the  smooth  tube 
array  (Figure  5.32). 
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Figure  5.37  Comparison  of  Lower  Tube  Heat  Flux  Settings  on 
Decreasing  Heat  Flux  of  TURBO-B  Tubes  with  P/D  of  1.5 
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Figure  5.38   Comparxson  of  Lower  Tube  Heat  Flux  Settings  on 
Decreasing  Heat  Flux  of  TURBO-B  Tubes  with  P/D  of  1.2 
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j.     TURBOS   Tubes   vith  P/D  of  1.2    (Rotated   30°) 

Figure  5.39  shows  the  decreasing  heat  flux  data  for 
the  TURBO-B  tubes  taken  at  five  lower  tube  heat  flux  settings. 
A  very  consistent  boiling  curve  is  again  seen  above  10  kW/m2 . 
Below  10  kW/m2,  the  data  seems  less  scattered  and  slightly  to 
the  left  of  any  other  configuration,  indicating  that  rotating 
the  upper  tube  reduces  the  enhancement  produced  from 
increasing  the  LTHFS  in  the  low  heat  flux  regions.  At  an 
upper  tube  heat  flux  of  2  kW/m2,  an  increase  in  performance  of 
1.7  is  found  going  from  a  unheated  lower  tube  to  a  LTHFS  of  20 
kW/m2,  compared  to  an  increase  of  3.5  for  the  smooth  tube 
array  (Figure  5.34). 

3.   The  Effect  of  a  Variable  Lower  Tube  Heat  Flux 
a.  Smooth   Tubes   with  P/D  of  2.0 

Figure  5.40  shows  the  heat  transfer  performance  of 
an  upper  tube,  held  at  a  constant  heat  flux,  as  the  lower  tube 
heat  flux  is  varied.  The  three  fixed  UTHFS  used  were  5,  10 
and  20  kW/m2.  The  labeling  in  the  plot  legend  describes  which 
data  were  taken  using  increasing  or  decreasing  lower  tube  heat 
flux.  The  sharp  transitions  in  the  increasing  curves  are  the 
points  where  the  lower  tube  nucleated  and  caused  the  upper 
tube  to  shift  from  natural  convection  to  nucleate  boiling 
(visually  verified  during  the  data  run) ,  significantly 
increasing  the  upper  tube  heat  transfer  coefficient,  as 
expected.   The  maximum  upper  tube  heat  transfer  coefficient 
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Figure  5.39   Comparison  of  Lower  Tube  Heat  Flux  Setting  on 
Decreasing  Heat  Flux  of  TURBO- B  Tubes  with  P/D  of  1.2 
(Rotated  30°) 
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Tube  at  a  Constant  Heat  Flux  for  a  P/D  of  2 . 0 
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for  all  three  UTHFS  was  1850  W/m2-K,  which  indicates  that  at 
a  lower  tube  heat  flux  of  80  kW/m2,  the  heat  transfer 
performance  of  the  upper  tube  is  independent  of  the  lower  tube 
heat  flux.  The  decreasing  curves  show  that  for  an  UTHFS  of  20 
kW/m2,  the  lower  tube  heat  flux  has  little  influence  on  the 
upper  tube  heat  transfer.  For  the  UTHFS  of  10  kW/m2,  as  the 
lower  tube  heat  flux  is  decreased,  the  upper  tube  remains  in 
nucleate  boiling,  but  for  the  UTHFS  of  5  kW/m2,  the  upper  tube 
returns  to  the  natural  convection  region  as  the  lower  tube 
heat  flux  is  decreased. 

b.    Smooth   Tubes   vith  P/D  of  1.8 

Figure  5.41  shows  the  heat  transfer  performance  of 
an  upper  tube,  held  at  a  constant  heat  flux,  as  the  lower  tube 
heat  flux  is  varied.  Similar  results  were  obtained  as  in  the 
P/D  of  2.0  case,  except  that  the  curves  did  not  merge  at  the 
highest  heat  flux  setting.  The  maximum  heat  transfer 
coefficient  for  the  UTHFS  of  5  and  10  kW/m2  converged  at 
approximately  1700  W/m2-K  and  for  the  UTHFS  of  20  kW/m2  the 
value  was  about  1900  W/m2-K.  It  is  not  known  why  the  curves 
did  not  merge  at  high  lower  tube  heat  fluxed.  Comparing  these 
results  to  the  P/D  of  2.0  results,  at  the  smallest  lower  tube 
heat  fluxes,  the  P/D  of  1.8  provides  better  heat  transfer. 
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Figure  5.41   Heat  Transfer  Performance  of  a  Smooth  Upper 
Tube  at  a  Constant  Heat  Flux  for  a  P/D  of  1 . 8 
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Figure  5 . 42   Heat  Transfer  Performance  of  a  Smooth  Upper 
Tube  at  a  Constant  Heat  Flux  for  a  P/D  of  1.5 
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c.  Smooth  Tubes  with  P/D  of  1.5 

Figure  5.42  shows  the  heat  transfer  performance  of 
an  upper  tube,  held  at  a  constant  heat  flux,  as  the  lower  tube 
heat  flux  is  varied.  The  maximum  heat  transfer  coefficient 
for  all  three  of  the  UTHFS  converged  at  approximately  1800 
W/m2*K,  similar  to  the  P/D  of  2.0  case.  The  decreasing  curves 
show  very  similar  behavior  found  at  a  P/D  of  2.0.  Comparing 
these  results  to  data  of  a  P/D  of  2.0  and  1.8,  a  P/D  of  2.0 
and  1.5  are  similar  in  the  amount  of  enhancement  achieved, 
whereas  a  P/D  of  1.8  outperforms  them  both,  in  the  mid  to  low 
heat  flux  regions. 

d.  Smooth  Tubes  with  P/D  of  1.2 

Figure  5.43  shows  the  heat  transfer  performance  of 
an  upper  tube,  held  at  a  constant  heat  flux,  as  the  lower  tube 
heat  flux  is  increased.  The  maximum  heat  transfer  coefficient 
for  all  three  of  the  UTHFS  converged,  as  did  the  P/D  of  2.0 
and  1.5,  at  approximately  1990  W/m2*K,  highest  value  of  all 
the  configurations  tested.  Therefore,  the  best  heat  transfer 
performance  was  given  by  a  P/D  ratio  of  1.8  when  the  lower 
tube  heat  flux  was  held  constant,  but  a  P/D  ratio  of  1.2  when 
the  upper  tube  heat  flux  was  held  constant.  No  decreasing 
curves  were  recorded  for  these  data  runs. 

e.  Smooth  Tubes  with  P/D  of  1.2    (Rotated  30°) 
Figure  5.44  shows  the  heat  transfer  performance  of 

an  upper  tube,  held  at  a  constant  heat  flux,  as  the  lower  tube 
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Figure  5 . 43   Heat  Transfer  Performance  of  a  Smooth  Upper 
Tube  at  a  Constant  Heat  Flux  for  a  P/D  of  1.2 
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Figure  5.44   Heat  Transfer  Performance  of  a  Smooth  Upper 
Tube  at  a  Constant  Heat  Flux  for  a  P/D  of  1.2  (Rotated  30°) 
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heat  flux  is  increased.  As  seen  for  the  P/D  of  1.8  case,  the 
curves  do  not  merge  at  the  highest  lower  tube  heat  flux.  The 
maximum  heat  transfer  coefficient  for  the  two  lower  curves 
converged  at  about  1500  W/m2-K  and  the  upper  curve  was  1725 
W/m2,K.  Comparing  Figures  5.43  and  5.44,  for  an  UTHFS  of  20 
kW/m2,  the  data  curves  are  identical  for  lower  tube  heat 
fluxes  up  to  8  kW/m2,  indicating  that  the  forced  convection 
effects  from  the  lower  tube  have  very  little  influence  on  the 
upper  tube  heat  transfer  behavior.  Above  lower  tube  heat 
fluxes  of  10  kW/m2,  the  in-line  configuration  provides  better 
heat  transfer  performance  because  the  lower  tube  is  in 
nucleate  boiling  and  all  of  the  bubbles  are  impacting  the 
upper  tube  providing  enhancement  by  bubble  pumping;  as 
opposed  to  the  offset  configuration,  where  most  of  the  bubbles 
coming  off  the  lower  tube  miss  the  upper  tube.  The  same  idea 
holds  true  for  the  other  data  curves.  The  in-line 
configurations  are  more  enhanced  due  to  bubble  pumping  from 
the  lower  tube.  No  decreasing  curves  were  recorded  for  these 
data  runs. 

f.    TUBBO-B   Tubes   with  P/D  of  2.0 

Figure  5.45  shows  the  heat  transfer  performance  of 
an  upper  tube,  held  at  a  constant  heat  flux,  as  the  lower  tube 
heat  flux  is  decreased.  At  high  lower  tube  heat  fluxes,  all 
three  of  the  UTHFS  appear  to  perform  almost  equally.  For  an 
UTHFS  of  20  kW/m2,  the  lower  tube  heat  flux  has  little 
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Figure  5 . 45   Heat  Transfer  Performance  of  a  TURBO-B  Upper 
Tube  at  a  Constant  Heat  Flux  for  a  P/D  of  2.0 


119 


influence  on  the  upper  tube  heat  transfer  and  for  an  UTHFS  of 
5  and  10  kW/m2,  as  the  lower  tube  heat  flux  is  decreased  so 
follows  the  performance  of  the  upper  tube  (more  for  5  than  10 
kW/m2)  . 

g.    TUPBO-B   Tubes   vith  P/D  of  1.8 

Figure  5.46  shows  the  heat  transfer  performance  of 
an  upper  tube,  held  at  a  constant  heat  flux,  as  the  lower  tube 
heat  flux  is  decreased.  Comparing  these  curves  to  the  ones 
for  a  P/D  of  2.0,  a  small  enhancement  in  performance  is  gained 
in  the  mid  and  low  heat  flux  regions  using  a  P/D  of  1.8,  which 
is  an  expected  result. 

h .  TUKBO-B   Tubes   vi  th  P/D  of  1.5 

Figure  5.47  shows  the  heat  transfer  performance  of 
an  upper  tube,  held  at  a  constant  heat  flux,  as  the  lower  tube 
heat  flux  is  decreased.  There  is  very  little  difference  in 
performance  between  a  P/D  of  1.5  and  1.8.  Comparing  the  P/D 
of  1.5  to  the  2.0  data,  a  small  decrease  in  performance  is 
noted  for  lower  tube  heat  fluxes  >  20  kW/m2  using  the  P/D  of 
1.5  configuration,  which  cannot  be  explained. 

i.    TUPBO-B   Tubes   vith  P/D  of  1.2 

Figure  5.48  shows  the  heat  transfer  performance  of 
an  upper  tube,  held  at  a  constant  heat  flux,  as  the  lower  tube 
heat  flux  is  decreased.  Similar  performance  is  noted  for  a 
P/D  ratio  of  1.2  and  1.5,  therefore,  the  comments  made 
previously  still  apply. 
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Figure  5.46  Heat  Transfer  Performance  of  a  TURBO-B  Upper 
Tube  at  a  Constant  Heat  Flux  for  a  P/D  of  1 . 8 
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Figure  5 . 47   Heat  Transfer  Performance  of  a  TURBO-B  Upper 
Tube  at  a  Constant  Heat  Flux  for  a  P/D  of  1.5 
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Figure  5.48   Heat  Transfer  Performance  of  a  TURBO-B  Upper 
Tube  at  a  Constant  Heat  Flux  for  a  P/D  of  1 . 2 
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j.     TURBO-B    Tubes    with  P/D  of  1.2    (Rotated  30°) 

Figure  5.49  shows  the  heat  transfer  performance  of 
an  upper  tube,  held  at  a  constant  heat  flux,  as  the  lower  tube 
heat  flux  is  decreased.  Comparing  Figures  5.48  and  5.49,  for 
an  UTHFS  of  20  kW/m2,  the  data  curves  are  almost  identical, 
indicating  that  the  lower  tube  orientation  and  heat  flux  have 
very  little  influence  on  the  upper  tube  heat  transfer 
behavior.  For  the  other  two  data  curves,  the  in-line 
configuration  provides  better  heat  transfer  performance 
because  of  bubble  pumping,  discussed  earlier.  For  the  offset 
configuration,  all  of  the  data  curves  do  not  merge  and  are 
almost  horizontal  lines,  indicating  there  is  no  enhancement  by 
bubble  pumping.  It  is  believed  that  most  of  the  bubbles 
coming  off  the  lower  tube  miss  the  upper  tube  because  the 
bubbles  depart  the  TURBO-B  tube  surface  at  the  upper  most 
section  of  the  tube,  which  is  different  than  how  bubbles 
depart  the  smooth  tube  surface.  This  phenomena  needs  to  be 
further  investigated.  No  decreasing  curves  were  recorded  for 
these  data  runs. 

C.   R-124/R-114  PERFORMANCE  COMPARISON 

Figure  5.50  compares  the  decreasing  heat  flux  data  for  R- 
124  found  in  this  study  with  the  R-114  data  of  Lake  [Ref.  3] 
for  smooth  tube  configurations  listed  in  the  plot  legend,  with 
the  lower  tube  unheated.   In  the  mid  to  high  heat  flux 
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Figure  5.49  Heat  Transfer  Performance  of  a  TURBO-B  Upper 
Tube  at  a  Constant  Heat  Flux  for  a  P/D  of  1.2  (Rotated  30°) 
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regions,  the  P/D  ratios  of  1.8  and  1.5  demonstrated  improved 
performance  with  the  R-124,  but  with  a  P/D  of  2.0,  little  or 
no  improvement  was  found.  In  the  low  heat  flux  region,  the 
uncertainty  is  too  large  to  make  any  definitive  conclusions. 
Figure  5.51  compares  the  decreasing  heat  flux  data  for  R- 
124  with  the  R-114  data  of  Lake  [Ref.  3]  for  the  smooth  tube 
configurations  listed  in  the  plot  legend,  with  10  kW/m2  on  the 
lower  tube.  In  the  high  heat  flux  region,  a  P/D  of  1.8  and 
1.5  again  provided  better  performance  using  R-124,  but  for  R- 
114,  a  P/D  of  2.0  displayed  slightly  better  performance.  In 
the  mid-range  heat  fluxes,  R-114  exhibited  improved 
performance  in  each  of  the  three  configurations.  Again,  in 
the  low  heat  flux  region,  there  was  too  much  uncertainty  in 
the  data  to  draw  any  viable  conclusions. 
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VI .    CONCLUSIONS 

The  interaction  effects  of  a  lower  heated  tube  on  pool 
boiling  of  R-124  from  an  upper  horizontal  tube  has  been 
studied  for  tube  two  types:  smooth  and  TURBO-B.  Based  on  the 
data  obtained,  the  following  conclusions  are  made: 


1.  In  natural  convection,  the  effect  of  a  lower  heated  tube 
on  the  heat  transfer  from  an  upper  tube  was  small. 

2.  For  a  smooth  tube  array,  a  vigorously  nucleating  lower 
tube  eliminated  upper  tube  hysteresis. 

3.  For  a  smooth  tube  array  in  nucleate  boiling,  P/D  of  1.8 
gave  the  best  upper  tube  heat  transfer  performance.  No 
optimal  P/D  was  found  for  the  TURBO-B  tube  array. 

4.  For  both  arrays  with  upper  tube  heat  fluxes  <  5  kW/m2, 
the  enhancing  effect  of  bubbles  from  a  lower  tube  was 
dramatic.  This  enhancement  increased  as  lower  tube  heat 
flux  increased.  When  upper  tube  heat  fluxes  were  >  20 
kW/m2,  all  enhancement  disappeared. 

5.  For  a  TURBO-B  tube  array,  a  30  degree  offset  of  the  upper 
tube  reduced  the  upper  tube  heat  transfer  performance 
(compared  to  the  in-line  configurations).  This  may 
indicate  bubbles  depart  TURBO-B  tubes  differently  than 
smooth  tubes. 

6.  For  a  smooth  tube  array  with  the  lower  tube  unheated  and 
an  upper  tube  heat  flux  of  >  3  kW/m2,  the  performance  of 
R-124  compared  to  R-114,  improved  as  the  P/D  ratio  was 
reduced.  With  a  nucleating  lower  tube  (at  10  kW/m2)  , 
again  the  performance  of  R-124  improved  as  the  P/D  ratio 
decreased,  but  only  for  upper  tube  heat  fluxes  of  >  40 
kW/m2 . 
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VII.    RECOMMENDATIONS 


1.  Conduct  further  studies  using  TURBO-B  tubes  to  verify 
that  bubbles  do  depart  the  TURBO-B  tubes  differently  than 
smooth  tubes. 

2.  Conduct  similar  tests  using  different  enhanced  tubes  to 
better  understand  the  influence  of  a  lower  heated 
enhanced  tube  on  the  heat  transfer  behavior  of  an  upper 
tube. 

3.  Modify  the  existing  evaporator  to  accept  a  larger  number 
of  tubes  (three  and  four  smooth  or  enhanced  tubes)  and 
conduct  tests  to  study  the  effect  of  multiple  heated 
lower  tubes  on  the  heat  transfer  performance  of  an  upper 
tube. 

4.  Conduct  similar  tests  using  multiple  tube  configurations 
with  R-124/oil  mixtures  to  determine  the  effect  of  oil  on 
the  heat  transfer  performance  of  an  upper  tube. 
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APPENDIX  A.   PROGRAM  SETUP JY 

Program  SETUPJY  was  used  to  evaluate  the  system 
performance  prior  to  commencing  a  data  collection  run.  This 
program  is  written  in  Hewlett-Packard  Basic  5.0  language  and 
enables  the  operator  to: 

1.  Monitor  the  coolant  sump  temperature. 

2.  Monitor  the  evaporator  liquid  pool  average  temperature. 

3.  Monitor  all  the  thermocouple  channel  output  temperatures. 

4.  Monitor  the  voltage,  current  and  power  supplied  to  the 
upper  tube  heater  (main  heater)  as  well  as  the  heater 
element  electrical  resistance. 

5.  Monitor  the  voltage,  current  and  power  supplied  to  the 
lower  tube  heater  (auxiliary  heater)  as  well  as  the 
heater  element  electrical  resistance. 
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101  PROGRAM:  SETUPJY 

20!  DATS:     MAY  15,  1993 

30!  PROGRAMMER:  LT  DEAN  SUGIYAMA 

40!  MODIFIED  BY  LANNIE  LAKE  JAN  22,  1992 

50!  MODIFIED  BY  GEORGE  D.  PERRY  10  APR  1993 

50!  MODIFIED  3Y  JOE  YUSICIAN  MAY  15,  1993 

70  COM  /Cc/  C(7) 

80  DATA  0.10086091, 25727. 94369,  -767345. 3295, 78025595. 31, -9247486589, 5. 97688E+ 

11,-2. 66192E+13 

90  DATA  3.94078E+14 

100  READ  C(») 

110  ON  KEY  1,15  GOTO  150 

120  PRINTER  13  1 

150  PRINT  USING  "2X, ""SELECT  OPTION""" 

160  PRINT  USING  "5X, ""0-MONITOR  SUMP""" 

170  PRINT  USING  "6X,  "."1=MONITOR  LIQUID""" 

180  PRINT  USING  "6X,  ""2=»CHECK  THERMOCOUPLES""" 

190  PRINT  USING  "6X, ""3-CHECX  MAIN  HEATER""" 

200  PRINT  USING  "6X, ""4-CHECX  AUX  HEATERS""" 

210  PRINT  USING  "6X,""5»EXIT  PROGRAM""" 

230  BEEP 

240  INPUT  Ido 

250  IF  Ido>5  THEN  Ido-5 

260  IF  ldo»0  THEN  340 

270  IF  Ido-1  THEN  520 

280  IF  Ido-2  THEN  690 

290  IF  Ido-3  THEN  870 

300  IF  Ido-4  THEN  370 

310  IF  Ido-5  THEN  1170 

340  PRINT 

350  PRINT  "SUMP  TEMPERATURE  DEG  C  " 

360  PRINT 

370  OUTPUT  709; "AR  AF19  AL19  VR5" 

380  OUTPUT  709; "AS  SA" 

390  Sub-0 

400  FOR  J-l  TO  5 

410  ENTER  709;E 

420  3um»Sum+S 

430  NEXT  J 

440  Save»Sum/5 

450  Temp»FNTvsv(Save) 

460  PRINT  USING  "4X,MDD.DD";Tamp 

470  BEEP 

480  PRINT 

490  WAIT  5 

500  GOTO  340 

5101 

520  PRINT 

530  PRINT  "LIQUID  TEMPERATURE  DBG  C» 

540  PRINT 

550  OUTPUT  709;"AR  AF16  AL18  VR5" 

560  3um>0 

570  FOR  1-1  TO  2 

580  OUTPUT  709;  "AS  3A." 

590  ENTER  709;S 

600  Sust^Suin+S 

610  NEXT  I 

820  Bave-Sum/2 

630  T«np"FNTvsv(Bave) 

640  PRINT  USING-  "4X,MDD.DD";Tamp 

6S0  HBXP 

660  WAIT  5 

670  GOTO   510 
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5801 

590  PRINT 

700  PRINT  "TEMPERATURE         DEG  C" 

710  OUTPUT  709; "AR  AFOO  AL19  VR5" 

720  FOR  I»l  TO  20 

730  OUTPUT  709; "AS  SA" 

740  Sum=0 

750  FOR  J«l  TO  5 

760  ENTER  709; S 

770  3um*3um+E 

780  NEXT  J 

790  Save=3um/5 

300  Temp=FNTv3v(Save) 

810  PRINT  TAB(3);I;TAB(15);Tamp 

320  NEXT  I 

830  3EEP 

840  WAIT  5 

850  GOTO  590 

360! 

870  PRINT 

880  OUTPUT  709; "AR  AF20  AL22  VR5" 

890  FOR  1=1  TO  3 

900  OUTPUT  709; "AS  3A" 

910  3um-0 

920  FOR  J-l  TO  5 

930  ENTER  709;E 

940  Sum»Sum+E 

950  NEXT  J 

960  IT  I«l  THEN  Volt-Sum/ 5 

970  IF  1-2  AND  Ido»3  THEN 

980  PRINT  "MAKE  SURE  VOLTAGE  BOX  IS  SET  TO  MAIN  HEATERS" 

990  Amp-Sum/  5 

1000  END  IF 

1010  IF  I»3  AND  Ido-4  THEN 

1020  PRINT  "MAKE  3URS  VOLTAGE  BOX  IS  SET  TO  AUX  HEATERS" 

1030  Amp-Sun/ 5 

1040  END  IF 

1050  NEXT  I 

1060  Amp- ABS( Amp* 1.9182) 

1070  Volt-ABS(Volt«25) 

1080  Power-Volt*Amp 

1090  Reslstence-Volt/Amp 

1100  PRINT 

1110  BEEP 

1120  PRINT  "VOLTAGE(V)  CURRENT(A)  RESISTENCS(ahms)  POWER(W)"' 

1130  PRINT 

1140  PRINT  USING  "1X,5(MDDDD.DD, 4X)  "; Volt, Amp, Resistance, Power 

1150  WAIT  5 

1160  GOTO  370 

1170  BEEP 

1180  PRINT 

1190  PRINT  "SETUP  COMPLETE" 

1200  END 

1210  DBF  FNTvsv(V) 

1220  COM  /Cc/  C(7) 

1230  T-C(0) 

1240  FOR  I«l  TO  7 

1250  T-T+C(I)*V*I 

1260  NEXT  I 

1270  T-T+8. 626897B-2+T*( 3 . 761199B-3-T»5 . 0689259B-5 ) 

1280  RSTURH  T 

1290  FHBHD 
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APPENDIX  B.   PROGRAM  DRPJY 

This  data  acquisition/reduction  program  is  written  in 
Hewlett-Packard  Basic  5.0  language  and  is  listed  on  the 
following  pages. 


134 


10   J  FILE  NAME:  DRPJY 

20   I  DATS:       MAY  15,  1993 

30   I  REVISED  VERSION  OF  DRP72  FOR  R-124  (2  INSTRUMENTED  TUBES) 

40   t  REVISED  BY  LT.  GEORGE  PERRY  &  LT.  JOE  YUSICIAN 

50  COM  /Idp/  Idp 

60  PRINTER  IS  1 

70  CALL  Salect 

80  BEEP 

90  PRINT 

100  PRINT  "DATA  COLLECTION  COMPLETE" 

110  END 

1201 

130  SUB  Main 

140  COM  /Iprop/  Ift 

150  COM  /Idp/  Idp 

160  COM  /Cc/  C(7),Ical 

170  COM  /Wil/  D2,Di,Do,L,Lu,Kcu 

180  DIM  Emf(20),T(20),Dla(13),D2a(13),Dia(13),Doa(13),La(13),Lua(13),Kcua(13), 

Bt(19),Tn$(4)[15] 

190  DATA  0.10086091,25727.94369,-767345.8295,78025595.81 

200  DATA  -9247486589, 6. 97688E+11, -2. 66192E+13, 3. 94078E+14 

210  READ  C(*) 

220  DATA  Smooth, High  Flux, Turbo-B, High  Flux  Mod, Turbo-B  Mod 

230  READ  Tn$(") 

240  PRINTER  IS  701 

250  PRINT 

260  PRINT  "  Date  : ",DATE$(TIMEDATE) 

270  PRINT 

280  PRINT  USING  "10X, ""NOTB:  Program  name  :  DRPJY""" 

290  BEEP 

3001  INPUT  "ENTER  INPUT  MODE  (0=3054A, 1-FILE) " , Im 

310  Im-0 

3201  INPUT  "SELECT  FLUID  (0-R-114, 1=R-124) ", Ift 

330  Ift»l 

3401  INPUT  "1  OR  2  TUBE  OPBRATION  (ENTER  1  OR  2)",Hwmntu 

350  Hwmntu»2 

3601  INPUT  "SELECT  HEATING  MODE  (0=ELEC;  1=WATER) " , Ihm 

370  Ihm-0 

3801  INPUT  "BNTER  THERMOCOUPLE  TYPB  ( 0=NEW, 1-OLD) ", leal 

390  Ical-0 

400  INPUT  "GIVE  A  NAME  FOR  THE  RAW  DATA  FILE" ,D2_file$ 

410  PRINT  USING  "16X,""New  file  name:  "", 15A";D2_£ile$ 

420  Sizel=»20 

430  CREATE  BDAT  D2_Cile$, Slzel 

440  ASSIGN  dFlla2  TO  D2_file$ 

450  t 

460  t  DUMMY  FILE  UNTIL  Nrun  KNOWN 

470  Dl_fileS«"DUMMY" 

480  CREATE  BDAT  Dl.fileS, Slzel 

490  ASSIGN  QFilel  TO  Dl_file$ 

500  OUTPUT  @Fi lei; Date$ 

510  BEEP 

520  INPUT  "BNTER  NR  OP  DEFECTIVE  UPPER  TUBE  TCS  (0-DBFAULT,  2  MAX.)",Idtc 

530  IF  Idtc-0  THEN 

540  Ldtcl-0 

550  Ldtc2-0 

560  PRINT  USING  "16X,""No  defective  UPPER  TUBE  TCs  exist""" 

570  END  IF 

580  IF  Idtc-1  THEN 

590  BEEP 

600  INPUT  "BNTBR  DBFBCTIVB  UPPER  TUBB  TC  LOCATION  (TEMP  l-8)",Ldtcl 

610  PRINT  USING  "16X," "UPPER  TUBB  TC  is  defective  at  LOCATION: "",DD";Ldtcl 

620  Ldtc2»0 

630  END  IF 

640  IF  Idtc-2  THEN 

650  BBS? 

660  INPUT  "BNTER  DBT8CTIVB  UPPBR  TUBE  TC  LOCATIONS  (TEMP  1-8) ",Ldtcl,Ldtc2 

670  PRINT  USING  "16X, ""UPPER  TUBB  TC  are  defective  at  LOCATIONS: "",DD,4X,DD";L 

dtcl,Ldtc2 
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680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

ldtcl 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

11001 

11101 

11201 

11301 

11401 

11501 

1160 

11701 

1180 

1190 

1200 

1210 

1220 

1230 

12401 

12501 

12601 

1270 

12801 

1290 

13001 

13101 

13201 

13301 

1340 

1350 

1360 

1370 


END  IF 

IF  Idtc>2  THEN 
BEEP 

PRINTER  IS  1 
BEEP 

PRINT  "INVALID  ENTRY" 
PRINTER  IS  701 
GOTO  510 
END  IF 

OUTPUT  aFilel;Ldtcl,Ldtc2 
IF  Hwmntu=l  THEN  GOTO  1190 

INPUT  "ENTER  NUMBER  OF  DEFECTIVE  LOWER  TUBE  TCS  (0=DEFAULT,  2  MAX)",Aidtc 
IF  Aidtc=0  THEN 
Aldtcl=0 
Aldtc2=0 

PRINT  USING  "16X,""No  defective  LOWER  TUBE  TC3  exist""" 
END  IF 

IF  Aidtc=l  THEN 
BEEP 

INPUT  "ENTER  DEFECTIVE  TC  LOWER  TUBE  LOCATIONS  (TEMP  9-16 ) " , Aldtcl 
PRINT  USING  "16X, ""LOWER  TUBE  TC  is  defective  at  LOCATION: "" ,DD"; Aldtcl 
Aldtc2=0 
END  IF 

IF  Aidtc=2  THEN 
BEEP 

INPUT  "ENTER  DEFECTIVE  LOWER  TUBE  TC  LOCATIONS  (TEMP  9-16) ", Aldtcl ,Aldtc2 
PRINT  USING  "16X, ""LOWER  TUBE  TC  are  defective  at  LOCATIONS: "" ,DD,4X,DD", -A 
,Aldtc2 
END  IF 

IF  Aidtc>2  THEN 
BEEP 

PRINTER  IS  1 
BEEP 

PRINT  "INVALID  ENTRY" 
PRINTER  IS  701 
GOTO  790 
END  IF 

OUTPUT  @Filel; Aldtcl, Aldtc2 
PRINTER  IS  1 
BEEP 

2X, ""Select  tube  number""" 

6X,""4   Wieland  Hard  8  inch""" 
HIGH  FLUX  8  inch""" 
GEWA-K  40  Fins/in""" 
GEWA-K  26  Fins/in""" 
GEWA-T  19  Fins/in""" 
GEWA-T  OR  GEWA-TY  26  Fins/ in""" 

6X,""10  THERMOEXCEL-E""" 

6X,""11  THERMOEXCEL-HE""" 

6X,""12  TURBO-B""" 

6X,""13  GEWA-K  19  Fins/in""" 


•6X, 
•6X, 
•6X, 


"5 
""6 
""7 
""8 


PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING  "6X, 

PRINT  USING  "6X, 

PRINT  USING 

PRINT  USING 

PRINT  USING 

PRINT  USING 

INPUT  Itt 

OUTPUT  @Filel;Itt 

PRINTER  IS  701 

IF  Itt=4  THEN  PRINT  USING  " 16X, ""Tubes : 

IF  Itt=12  THEN  PRINT  USING  "16X," "Tubes 

IF  Itt=5  THEN  PRINT  USING  "16X, ""Tubes: 

IF  Itt<10  THEN  PRINT  USING  "16X,""Tube  Number 

IF  Itt>9  THEN  PRINT  USING  "16X,""Tube  Number 

INPUT  "ENTER  OUTPUT  VERSION  (  0=LONG,  l=*SHORT,  2=NONE)  "  ,  Iov 

lov=0 

INPUT  "SELECT  (0=LIQ, 1-VAP, 2=(LIQ+VAP) /2 ) ", Ilqv 

Ilqv=2 

DIMENSIONS  FOR  TESTED  TUBES 

Dl-Diameter  at  thermocouple  positions 

DATA  .0111125, .0111125, .0111125, .0129540, .012446, .0129540,-0100965 

DATA  .0100965,-01157, .01157, .01157, .01157, .01157, .0100965 

READ  Dla(*) 

Dl=Dla(Itt) 


Wieland  Hard  8  inch' 

Turbo-B""" 
High  Flux  8  inch""" 
"",D";Itt 
"",DD";Itt 
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13801 

13901  D2=Diameter  of  test  section  to  the  base  ot    fins 

DATA  .01587  5, .015875, .015875, .015824, .015875, .015824, .01270 
DATA  .0127, .0138, .0138, .0138, .0138, .0138, .0127 
READ  D2a(«) 
D2=D2a(Itt) 


Di=Inside  diameter  of  unenhanced  ends 

DATA  .0127, .0127, .0127, .0132, .0127, .0132, .0111125,-0111125 

DATA  .01 18,. 01 18,. 0118,. 01 18,. 01 18,. 0111 125 

READ  Dia(*) 

Di=Dia(Itt) 

Do-Outside  diameter  of  unenhanced  ends 

DATA  .015875, .015875, .015875, .015824, .015875, .015824, .01270, .01270 

DATA  .01331, .01331, . 0133 1 ,. 01331, . 0158, .0127 

READ  Doa(») 

Do=Doa(Itt) 


1400 

1410 

1420 

1430 

14401 

14501 

1460 

1470 

1480 

1490 

15001 

15101 

1520 

1530 

1540 

1550 

15601 

15701  L-Length  of  enhanced  surface 

1580   DATA  .1016, .1016, .1016, .1016, .2032, .2032, .2032, .2032, .2032, .2032, .2032, .20 

32, .2032, .2032 

1590   READ  La(«) 

1600   L=La(Itt) 

16101 

16201  Lu=Length  of  unenhanced  surface  at  the  ends 

1630   DATA  .0254, .0254, .0254, .0254, .0762, .0762, .0762, .0762, .0762, .0762, .0762 ,.07 

62, .0762, .0762 

1640   READ  Lua(*) 

1650  Lu=Lua(Itt) 

16601 

16701 

1680 

1690 

1700 

17101 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 


Kcu=Thermal  Conductivity  of  tube 

DATA  398,344,344,45,344,45,344,344,398,398,398,398,398,344 

READ  Kcua(*) 

Kcu=Kcua(Itt) 


Xn=.8 

Fr=.3 

Cf=3.7037E+10 

A=PI*(Do~2-Di~2)/4 

P=PI«Do 

Naub=0 

Ntc=20 

J=l 

1800  Repeat:  1 
1810  Dtld=2.22 

Ido=2 

ON  KEY  1,15  RECOVER  1800 

PRINTER  IS  1 

PRINT  USING  "2X, ""SELECT  OPTION""" 

PRINT  USING  "6X,""0=TAKE  DATA""" 

PRINT  USING  "6X,""1=SET  UPPER  TUBE  HEAT  FLUX  (MAIN)""" 
6X,""2=»SET  Tsat""" 
6X,""3=SET  LOWER  TUBE  HEAT  FLUX  (AUX)""" 


1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1920 
1930 


PRINT  USING 

PRINT  USING 

BEEP 

INPUT  I do 

IF  Ido>3  THEN  Ido-3 

IF  ldo=0  THEN  3680 


19401 

19501  LOOP  TO  SET  UPPER  TUBE  HEAT  FLUX 

1960  IF  Ido=l  THEN 

1970  OUTPUT  709; "AR  AF20  AL21  VR5" 

1980  PRINT  "SET  VOLT  BOX  TO  MAIN" 

1990  BEEP 

2000  INPUT  "ENTER  DESIRED  Qdp",Dqdp 

2010  PRINT  USING  "4X. ""DBSIRBD  Qdp  ACTUAL  Qdp" 
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2020 

Err=1000 

2030 

FOR  1=1  TO  2 

2040 

OUTPUT  709; "AS  SA" 

2050 

Sum=0 

2060 

FOR  Jl=l  TO  5 

2070 

ENTER  709; E 

2080 

Sum=Sum+E 

2090 

NEXT  Ji 

2100 

IF  1=1  THEN  Volt=Sum/5 

2110 

IF  1=2  THEN  Amp=Sum/5 

2120 

NEXT  I 

2130 

Amp= ABS ( Amp  *  1 . 9 1 82 ) 

2140 

Volt=ABS(Volt*25) 

2150 

Aqdp=Volt*Amp/(PI*D2*L) 

2160 

IF  AflS(Aqdp-Dqdp)>Err  THEN 

2170 

IF  Aqdp>Dqdp  THEN 

2180 

BEEP  4000, .2 

2190 

BEEP  4000, .2 

2200 

BEEP  4000, .2 

2210 

ELSE 

2220 

BEEP  250, .2 

2230 

BEEP  250,. 2 

2240 

BEEP  250,. 2 

2250 

END  IF 

2260 

PRINT  USING  "4X/MZ.3DE,2X,MZ.3DE";Dqdp,Aqdp 

2270 

WAIT  2 

2280 

GOTO  2030 

2290 

ELSE 

2300 

BEEP 

2310 

PRINT  USING  "4X,MZ.3DE,2X,MZ.3DE";Dqdp,Aqdp 

2320 

Err=500 

2330 

WAIT  2 

2340 

GOTO  2030 

2350 

END  IF 

2360 

END  IF 

23701 

23801 

LOOP  TO  SET  Taat 

2390 

IF  Ido=2  THEN 

2400 

IF  Ikdt=l  THEN  2440 

2410 

Dtld=2.22 

24201 

PRINT  USING  "4X,""  DTsat   ATsat    Rate    T' 

2430 

Ikdt=l 

2440 

Oldl=0 

2450 

Old2=0 

2460 

Nn=l 

2470 

Nrs=Nn  MOD  15 

2480 

Nn=Nn+l 

2490 

IF  Nra=l  THEN 

2500 

PRINT  USING  "iX,""      Taat    Tldl    Tld2 

2510 

END  IF 

2520 

OUTPUT  709; "AR  AF16  AL19  VR5" 

2530 

FOR  1=1  TO  6 

2540 

IF  I>4  THEN  2700 

2550 

Sum=0 

2560 

OUTPUT  709; "AS  SA" 

2570 

FOR  Ji=l  TO  20 

2580 

ENTER  709;Ellq 

2590 

Sum=3um+Eliq 

2600 

NEXT  Ji 

2610 

Eliq-Sum/20 

2620 

Tld=FNTvav(Bliq) 

2630 

IF  1=1  THEN  Tldl=Tld 

2640 

IF  1=2  THEN  Tld2=Tld 

2650 

IF  1=3  THEN  Tv»Tld 

2660 

IF  1=4  THEN  Taump=Tld 

2670 

IF  1=5  THEN  Tinlat=Tld 

2680 

IF  1=6  THEN  Tout=Tld 

Rata"»" 


Tv  Tsump' 
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2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

3150 

3160 

3170 

31801 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 

3270 

3280 

3290 

3300 

3310 

3320 

3330 

3340 

33.50 


NEXT  I 

Atld=(Tldl+Tld2)«.5 

IF  AflS(Atld-Dtld)>.2  THEN 

IF  Atld>Dtld  THEN 

BEEP  4000,. 2 

BEEP  4000, .2 

BEEP  4000,. 2 

ELSE 

BEEP  250, .2 

BEEP  250,  .2 

BEEP  250, .2 

END  IF 

Errl=Atld-01dl 

Oldl-Atld 

Err2=Tv-01d2 

01d2=Tv 

IF  Tldl>100.  THEN  2880 

PRINT  USING  "4X,5(MDD.DD,2X)";Dtld,Tldl,Tld2,Tv,Tsunip 

WAIT  2 

GOTO  2470 

ELSE 

IF  ABS(Atld-Dtld)>.l  THEN 

IF  Atld>Dtld  THEN 

BEEP  3000, .2 

BEEP  3000, .2 

ELSE 

BEEP  800, .2 

BEEP  800,. 2 

END  IF 

Errl=Atld-Oldl 

01dl=Atld 

Err2=Tv-01d2 

01d2=Tv 

PRINT  USING  "4X,5(MDDD.DD,2X)";Dtld,Tldl,Tld2,Tv,T3ump 

WAIT  2 

GOTO  2470 

ELSE 

BEEP 

Errl=Atld-01dl 

01dl=Atld 

Err2=Tv-01d2 

01d2=Tv 

PRINT  USING  "4X,5(MDDD.DD,2X)";Dtld,Tldl,Tld2,Tv,T3ump 

WAIT  2 

GOTO  2470 

END  IF 

END  IF 

END  IF 

I 

LOOP  TO  SET  LOWER  TUBE  HEAT  FLUX 

IF  Ido=3  THEN 

PRINT  "  SET  VOLT  BOX  TO  AUX" 

OUTPUT  709; "AR  AF20  AL22  VR5" 

BEEP 

INPUT  "ENTER  DESIRED  AuQdp" ,Duxqdp 

PRINT  USING  "2X, ""DESIRED  AuxQdp  ACTUAL  AuxQdp""" 

Err=1000 

FOR  1=1  TO  3 

OUTPUT  709; "AS  SA" 

Sum=0 

FOR  Ji=l  TO  5 

ENTER  709; E 

Sum=Sum+E 

NEXT  Ji 

IF  1-1  THEN  Volt=Sum/5 

IF  1=3  THEN  Amp=3um/5 

NEXT  I 
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3360 

3370 

3380 

3390 

3400 

3410 

3420 

3430 

3440 

3450 

3460 

3470 

3480 

3490 

3500 

3510 

3520 

3530 

3540 

3550 

3560 

3570 

3580 

3590 

3600 

36101 

3620 

3630 

3640 

3650 

3660! 

36701 

3680 

3690 

3700 

3710 

3720 

3730 

3740 

3750 

3760 

3770 

3780 

3790 

3800 

3810 

3820 

3830 

3840 

3850 

3860 

3870 

3880 

3890 

3900 

3910 

3920 

3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 


1 4X , MZ . 3DE , 2X , MZ . 3DE " ; Duxqdp , Auxqdp 


'4X,MZ.3DE,2X,MZ.3DE**;Duxqdp, Auxqdp 


Amp=ABS ( Amp* 1 . 9 182> 

Volt=ABS(Volt*25) 

Auxqdp=Volt*Amp/ ( PI *D2*L ) 

I?  ABS(Auxqdp-Duxqdp)>Err  THEN 

IF  Auxqdp>Duxqdp  THEN 

BEEP  4000, .2 

BEEP  4000,. 2 

BEEP  4000,  .2 

ELSE 

BEEP  250, .2 

BEEP  250,. 2 

BEEP  250,  .2 

END  IF 

PRINT  USING 

WAIT  2 

GOTO  3260 

ELSE 

BEEP 

PRINT  USING 

Err=500 

WAIT  2 

GOTO  3260 

END  IF 

GOTO  1850 

END  IF 

ERROR  TRAP  FOR  I do  OUT  OF  BOUNDS 

IF  Ido>3  THEN 

BEEP 

GOTO  1850 

END  IF 

TAKE  DATA  LOOP 

IF  Ikol=l  THEN  3730 

BEEP 

Ikol-1 

INPUT  "ENTER  BULK  OIL  %",Bop 

Ikol=l 

OUTPUT  709; "AR  AFOO  AL19  VR5" 

Ntc=20 

FOR  1=1  TO  Ntc 

OUTPUT  709; "AS  SA" 

Sum=0 

FOR  Ji=l  TO  20 

ENTER  709;E 

Sum=Sum+E 

IF  I=(17-Nsub)  OR  I=(18-Nsub)  THEN  Et(Ji-l)=E 

NEXT  Ji 

Kdl=0 

IF  I=(17-Nsub)  OR  I=(18-Naub)  THEN 

Eave=Sum/20 

Sum=0. 

FOR  J)c=0  TO  19 

IF  ABS(Et(Jk)-Eavo)<5.0E-6  THEN 

Sum=Sura+Et( Jk) 

ELSE 

Kdl=Kdl+l 

END  IF 

NEXT  Jk 

IF  I=(17-Nsub)  OR  I=(18-Nsub)  THEN  PRINT  USING  "4X,""Kdl  -  "",DD";Kdl 

IF  Kdl>20  THEN 

BEEP 

BEEP 

PRINT  USING  "4X,""Too  much  scattering  in  data  -  repeat  data  set""" 

GOTO  1840 

END  IF 

END  IF 

Emf(I)=»Sum/(20-Kdl) 
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4030  NEXT  I 

4040  Coun=0. 

4050  OUTPUT  709; "AR  AF20  AL22  VR5" 

4060  FOR  1=1  TO  3 

4070  OUTPUT  709; "A3  SA" 

4080  Sum=0 

4090  FOR  Ji=l  TO  5 

4100  ENTER  709;E 

4110  Sum=Sum+E 

4120  NEXT  Jl 

4130  IF  Coun=0.  THEN 

4140  IF  1=1  THEN  Vr=Sum/5 

4150  IF  1=2  THEN  Ir=Sum/5 

4160  ELSE 

4170  IF  1=1  THEN  Avr=3um/5 

4180  IF  1=3  THEN  Air=Sum/5 

4190  END  IF 

4200  NEXT  I 

4210  IF  Coun=0  THEN 

4220  PRINT  "SHIFT  VOLT  BOX  TO  AUX" 

4230  END  IF 

4240  IF  Coun=l  THEN  GOTO  4290 

4250  INPUT  "TAKE  AUX  READINGS ( 1=YES) ?" ,Ccon 

4260  IF  Ccon>l  THEN  GOTO  4250 

4270  Coun=Coun+Ccon 

4280  GOTO  4060 

42901 

43001  CONVERT  emf'S  TO  TEMP, VOLT, CURRENT 

4310  Twa=0 

4320  Atwa=0 

4330  FOR  1=1  TO  Ntc 

4340  IF  IdtoO  THEN 

4350  IF  I=Ldtcl  OR  I=Ldtc2  THEN 

4360  T(I)=-99.99 

4370  GOTO  4470 

4380  END  IF 

4390  END  IF 

4400  IF  AidtoO  THEN 

4410  IF  I=Aldtcl  OR  I=Aldtc2  THEN 

4420  T(I)=-99.99 

4430  GOTO  4470 

4440  END  IF 

4450  END  IF 

4460  T(I)=FNTvav(Emf(I)) 

4470  NEXT  I 

4480  Twa=Twa/(4-Idtc) 

4490  FOR  1=1  TO  8 

4500  IF  I=Ldtcl  OR  I=Ldtc2  THEN 

4510  Twa=Twa 

4520  ELSE 

4530  Twa=Twa+T(I) 

4540  END  IF 

4550  NEXT  I 

4560  Tw=Twa/ ( 8-Idtc ) 

4570  FOR  1=9  TO  16 

4580  IF  I=Aldtcl  OR  I=Aldtc2  THEN 

4590  Atwa=Atwa 

4600  ELSE 

4610  Atwa=Atwa+T(I) 

4620  END  IF 

4630  NEXT  I 

4640  Atw=Atwa/(8-Aidtc) 

4650  Tld=T(17-Nsub) 

4660  Tld2=T(18-Nsub) 

4670  Tlda=»{Tld+Tld2)*.5 

4680  Tv=T(19-Naub) 

4.690  Taump>T(20-N8Ub) 
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4700  Amp=ABS(Ir*1.9182) 

4710  Volt=ABS(VrJ*25 

4720  Q=Volt*Amp 

4730  Auanip=ABS(Air*1.9182) 

4740  Auvolt=ABS(Avr)*25 

4750  Auq=Auvolt*Auamp 

4760  Kcu=Kcua(Itt) 

47701 

47801  FOURIER  CONDUCTION  EQUATION  WITH  CONTACT  RESISTANCE  NEGLECTED 

4790  Tw=Tw-Q*LOG(D2/Dl)/(2*PI*Kcu*L) 

4800  Atw=Atw-Auq*LOG(D2/Dl)/(2*PI*Kcu*L) 

4810  Tsat=(Tlda+Tv)*.5 

48201 

4830!  COMPUTE  THETHAB 

4840  Thetab=Tw-Tsat 

4850  Athetab=Atw-Tsat 

4860  IF  Thetab<0  THEN 

4870  BEEP 

4880  INPUT  "UPPER  TWALL<TSAT  (0=CONTINUE,  l=END)",Iev 

4890  IF  lev=0  THEN  GOTO  3150 

4900  IF  Iev=l  THEN  1850 

4910  END  IF 

4920  IF  Athetab<0  THEN 

4930  BEEP 

4940  INPUT  "LOWER  TWA1L<TSAT  (0=CONTINUE,  1=END ) " , Aiev 

4950  IF  Alev=0  THEN  GOTO  3150 

4960  IF  Alev=l  THEN  1850 

4970  END  IF 

4980! 

4990!  COMPUTE  VARIOUS  PROPERTIES 

5000  Tfilm=(Tw+Tsat)*.5 

5010  Rho=FNRho(Tfilm) 

5020  Mu=FNMu(Tfilm) 

5030  K=FNK(Tfilm) 

5040  Cp=FNCp(Tfilm) 

5050  Beta=FNBeta(Tfilm) 

5060  Hfg=FNHfg(Tsat) 

5070  Ni=Mu/Rho 

5080  Alpha=K/(Rho*Cp) 

5090  Pr=Ni/Alpha 

5100  Psat=FNPsat(Tsat) 

5110! 

5120!  COMPUTE  VARIOUS  PROPERTIES  FOR  AUX  TUBE 

5130  Atfilm=(Atw+Tsat)*.5 

5140  Arho=FNRho(Atfllm) 

5150  Amu=FNMu(Atfilm) 

5160  Ak=FNK(Atfilm) 

5170  Acp=FNCp(Atfllm) 

5180  Abeta=FNBeta(Atfilm) 

5190  Ani=Amu/Arho 

5200  Aalpha=Ak/(Arho*Acp) 

5210  Apr=Ani/Aalpha 

5220! 

5230!  COMPUTE  NATURAL-CONVECTIVE  HEAT-TRANSFER  COEFFICIENT 

52  40!  FOR  UNENHANCED  END(S) 

5250  Hbar=190 

5260  Fe=(Hbar*P/(Kcu*A))~.5*Lu 

5270  Tanh=FNTanh(Fe) 

5280  Theta=Thetab*Tanh/Fe 

5290  Xx=(9.81*Beta*Thetab*Do~3*Tanh/(Fe*Ni*Alpha) )*. 166667 

5300  Yy»(l+(.559/Pr)*( 9/16) )~( 8/27) 

5310  Hbarc=K/Do*( .6+.387*Xx/Yy)~ 2 

5320  IF  ABS((Hbar-Hbarc)/Hbarc)>.001  THEN 

5330  Hbar=(Hbar+Hbarc)*.5 

5340  GOTO  5260 

5350  END  IF 

53601 

5370!  COMPUTE  NATURAL-CONVECTIVE  HEAT-TRANSFER  COEFFICIENT 

5380!  FOR  UNENHANCED  END(S)  FOR  AUX  TUBE 

5  390  Ahbar=190 

5400  Fe=(Ahbar*P/(Kcu*A) )~.5*Lu 

5410  Atanh=FNTanh(Fe) 
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5420  Atheta=Athetab*Atanh/Fe 

5430  Axx=(9.81*Abeta*Athetab*Do~3*Atanh/ ( Fe*Ani*Aalpha) )*. 166667 

5440  Ayy=(lt( . 559/Apr )~( 9/ 16) )~(8/27) 

5450  Ahbarc=Ak/Do*( . 6+ . 387*Axx/Ayy ) ~2 

5460  IF  ABS( ( Ahbar-Ahbarc ) /Ahbarc ) > . 001  THEN 

5470  Ahbar=(Ahbar+Ahbarc)*.5 

5480  GOTO  5400 

5490  END  IF 

55001 

5510!  COMPUTE  HEAT  LOS3  RATE  THBOUGH  UNENHANCED  END(S) 

5520  Ql=(Hbar«P«Kcu*A)~.5*Thetab«Tanh 

5530  Qc=Q-2«Ql 

5540  As=PI*D2*L 

5550  I 

5560  ICOMPUTE  HEAT  LOSS  RATE  THROUGH  UNENHANCED  END(S)  OF  AUX  TUBE 

5570  Aql=(Ahbar*P*Kcu*A)~.5*Athetab*Atanh 

5580  Aqc=Auq-2*Aql 

55901 

56001  COMPUTE  ACTUAL  HEAT  FLUX  AND  BOILING  COEFFICIENT 

5610  Qdp=Qc/Aa 

5620  Htube=Qdp/Thetab 

5630  Csf=(Cp*Thetab/Hfg)/(Qdp/(Mu«Hfg)*( .014/ ( 9. Bl*Rho)~. 5)~( 1/3. ) *Pr~1.7) 

56401 

56501  COMPUTE  ACTUAL  HEAT  FLUX  AND  BOILING  COEFFICIENT  FOR  AUX  TUBE 

5660  Auqdp-Aqc/Aa 

5670  Ahtube=Auqdp/Athetab 

5680  Acsf=(Acp*Athetab/Hfg)/(Auqdp/(Amu''Hfg)*( .014/ (9. 81*Arho)~. 5)~( 1/3. )*Apr~l 

•  7) 

56901 

57001  COMPUTE  CHURCHILL-CHU  CORRELATION  FOR  UPPER  TUBE 

5710  Cc=Hbar*Thetab 

57201 

57301  RECORD  TIME  OF  DATA  TAKING 

5740  OUTPUT  709; "TD" 

5750  ENTER  709;Told$ 

57601 

57701  OUTPUT  DATA  TO  PRINTER 

5780  PRINTER  IS  701 

5790  PRINT 

5800  PRINT  USING  "10X, ""Data  Set  Number  =  "",DDD, 2X, ""Bulk  Oil  %  =  "" ,DD.D" ; J,B 

op 

5810  PRINT  "  TIME:", TIME$ (TIMEDATE) 

5820  PRINT 

5830  PRINT  USING  "10X,""TC  No:     1       2       3       4       5       6       7 

a  ft  it  n 

5840  PRINT  USING  " 10X, ""Temp  : "" , 8( 1X,MDD.DD) ";T( 1) ,T( 2 ) ,T( 3 ) ,T( 4) ,T( 5 ) ,T( 6) ,T( 

7),T(8) 

5850  PRINT  USING  "10X,""TC  No:     9       10      11      12      13      14      15 

16""" 

5860  PRINT  USING  "10X,""Temp  : "" , 8(  1X,MDD.DD) ";T( 9) ,T( 10) ,T( 11 ) ,T( 12 ) ,T( 13 ) ,T( 1 
4),T(15),T(16) 

5870  PRINT 

5880  PRINT  USING  "10X,""   Twa   ATwa   Tllqd   Tliqd2   Tvapr   Psat   Tsat   Tau 

jjjpii  n  n 

5890  PRINT  USING  " 10X, 3 (MDD.DD, IX) , IX, MOD .DD, IX, 3( 1X,MDD .DD ) , 2X,MDD .D" ;Tw, Atw,T 

ld/Tld2,Tv,Psat/Taat,Taump 

5900  PRINT  USING  "10X,""      Vr         Ir        Cc""" 

5910  PRINT  USING  "10X, 2 (MDDD.DD, 3X) , 1X,MZ. 3DE" ; Vr, Ir , Cc 

5920  PRINT  USING  "10X,""  Thetab    Htube      Qdp      Athetab  Ahtube     AuQd 

qN  n  ti 

5930  PRINT  USING  "10X,MDD. 3D, 1X,MZ.3DE, IX, MZ. 3DE, 1X,MDD. 3D, IX, MZ. 3DE, IX, MZ. 3DE" 

; Thetab, Htube, Qdp, Athetab, Ahtube, Auqdp 

5940  PRINT 

5950  BEEP 

5960  INPUT  "OK  TO  STORE  THIS  DATA  SET  ( 1=Y, 0=N) ?" ,Ok 

5970  IF  Ok=l  THEN  J=J+1 

5980  IF  Ok=l  AND  lm=0  THEN 

5990  OUTPUT  @Filel;Bop,Told$,Emf(*),Vr,Ir,Avr,Air 

6000  END  IF 

6010  IF  Iuf»l  THEN  OUTPUT  @Ufile;VW,Uo 

6020  IF  Ire-1  THEN  OUTPUT  @Ref lle;Fms,Rew 

6030  BEEP 
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6040 

6050 

6060 

6070 

6080 

5090 

6100 

1,D2_ 

6110 

6120 

6130 

6140 

6150 

6160 

6170 

6180 

6190 

6200 

6210 

6220 

6230 

6240 

6250 

6260 

6270 

6280 

6290 

63001 

63101 

6320! 

6330 

6340! 

6350 

6360 

6370 

6380 

6390! 

6400 

6410 

6420 

6430 

5440! 

6450 

6460 

6470 

6480 

6490! 

6500 

6510 

6520 

6530 

6540 

6550! 

6560 

6570 

6580 

6590! 

5600 

6610 

6620 

6630 

6640 

6650! 

6660 

T-5 

6670 

6680 

6690 

6700 


, ZZ,""  data  runs  were  stored  in  file  "" , 15A"; J- 


INPUT  "WILL  THERE  BE  ANOTHER  RUN  ( 1=Y, 0=N ) ?" ,Go_on 
Nrun=J 

IF  Go_on=0  THEN  6080 
IF  Go_on<>0  THEN  Repeat 
BEEP 
PRINT 

PRINT  USING  "10X,""NOTE: 
file$ 
ASSIGN  @Filel  TO  * 
OUTPUT  @File2;Nrun-l 
ASSIGN  @Filel  TO  Dl_Clle$ 
ENTER  @Fi lei ; Date$ , Ldtcl , Ldtc2 , 1 tt 
ENTER  @Filel;Aldtcl,Aldtc2 
OUTPUT  @File2 ; DateS , Ldtcl , Ldtc2 , I tt 
OUTPUT  @File2;Aldtcl,Aldtc2 
FOR  1=1  TO  Nrun-1 

ENTER  @Filel;Bop,Told$,Enif  ( *  )  ,Vr,  Ir 
ENTER  @Filel;Avr,Alr 

OUTPUT  @File2;Bop,Told$,Emf(*) ,Vr,Ir 
OUTPUT  @File2;Avr,Air 
NEXT  I 

ASSIGN  @Filel  TO  * 
PRINTER  IS  1 
PURGE  "DUMMY" 
PRINT 

PRINT  "DUMMY  FILE  PURGED" 
SUBEND 

CURVE  FITS  OF  PROPERTY  FUNCTIONS 

DEF  FNKcu(T) 

OFHC  COPPER  250  TO  300  K 

Tk=T+273.15      IC  TO  K 

K=434-.112*Tlc 

RETURN  K 

FNEND 

DEF  FNMu(T) 

COM  /Iprop/  Ift 

IF  Ift=l  THEN 

!223  TO  373  K 

170  TO  360  K  CURVE  FIT  OF  VISCOUSITY 

Tk=T+273.15      IC  TO  K 

Mu=EXP(-4.4636+(1011.47/Tk) )*1.0E-3 

END  IF 

IF  Ift=l  THEN 

223  TO  373  K 

Mu=3. 4593-4. 26E-2*T+3 . 9485E-4*T~2-4 . 193E-6*T~3+2 . 0709E-8*T~4 

Mu=Mu*1.0E-4 

END  IF 

RETURN  Mu 

FNEND 

DEF  FNCp(T) 

COM  /Iprop/  Ift 

IF  lft=0  THEN 

180  TO  400  K  CURVE  FIT  OF  Cp 

Tk»T*273.15      IC  TO  K 

Cp=.40188+1.65007E-3*T)c+1.51494E-6*Tk~2-6.67853E-10*TJc~3 

Cp=Cp*1000 

END  IF 

IF  Ift=l  THEN 

223  TO  373  K  CURVE  FIT  Cp(R124) 

Cp»1.0542+2.1405E-3*T+1.0709E-5*T"2-6.4721E-8*TA3-1.4324E-9*T'*4+4.136E-ll* 

Cp=Cp*1000 
END  IF 
RETURN  Cp 
FNEND 
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67101 

6720 

6730 

6740 

6750 

6760 

6770 

6780 

6790 

6800 

68101 

6820 

6830 

6840 

6850 

68601 

6870 

6880 

6890 

69001 

6910 

6920 

69301 

6940 

6950! 

6960 

6970 

6980 

6990 

70001 

7010 

7020 

7030 

7040 

7050 

7060 

7070! 

7080 

7090 

7100 

7110 

7120 

7130 

7140 

7150 

7160 

7170 

7180 

7190 

7200 

7210! 

7220 

7230 

7240 

7250 

7260 

7270 

7280! 

7290 

7300 

7310 

7320 

7330 

7340 

7350! 

7360 

7370 


DEF  FNRho(T) 

COM  /Iprop/  Ift 

IF  lft=0  THEN 

Tk=T+273.15      !C  TO  K 

X=l-(1.8*Tk/753.95)   IK  TO  R 

Ro=36.32+61.146414*X~(l/3)+16.418015*X+17.476838*X~.5+1.119828*X~2 

Ro=Ro/. 062428 

END  IF 

IF  Ift=l  THEN 

223  TO  373  K  (R124) 

Ro=1434.8-2.8619*T-6.7267E-3*T~2-7.2852E-5*T~3 

END  IF 

RETURN  Ro 

FNEND 

DEF  FNK(T) 

COM  /Iprop/  Ift 

IF  lft=0  THEN 

T<360  K  WITH  T  IN  C 

K=.071-.000261*T 

END  IF 

RETURN  K 

IF  Ift=l  THEN 

223  TO  373  K  (R-124) 

K=7.5191E-2-3.5436E-4*T-1.9545E-7*T~2+3.1835E-9*T~3 

END  IF 

RETURN  K 

FNEND 

DEF  FNTanh(X) 

P=EXP(X) 

<2=1/P 

Tanh=(P-Q)/(P+Q) 

RETURN  Tanh 

FNEND 

DEF  FNTvsv(V) 

COM  /Cc/  C(7),Ical 

T=C(0) 

FOR  1=1  TO  7 

T=T«-C(I)*V~I 

NEXT  I 

IF  Ical=l  THEN 

T=T-6.7422934E-2+T*(9.0277043E-3-T*(-9.3259917E-5) ) 

ELSE 

T=T+8.626897E-2+T*(3.76199E-3-T*5.0689259E-5) 

END  IF 

RETURN  T 

FNEND 

DEF  FNBeta(T) 

Rop=FNRho(T+.l) 

Rom=FNRho(T-.l) 

Beta=-2 / ( Rop+Rom ) * ( Rop-Rom ) / . 2 

RETURN  Beta 

FNEND 

DEF  FNHfg(T) 

COM  /Iprop/  Ift 

IF  lft=0  THEN 

Hfg=»1.3741344E+5-T*(  3.  3094361E  +  2+T*  1.2165143) 

END  IF 

IF  Ift=l  THEN 

223  TO  373  K  CURVE  FIT  (R124) 

Hfg=159.93-.4609*T-1.458E-3*T~2-1.3715E-5*T*'3 

Hfg=Hfg*1000 
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7380  2ND  IF 

7390  RETURN  3fg 

7400  FNEND 

74101 

7420  DEF  FNPsat(Tc) 

7430  COM  /Iprop/  I£t 

7440  IF  lft=0  THEN 

7450!  0  TO  80  dag  F  CURVE  FIT  OF  Psat 

7460  Tf=1.8*Tc+32 

7470  ?a=5.945525+Tf*( . 15352082+Tf * ( I. 4840963E-3+Tf "9 . 6150671E-6) ) 

7480  ?g=Pa-14.7 

7490  END  IF 

7500  IF  ICt-1  THEN 

75101  223  TO  373  K   CURVE  FIT  (R-124) 

7520  Pa»( 162. 68+5 . 946*Tc+9. 2081S-2*Tc~2+6. 9359E-4*Tc~3 ) * . 14504 

7530  ?g=Pa-14.7 

7540  END  IF 

7550  IF  Pg>0  THEN      !  *=PSIG,-=in  Hg 

7560  Psat=Pg 

7570  ELSE 

7580  P3at=Pg*29.92/14.7 

7590  END  IF 

7600  RETURN  ?3at 

7610  FNEND 

7620  : 

7630  ! 

7640  SUB  Purg 

7650  3EEP 

7660  INPUT  "ENTER  FILE  NAME  TO  3E  DELETED " ,FiieS 

7570  PURGE  FileS 

7680  GOTO  7650 

7690  SUBEND 

7700  » 

7710  ! 

7720  SUB  Select: 

7730  COM  /Idp/  Idp 

7740  3EEP 

7750  PRINTER  IS  1 

7760  PRINT  USING  "2X, ""Select  option:""" 

7770  PRINT  USING  "6X,""  0  TAKE  DATA""" 

7780  PRINT  USING  "5X,""  1  PURGE  FILES""" 

7790  INPUT  Idp 

7800  IF  Idp-0  THEN  CALL  Main 

7810  IF  Idp-1  THEN  CALL  Purg 

7820  SUBEND 
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APPENDIX  C.   REPRESENTATIVE  DATA  SET 


Oats    :  7    Sep    1333 

MOTE:    Program    name    :    ORPJY 

Neu    Pile    nane:    O0907SC2 

UPPER   TUBE    TC    sr-   defective   at    LOCATIONS:    7  3 

LOWER    TUBE    TC    :s    defective   at    LOCATION:    9 
Tubea:      Wieiand   Hard   3    men 

Oata    Set    Numoer    =•         I       3ulk    Oil    ".    =      0.0 
TIME:  13:13:35 

TC    Mo:  I  2               3  4               5               S               7               3 

Temp    :  3.35  3.73  3.75  3.55  3.37  3.39   -99.39   -39.39 

TC    No :  9  10             1!  12              13             14             IS             IS 

Temp    :  -99.99  20.20  1 3 . 3S  20.46  19.34  20.49      13.53      17.30 

Tuia        ATua         Tliqd        Tliqd2      Tvapr      Paat         Taat         Taump 
3.79      19.36        2.20  2.22  2. II       10.32        2.15      -10.3 

Ur  Ir  Co 

.34  .17  2.388E*02 

Thetao  Htube  Odp  Athetab      Ahtube  AuQdP 

1.S2S      3.408E+02      5.542E+02      17.194      l.332E*02      3.1S0E-M33 


Oata   Set    Number   »        2      9ulk    Oil    7.   *      0.0 
TIME:  13:14:38 

TC    No:  1                 2               3               4               5               S  7  9 

Temp    :  4.37  4.31  4.92  4.32  4. 96  5.10  -99.39  -39.39 

TC    Mo:  9                10             II              12              13             14  15  IS 

Temp    :  -39.99  20.49  19.22  20.57  20.07  20.32  13.38  17.53 

Tua  ATua        Tliqd        Tliqd2      Tvapr      Paat         Taat         Taump 

4.94  19.55        2.29  2.36  2.18      10.30        2.25     -10.3 

Vr  Ir  Cc 

1.13  .23  4.530E+02 

Thetab  Htube  Qdo  Athetab      Ahtube  AuOdp 

2.538  3.3S3E-M32      I.340E+03      17.397      1 .793E^02      3.120E+03 


Oata   Set    Number   =•        3      Bulk   Oil    7.   =■      0.0 
TIME:  13:16:14 

TC   No:  I                2               3               4               5               5  7  9 

Temp    :  7.38  7.12        7.12  S.32        7.19  7.43  -39.39  -99.39 

TC    Mo:  9                10             II              12              13             14  15  16 

Temp    :  -99.99  20.08  13. 7S  20.32  19.41  20.16  18.35  16.36 

Tua        ATua         Tliqd        Tliqd2      Tvapr      Paat        Taat         Taump 
7.17      19.12        2.07  2.13  2.02      10.73        2.06      -10.3 

Ur  Ir  Cc 

1.49  .31  I.32SE+03 

Thetao  Htube  Qdp  Athetab      Ahtuoe  AuQdp 

5.104      3.410E+02      1.740E+03      17.057      I . 360E+02      3.I7ZE+03 
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Data  Set  Number  =■   4   3ulk  Oil  X  -   3.3 
TIME:      13:17:43 


TC  No 

Tamp 
TC  No 
Temp 


I  2               3  4-               5               5  7  3 

13.40  10.! a  10.12  3.54  10.01  13.53  -39.39  -39.99 

3  10             11  12             13             14  ]5  15 

•39.39  23.32  13.77  20.12  19.31  20.29  13.32  17.29 


Tuia        rtTwa        Tliqd        TIiqd2      Tvapr      Paat         Taat 


Tsump 


13. 


32   13.55    1.39   -13.7 


19.13    1.35    2.15 
\Jr  Ir  Cc 

1.37       .39     I.8S5E+03 
Thetao    Htube      QdD      nthetab   fthtube     nuOdo 
3.121   3.319E*02  2.59SE*03   17.139   I .834E+02   3.1S2E-33 


Data  Set  Number  -   5  Bulk  Oil  7.    =*  3.3 
TIME:      13:13:54 

TC  No:     1       2      3  4  5  5  7  3 

Temp  :  14.37  14.17  14.12  13.47  13.75  14.50  -39.39  -39.39 

TC  No :    9      10     11  12  13  1 4  15  IS 

Temp  :  -99.39  20.49  19.23  20.42  20.32  23.52   18.47   17.57 

Twa  rtTwa    Tliqd   Tliqd2   Tvapr  Psat    Taat    Tsump 

14.35  19.52   2.31    2.25    2.34   10.76   2.39   -13.3 
v"r        Ir       Cc 
2.3S       .53     3.346E-H33 
Thetab    Htube      Qdp      nthetab   Ahtube     AuQdp 
11.353   3.530Et02   4.342Er33   17.429   1.780E+02   3.I02E-03 


Oata  Set  Number  *   3  9ulk  Oil  %    =■  3.3 
TIME:      13:20:03 

TC  No:     1      2      3      4  5  3      7      3 

Temo  :   17.59   17.40  17.30   1G.51  15.54.   17.53  -99.39  -"99.39 


TC  No: 


0 


1 


12 


13 


IS 


IS 


Temp    :    -99.39      20.21  13.36      20.32      19.39      20.48  18.47      17.08 

Tya  ATwa         Tliqd  TL iqd2      Tvapr      Psat         Taat  Taump 

17. IS  19.27        2.34  2.32          2.39      10.30        2.13      -10.7 

Vr  Ir  Cc 

2.58  .55  4.084E+03 

Thetao  Htube  Qdo               nthetab      nhtube  nuQdo 

15.024  3.379E-H0Z     S.377E+03      17.139      1.334E+02  3.143E+03 
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Oata  Set  Numoer  =■    7   Bulk  Oil  X    -   3.3 
TIME:     13:21:13 


rc  mo 

Teno 

rc  no 

Tamo 


1      2      3      4      5  S  7  8 

19. 40  19.13  13.35  13.35  13.12  13.37  -99.39  -39.39 

9      13     II      12      13  14  IS  IS 

■39.39  13.39  13.74  19.34  13.53  23.11   13.38  17.07 


Twa   ftTua         Tliqd   Tliqd2   Tvapr   Psat    Taat    T3umo 
19.79   19.33    1.39    2.21     1.36   13.53   2.30   -13.7 
<Jr  Ir       Cc 

2.30       .50     4.7075+03 
Thetab    Htube      Qdo      Athetab   Ahtube     nuGdp 
15.733  3.5I25-32  S.36IE+33   17.327   I.348E+32   3 . 1 46E+33 


□ata  Set  Number  *   3   Sulk  Oil  %   -  0.0 
TIME:      13:22:11 

TC  No:  1  2      3      4      5      5  7  3 

Temp  :  22.75  22.21  21. 9S  20.86  21.04  22.51  -99.99  -99.99 

TC  No:  9  10      11      12      13      14  15  16 

Temp  :  -99.99  19.75  18.48  19.79  19.35  20.35  19.34  15.35 

Twa   ATua    Tliqd   Tliqd2   Tvapr   Psat    Tsat    "sump 
21.34   18.37   1.31    2.19    1.90   10.53   1.95  -10.7 
{Jr  Ir  Cc 

3.37       .55     5.3S5E-M33 
Thetab    Htube      Qdp      Athetab   Ahtube     AuQdp 
19.389  3.S52E-32  7.2S3E+03   15.925   I . 364E-32   3.155E-33 


Data  Set  Number  -   9  Bulk  Oil  X  -  0.3 
TIME:     13:23:55 

TC  No:     1  2  3  4  5      5      7  3 

Temp  :   26.31  25.51  25.32  25.12  24.95  25.71  -99.39  -"99.39 

TC  No:    9  10  II  12  13     14     15  16 

Temp  :  -99.39  20.53  19.25  20.45  23.10  20.75   13.59  17.71 

Twa   ATwa    Tliqd   TLiqd2   Tvapr  Psat    Taat    Tsump 
25.33   19.50   2.14    2.47    2.23   10.32   2.27  -10.7 
vr        Tr       Cc 
3.23       .70     7.369E+03 
Thetab    Htube      Qdp      nthetab   Ahtube     muOdp 
23.757  3.391E+02  3.0S6E+03   17.329   1.78SE+02  3.394£f03 
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Da; a    Set    Numaer    =■      !3      3ulk    Oil    X    -      3.3 

TIME:  1 3 : 27 : ZS 


TC    No 

Temp 
TC    No 

Temp 


10.92       1 ! .S3 
3  10 

-99.39       17.36 


3                4                5  5  7  3 

11.27  11.74.  11.29  13.35    -39.39    -99.39 

II              12              13  14.  IS  IS 

15.38  17.50  17.37  13.35       13.48       14.36 


Twa        ATwa         Tliqd 
I  I  .  IS       15.34         2.25 
Vr  Ir 

3.79  .31 

Thetab  Htube 


Tliqd2      Tvagr      Paat         Tsat         Tsump 
2.25  2.39       13.33        2.17      -10.3 


2.  1  UE-H33 

Qdo      mthetao   nhtube 


AuQdp 


3.391   I  .StaE+03   I.3S7E+04   14.4-71   2.4S5E*02   3.553Er03 


Oata  Sat  Number  =■   11   3uik  Oil  X  -  0.3 
TIME:     13:23:22 


TC  No 
Temp 
TC  No 

Temp 


I       2      3      4      5  S  7  3 

11.77  12.52  12.23  12.73  12.40  11.31  -39.39-99.99 

9      10     II     12     13  14  15  IS 

-99.39  IS. 30  15.33  I S . 38  16.50  17.31  16.36   13.76 


Twa  ATwa    Tliqd   Tliqd2   Tvapr   Psat    Tsat    Tsump 

12.14  IS. 21    2.23    2.22    2.31   10.73   2.12   -10. 5 
Vr        Ir       Cc 
4.48       .36     2.430E-03 
Thetab    Htube      Qdo      Athetab   nhtube     AuOdo 
13.327   1.324E-33   I.929E-H34   14.392   4.2I3E^02   5.944E+03 


Oata  Set  Nunoer  -   12   3ulk  Oil  7.  =   0.3 
TIME:      13:30:20 

TCNo:l  2  3  4  S  S  7-3 

Temp    :       13.30      14.45  14.14      14.73      14.52      13.94-99.39-99.39 
TC   No :           9                10  11  12  13  14  IS  16 

Temp    :    -99.99      15.21  14.70      16.25      15.95      17.37      15.51       13.21 

Twa   nTwa   Tliqd   Tliqd2  Tvapr  Psat   Taat   Taump 
14.37   15.53   2.11    2.14    1.33   10.70   2.33  -I0.S 
Ur        Ir       Cc 

5.38      1.25  3.372E+03 

Thetab    Htube-  Qdp      nthetao   nhtube     nuQdp 
12.039  2.792E+03  3.361E+04   13.502   2.7S5E+02  3.720E+03 
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Qata  Set  Nunoer  ■   13   3uik  Gil  % 
TIME:      13:31 : 43 


3.3 


rc  no 

rc  no 

Tamo 


!      2      3      4.      5      5      7  3 

IS. 12  IS. 75  IB.+S  17.32  15.37  IS. 17  -99.99  -99.99 

3      10     M      12     13     14     15  15 

-39.39  15. "3  '4.44  15.37  15.71  15.33   15.49  13.23 


Tua  ATwa    Tliqd 

15.24  15.34    2.21 
[Jr  Ir 

7.71       I .54 
Thetab    Htuae 

14.333  4. 1525^33  5.333£*04   13.125   2.363E+0Z   3.7SSE+03 


Tliqd2      Tvapr      Psat         T3at 

T3unio 

2.32          2.15       10.37        2.21 

-10.3 

Co 

3.739E+33 

Qao               Athetab      Ahtube 

AuQao 

Oata  Set  IMunoer  »   1 4  Bulk  Oil  X 
TIME:      13:32:43 


3.3 


TC  No 
Tenp 

rc  no 

Teno 


1       2      3      4      5      5  7  3 

17.77  13.27  17.33  13.48  13.53  17.70  -39.39  -39.39 

9      10     11     12     13     14  IS  IS 

-99.39  15.33  14.22  15.55  15.39  15. Si  14.39  12.90 


Tuia   ATwa    Tliqd   Tliqd2   Tvapr   Psat    Tsat    Tsurtp 
17.34.   15.33   2.29    2.3S    2.21   10.92   2.27   -10.5 
Vr         Ir       Cc 
3.23      I.9S     4.237E+03 
Thetaa    Htube      Qdo      Athetab   Ahtube     AuQdp 
15.375  5.4S3E+03  3.334E-04   12.753   2 . 372E-32   3.79ZE+03 


Data  Set  Number  -   15  9ulk  Oil  X    -   3.3 
TIME:      13:34:20 

TC  No:     1       2      3      4      5      5      7-3 

Temp  :   IS. 96   17.52  17.20   17.74   17.59   I S . 30  -39.39  -99.99 
TCNo:    9      10     II      12     13     14     15     IS 

Temp  :  -99.99   15.59  14.13   15.33   15.39   IS. 54   15.18   12.94 

Twa   ftTwa   Tliqd  Tliqd2   Tvapr   Psat    Tsat    Tsumo 
IS. 35   IS. 35   2.40    2.54    2.31   11.33   2.39  -13.4 
Ur        Ir       Cc 

3.23      1.75  3.32SE+03 

Thetad    Htube  Qdo      Athatab   Ahtupe     ftuQdp 
I4.5S3  4.572E-H33  S.353E+04   12.362   2.387E+02   3.792E+03 
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Oata  Sat  Number  ■   15   3ulk  Oil 
TIME:     13:35:45 


'3.3 


3  15  5  7  3 

15.38  15.59      15.32  15.53-39.39-39.39 

II      12      13     14.     15     IS 

14.3!  15.51   15.31  15.33   14.36   12.72 

TLiqd2   Tvapr   Paat    Tsat    Tsumo 
2.35    2.25   13.34   2.29  -13.4 
Cc 
3.550E+03 

Qdp      nthetab   Ahtube     nuGdp 
13.503  3.431E+03  4.334E+04   12.553   3.338E+32   3.3I3E-03 


rc  no 

I               2 

Tamp 

15.51       15.29 

7C    No 

9               10 

Temo 

-39.39      15.49 

Twa 

ATwa        TLiqd 

15.30 

14.36        2.33 

Vr 

Ir 

5.38 

1  .47 

fhete 

lb 

Htube 

Oata  Sat  Number  -   17   Bulk  Oil  5  -   0.3 
TIME:      13:36:51 

TC  No:  1       2      3      4      5      5  7  3 

Temp  :  13.75  14.46  14.18  14.31  14.21  13.55  -39.99  -99.39 

TC  No :  9      10     11      12     13     14  15  IS 

Temp  :  -99.99  15.21  13.75  15.17  15.10  16.15  14.30   12.42 

Twa  ATuta    Tliqd   TliqdZ   Tvapr   Psat    Tsat    Tsump 

14.02  14.50   2.25    2.29    2.08   10.33   2.17   -10.2 
Up         Ir       Cc 
5.28      1.13     3.010E+03 
Thetab    Htube      Qdp      Athetab   Ahtuoe     AuQop 

11.348  2.271E^03  2.590E+04   12.427   3.39SE*32   3.347E<-33 


Oata  Set  Number  =•   13   Bulk  Oil 
TIME:      13:38:09 


3.3 


TCNo:!       2  3      4      5      S      7      8 

Temp  :   11.86   12.34  12.42   12.38   12.20   11.69-99.39-39.99 

TCNo:    9      10  11     12     13     14     15     16 

Tamp  :  -99.99   15.29  13.32   15.23   15.25   16.45   15.37   12. SI 

Twa   ATwa    Tliqd  TLiqd2   Tvapr  Psat   Tsat    Tsump 

12.21   14.30   2.23  2.24    2.36   10.32   2. IS   -10.! 

'Jr                   Ir  Cc 

4.06       .37  2.437E+03 

Thetab    Htube  Qdp      Athetab  nhtube     AuQdp 
10.048   1.SS6E+03   I.5S3E+04   12.338  3.305E+02  3.797E-M33 
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Data  Set  Number  -   19   Sulk  Oil  X  ■   3.3 
TIME:      13:39:23 

TC  .No:  I       2      3  4      5      5  7  9 

7emp  :  9.20  9.53  3.-15  9.30    9.35   9.39  -39.99  -99.93 

TC  No:  3      13     M  12      13     14  15  IS 

Tamo  :  -99.99  15.53  14.21  15. 46  15.55  15.53  15.32  13. IS 

Twa  ATua    Tliqd   TliqdZ   Tvapr   Paat    T3at    Tsump 

9.50  15.13   2.23    2.35    2.17   10.90   2.24   -|0.l 
Vr        lr  Cc 

2.72       .53     I.S09E+03 
Thetab    Htube      Qdo      Athetab   Ahtube     AuQdo 

7.2S1  9.245E+02  S.712E+03   12.390   2.363E+02   3.S36E+33 


Data  Set  Number  -  20  Bulk  Oil  S  -   0.3 
TIME:      13:40:13 

TC  No:  1       2      3      4      5      3  7  3 

Temp  :  7.12  7.13  7.33   7.10   7.79   7.23  -99.39  -99.39 

TC  No:  9      10     11      12      13     14  IS  15 

Temp  :  -99.99  15.35  14.29  15.45  15.53  16.57  15.31  13.21 

Tua   frTua        Tliqd   Tliqd2   Tvapr  Paat    Tsat    Tsump 
7.21   15.15   2.25    2.32    2.22   10.31    2.2S  -10.1 
<Jr  lr  Cc 

2.34     .  .43     3.373E+02 
Thetab    Htube      Qdo      nthetab   Ahtuoe     AuQdo 
4.955   7.498Er02  3.717E+03   12.900   2.9225^02   3.7S9E-03 


Oata  Set  Number  *  21   Sulk  Oil  v.  *  0.3 
TIME:      13:41:00 

6      7    _  3 
S.33  -99.99  -19.99 

14     15     IS 
Temp  :  -99.39   15.78   14.52   15.53   15.39   16.76   15.50   13.48 

Tua        ATua   Tliqd   Tliqd2   Tvapr   Paat    Tsat    Tsump 
8.23   15.34   2.29    2.36    2.25   10.94   2.29  -10.3 
Vr        lr       Cc 
1.74       .37     7.367E-H32 
Thetab    Htube      Qdp      Athetab  Ahtube     AuQdp 
3.936  S. 301 £+02  Z.S77E+03   13.353  2.368E+02  3.744E+03 


TC    No: 

I 

2 

3 

4 

5 

Temp    : 

S.  17 

6.12 

5.34 

5.32 

5.7' 

TC    No: 

9 

10 

1  1 

12 

13 
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Data  Sat  Number  =■   22   3ulk  Gil  ". 
TIME:      13:42:13 


a. a 


TC  No 
Temp 
TC  'Mo 
Temp 


I       3      3      4      3      3  7  3 

4.33    4.40   4.39  4.3?  4.33    4.31  -39.33  -39.39 

9      13     II      12      13     14  IS  IS 

-39.39  15.38  14.70  15.33  15.37  1 S . 36  15. S3   13.55 


Tua  ATua    TLiqd   TLiqd2   Tvapr   Psat    Tsat    Tsump 

4.47  15.51    2.45    2.51    2.35   11.36   2.41   -10.1 
Vr        Ir       Cc 
1.33       .27     3.2255+02 
Thetao    Htube      Qap      Athetao   rtntube     nuQdp 
2.353  7.S40E+02   I .S72E+03   13.395   2.363E+02   3."S5E+03 


Oata  Set  Number  =■  23  Bulk  Oil  %   =■  0.3 
TIME:      13:43:07 

TC  No:  I       2      3      4      5      5  7  9 

Temp  :  3.53   3.53   3.52   3.47  3.33  3.50  -99.39  -99.99 

TC  No:  9      10     II      12      13     14  15  IS 

Temp  :  -39.39  15.38  14.54  15.38  15.98  16.39  15.54  13.42 

Tuia  ATwa    Tliqd   Tliqd2   Tvapr   Psat    Tsat    T3ump 

3.53  15.42   2.3!    2.37    2.19   10.92   2.27  -10.3 
Vr        Ir       Cc 
.99       .20     1.3265*02 
Thetao    Htube      Qdp      Athetao   Ahtube     rtuGdp 

1.317  3.523E+02   3.597E+02   13.  1 53   2.308E+32   3.S33E-03 


Data  Set  Number  -  24   Bulk  Oil  X  -  0.3 
TIME:      13:44:32 

TC  No:  1      2      3      4.      5      5  7  _  3 

Temp  :  2.35  2.34  2.95  2.33  3.12  2.39  -39.39  -99.99 

TCNo:  9      10     11      12      13     14  15     15 

Temp  :  -99.39  15.99  14.63  15.37  15.08  17.01  15.58   13.5(3 

Twa  rtTya    Tliqd   Tliqd2   Tvapr  Psat    Tsat    Tsump 

2.38  15.52   2.46    2.48    2.34   11.35   2.40  -10.3 
\lr  Ir  Cc 

.70       .14     5 . 344E+0 I 
Thetab    Htube      Qdp      Athetao   AhtuPe     AuOdo 
.574  7.3935-02   4.2415+02   13.123  2.347E+02   3.735E+03 


NOTE:  24  data  runs  jere  3tored  in  fiLe  O0907SC2 
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APPENDIX  D.   SAMPLE  CALCULATIONS3 


Data  file  D0921ASC4  was  used  for  sample  calculations.  The 
saturation  temperature  was  2.2  °C  with  a  heat  flux  of  59,130 
W/m2  on  the  smooth  upper  test  tube. 

A.  SMOOTH  TEST  TUBE  DIMENSIONS 

Do  =  0.015875  m 
Di  =  0.0127  m 
Dl  =  0.012446  m 
L  =  0.2032  m 
Lu  =  0.0762  m 

B.  MEASURED  PARAMETERS 

V  =  Vr-25  =  (7.76  V) (25)  =  194.0  V 
I  =  Ir-1.9182  =  (1.66  A)  (1.9182)  =  3.18  A 
Tl  =  19.44  °C  T2  =  20.03  °C 

T3  =  19.75  °C  T4  =  19.63  °C 

T5  =  19.46  °C  T6  =  19.38  °C 

T7  and  T8  not  used,  defective 


3This  analysis  is  essentially  identical  to  the  analysis 
done  by  Perry  [Ref.  2]  and  therefore  his  procedure  is  largely 
reproduced. 
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Tsat  =  2.20  °C 
kc  =  344  W/m-K 

C.   OUTER  WALL  TEMPERATURE  OF  THE  UPPER  TEST  TUBE 

1 .  Surface  Area  of  Smooth  Tube 

P  =  ic-Do  =  7t(0. 015875  m)  =  0.04987  m 

Ac  =  *(Do2-Di2)  /4  =  Jl((0.  015875  m) 2-  (0  .  0127  ) 2)  /4 

Ac  =  7.13  x  10"5  m2 

2 .  Heater  Power 

Qh  =  V-I  =  (194.0  V)  (3.18  A)  =  616.92  W 

3 .  Average  Inside  Tube  Temperature 

Tavg  =  Sim  Tn/m  =  V  Tn/6  =  (T1+T2+T3+T4+T5+T6) /6 
Tavg  =  (19. 44+20. 03+19. 75+19. 63+19. 46+19. 38)/6 
=  19.62  °C 

4 .  Outside  Wall  Tube  Temperature 

Two  is  calculated  using  Fourier's  Conduction  Law, 
assuming  uniform  radial  conduction  and  using  a  known  Tavg: 

Two  =  Tavg  -  {Qh  (In  (Do/Di)  )  /  (2it  -L  •  kc)  } 
Two  =  19. 62-{616.92(ln(0. 015875/0. 01245)  )/ 

(2ti-0.203-344  )  } 
Two  =  19.28  °C 

5.  Tube  Wall  Superheat 

6  =  Two  -  Tsat  =  19.28  -  2.20  =  17.08  °C 
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D.   R-124  PROPERTIES  AT  FILM  TEMPERATURE 

The  thermophysical  properties  of  R-124  were  computed  at 
the  film  temperature  using  the  REFPROP  equations  provided  by 
Bertsch  [Ref .  1]  : 

1 .  Film  Temperature 

Tf   =    (Two   +    Tsat)/2    =    (19.28    +    2.20)/2    =    10.74    °C 

2.  Dynamic  Viscosity 

u  =  3.4593-  (0.042  6)  Tf+  (3  .  94  85x1 0"4 )  Tf2-  (4  .  193xl0"6)  Tf  3+ 

(2.0709xlO-8)Tf4 
u  =    {3.  4593-(0.  0426)10.  74+ (3  .  9485xl0'4 )  (10. 74)2- 

(4.193xl0"6)  (10.74)3+(2.0709xl0-8)  ( 10  .  74  ) 4  }xl0"4 
u   =    3.0424X10'4    N-s/m2 

3 .  Density 

p  =  1434.8-(2.8619)Tf-(6.7267xl0-3)Tf2-(7.2852xl0"5)Tf3 
p   =    1434.  8-(2.  8619)10.  74- (  6  .  7267xl0'3)  (10.  74)2- 

(7.2852xl0"5)  (10. 74)3 
p  =   1403.2    kg/m3 

4.  Kinematic  Viscosity 

v    =   u/p   =   3.0422xl0"4  /1403.2   =   2.17xl0"7   m/s 

5.  Thermal  Conductivity 

k   =    7.5191xl0-2-(3.5436xl0'4)Tf-(1.9545xl0-7)Tf2  + 

(3.1835xl0"9)Tf3 
k  =   7.5191xl0-2-(3.5436xl0"4)  ( 10  .  74  )  -  ( 1 .  9545xl0~7 ) 

(10.74)2+(3.1835xl0-9)  (10.  74)3 
k   =    7.137xl0"2   W/m-K 
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6.  Specific  Heat 

Cp   =    1.0542+  (2.  1405xlO'3)Tf  +  (1.  0709xl0"5)  Tf2- 

(  6.  4 72 lxl 0"8)  -Tf3-  (1.4  324xl0"9)Tf4  +  (4  .  13  6X10"11 )  Tf5 

Cp  =    {1.0542+(2.1405xlCT3)  ( 10  .  74 )  +  ( 1 .  0709xl(T5)  (10. 74)2- 
(6.4721xl0~8)  (10.74)3-(1.4324xl0-9)  (10.74)4  + 
(4.136X10"11)  (10.74)5}xl03 

Cp   =    1078.35    J/kg-K 

7 .  Thermal   Dif f usivity 

a   =   k/(p-Cp)    =      7.137xl0'7(1403.2-1078.35) 
a   =    4.717xl0"8   m2/s 

8 .  Volumetric  Thermal  Expansion  Coefficient 

P   =    -Ap/(p-ATf)    =    -(1402.89-1403.50)/ (1403.2-0.2) 
P   =    2.174xl0"3       (1/K) 

9.  Prandtl  Number 

Pr  =   v/o   =   2.17xl0-7/4.717xl0"8   =   4.60 

E.   HEAT  FLUX  THROUGH  NON-BOILING  TUBE  LENGTHS 

Using  the  Churchill-Chu  correlation  [Ref.  9]  for  a  smooth 
cylinder  in  natural  convection,  the  average  natural  convection 
heat  transfer  coefficient  of  one  non-boiling  tube  length  is: 

h  =  k/Do{0.6+0.387[I7(l+(0.559/Pr)9/16)8/27]  }2 

where   T=    {  [g  •  P  -Do3  -0  •  tanh  (m-Lu)  ]  /  ( va  -Lu  -m)  }1/6 
m  =  {  (h-P)/  (kc«Ac)  }1/2 
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The  program  DRPJY  computes  the  value  of  h  by  assuming  a 
starting  value  of  190  W/m2-K  and  iterating  to  the  final  value. 
The  resulting  values  of  h  and  m  are: 

h  =  281.93  W/m2-K 

m  =  23.94  (1/m) 

therefore, 

Qf  =  (h-p-kc-Ac)1/2  -0-tanh(m-Lu)  =  9.52  W 

F.   HEAT  FLUX  THROUGH  BOILING  TUBE  LENGTH 

Q  =  Qh  -  2-Qf  =  616.92  -  2(9.52)  =  597.88  W 

Ab  =  n-Do-L  =  71(0.015875)  (.2032)  =  1.013xl0"2  m2 

q"  =  Q/Ab  =  597. 88/1. 013xl0"2  =  59,021  W/m2 

h  =  q"/0  =  59021/17.08  =  3456  W/m2-K 

The   following  results  were  calculated  by  the  data 
acquisition  and  reduction  program  DRPJY: 

q"    =   59,130   W/m2 

h  =   3463   W/m2-K 

0   =    17.08   °C 
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APPENDIX  E.   UNCERTAINTY  ANALYSIS4" 

The  same  data  run  used  for  the  sample  calculations 
(Appendix  D)  was  used  for  the  uncertainty  analysis,  thus  all 
the  measured  and  most  of  the  calculated  results  still  apply. 
A  lower  heat  flux  of  1561  W/m2  was  also  analyzed  from  the 
same  data  run  to  show  the  difference  between  the  two  heat 
fluxes.  All  uncertainties  are  presented  as  a  percentage  of 
the  calculated  parameter,  i.e.  a  relative  uncertainty.  The 
uncertainty  analysis  method  suggested  by  Kline  and  McClintock 

[Ref.  29]  was  used.   For  example:  If  R  =  R(x1,x2,x3, xn)  , 

then: 


bR  = 


bR 


bxx 


\2 


fix, 


iR  5*3X 


6x 


/  J 


where 


6R  =  uncertainty  of  the  desired  dependent  variable 

xn  =  measured  variable 

6xn  =  uncertainty  in  the  measured  variable 


"This  analysis  is  essentially  identical  to  the  analysis 
done  by  Perry  [Ref.  2]  and  therefore  his  procedure  is  largely 
reproduced. 
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A.   HEAT  TRANSFER  RATE  UNCERTAINTY 


Ir  =  1.66  amps 
Vr  =  7.76  volts 


6lr  =  0.025  amps 
8Vr  =  0.05  volts 


so  we   have 


hQh 
Qh 


bvi\2  +  /  5 IrV' 
Vr  j        \  Ir  ) 


bQh 
Qh 


0.025\2 


\l  .76 j         \   1.66  J 


,  i 


2      =     1.64% 


B.   SURFACE  AREA  UNCERTAINTY 

The    following   dimensions   were    taken    from   the 
manufacturer's  data  sheet,  and  from  actual  measurements: 


Do  =  15.88  (mm; 
L  =  203.2  (mm; 
Ab  =  Do  x  L 


8Do   =0.1     (mm) 
6L      =    0.1     (mm) 


Sad 

Ab 


I  bDo\2  +  /6L\2 
\  Do)         \  L  )  J 


,  i 
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6Ab 


™»    =    /     O-1    Y   W     Q-1    p      =     0.63% 
Ab         [\  15.88/         \203.2J  J 


WALL    SUPERHEAT   UNCERTAINTY 

The    following   temperatures   were  measured: 
Tl      =    19.44   °C 
T2      =   20.03   °C 
T3      =    19.75   °C 
T4      =    19.63    °C 
T5      =    19.46   °C 
T6      =    19.38    °C 
Tavg  =    19.62      °C 
Standard  Deviation  =   0.245     °C 
Tsat  =   2.20  °C      and     STsat   =   0.01   °C 


Two  =  Tavg  -  Qh 


ln(^) 

Dl 

2  K  -L-kc 


Note:  The  second  term  on  the  right  hand  side  of  the  above 
equation  is  the  Fourier  heat  transfer  conduction  term,  and  can 
be  neglected  for  the  uncertainty  analysis  since  its  value  is 
less  than  the  standard  deviation,  so: 

Two  =  Tavg  =  19.62  °C 

AT  =  Two  -  Tsat  =  19.62  °C  -  2.20  °C  =  17.42  °C 
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therefore, 


bAT  =  ["/  6A7Vo\2  +  ibTsat\2' 
AT        [\      ArJ         \      AT    ) 


SAT 


7  0.24  5\2  +  /   0.01  \2" 
\17.42J        \ 17. 42  J 


=     1.41% 


D.       HEAT   FLUX  UNCERTAINTY 


/,_     (Qh  -  2  x  Qf) 

y  Ab 


where:   Qh  =  V  x  I   =  (194  V) (3.18  A)  =  616.92  W 

and     6Qh  =  (616.92)  (0.0164)  =  10.12  W 
q"   =   59130  W/m2   (measured) 
Ab   =   0.010137  m2 
Qh  -  2.Qf  =  (q)  (Ab)  =  (59130)  (0.010137)  =  599.40  W 

solving  for  Qf: 

Qf  =  (616.92  -  599.40)/2  =  8.76  W 

assuming  the  same  proportion  uncertainty  for  Qf,  we  have: 
6Qf  =  (8.76)  (0.0164)  =  0.144  W 
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and  for  the  uncertainty: 


,g"    . 


bQh      \2  +  /      2bQf     \2  +  /6Ab\2 
<#  -  2Qf)         \Qh  -  2Qf)         \  Ab) 


,  1 

2 


bq"   m 


(1Q-12f  +  (2^°-144f  +  (0.0063)2  1 
\599.4J         \       599.4      j         V  ] 


2 


=     1.80% 


E. 


HEAT  TRANSFER  COEFFICIENT  UNCERTAINTY 


therefore, 


h  =  q'7AT 


bh 

h 


bq 


//V 


Q"  ) 


6  A  TV 


*m 


-^  =  [(0.018)2  +  (0.0141)2]2   =  2.29% 


The  above  series  of  calculations  was  also  conducted  with  a 
lower  heat  flux  for  the  same  data  run  and  with  a  high  and  low 
heat  flux  from  a  TURBO-B  data  run  (D1022FTBC4)  with  the 
results  tabulated  in  Table  5. 
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TABLE  5.   RESULTS  OF  UNCERTAINTY  ANALYSIS 


PARAMETER 

SMOOTH 
1561    W/m2 

SMOOTH 
59130   W/m2 

TURBO-B 
1703    W/m2 

TURBO-B 
71270    W/m2 

6Qh/Qh     (%) 

8.49 

1.64 

10.47 

1.60 

5Ab/Ab    (%) 

0.63 

0.63 

0.63 

0.63 

6AT/AT    (!) 

3.55 

1.41 

26.62 

12.44 

5q'7q"     (%) 

13.05 

1.80 

9.26 

1.55 

6h/h         (%) 

13.52 

2.29 

28.18 

12.54 
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